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The status of the research on muon colliders is discussed and plans are outlined for future theo- 
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retical and experimental studies. Besides continued work on the parameters of a 3-4 and 0.5 TeV 
center-of-mass (CoM) energy collider, many studies are now concentrating on a machine near 0.1 TeV 
(CoM) that could be a factory for the s-channel production of Higgs particles. We discuss the re- 
search on the various components in such muon colliders, starting from the proton accelerator 
needed to generate pions from a heavy-Z target and proceeding through the phase rotation and 
decay (n — > channel, muon cooling, acceleration, storage in a collider ring and the collider 

detector. We also present theoretical and experimental R&D plans for the next several years that 
should lead to a better understanding of the design and feasibility issues for all of the components. 
This report is an update of the progress on the R&D since the Feasibility Study of Muon Colliders 
presented at the Snowmass'96 Workshop [R. B. Palmer, A. Sessler and A. Tollestrup, Proceedings of 
the 1996 DPF/DPB Summer Study on High-Energy Physics (Stanford Linear Accelerator Center, 
Menlo Park, CA, 1997)]. 

13.10.+q,14.60.Ef,29.27.-a,29.20.Dh 
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I. INTRODUCTION 



The Standard Model of electroweak and strong interactions has passed precision experimental tests at the highest 
energy scale accessible today. Theoretical arguments indicate that new physics beyond the Standard Model associated 
with the electroweak gauge symmetry breaking and fermion mass generation will emerge in parton collisions at or 
approaching the TeV energy scale. It is likely that both hadron-hadron and lepton-antilepton colliders will be required 
to discover and make precision measurements of the new phenomena. The next big step forward in advancing the 
hadron-hadron collider energy frontier will be provided by the CERN Large Hadron Collider (LHC), a proton-proton 
collider with a center-of-mass (CoM) energy of 14 TeV which is due to come into operation in 2005. Note that in a 
high energy hadron beam, valence quarks carry momenta which are, approximately, between g to ^ of the hadron 
momentum. The LHC will therefore provide hard parton-parton collisions with typical center of mass energies of 
2.3 - 1.5 TeV. 

The route towards TeV-scale lepton-antilepton colliders is less clear. The lepton-antilepton colliders built so far 
have been e + e~ colliders, such as the Large Electron Positron collider (LEP) at CERN and the Stanford Linear 
Collider (SLC) at SLAC. In a circular ring such as LEP the energy lost per revolution in keV is 88.5 x E 4 /p, where 
the electron energy E is in GeV, and the radius of the orbit p is in meters. Hence, the energy loss grows rapidly as 
E increases. This limits the center-of-mass energy that would be achievable in a LEP-like collider. The problem can 
be avoided by building a linear machine (the SLC is partially linear), but with current technologies, such a machine 
must be very long (30-40 km) to attain the TeV energy scale. Even so, radiation during the beam-beam interaction 
(beamstrahlung) limits the precision of the CoM energy |0. 



VLHC 100 TeV pp 
(5 - 17 TeV) 




Pipetron 100 TeV pp 
(5 - 17 TeV) 



FIG. 1. Comparative sizes of various proposed high energy colliders compared with the FNAL and BNL sites. The energies 
in parentheses give for lepton colliders their CoM energies and for hadron colliders the approximate range of CoM energies 
attainable for hard parton-parton collisions. 

For a lepton with mass m the radiative energy losses are inversely proportional to m 4 . Hence, the energy- loss 
problem can be solved by using heavy leptons. In practice this means using muons, which have a mass « 207 times 



G 



that of an electron. The resulting reduction in radiative losses enables higher energies to be reached and smaller 
collider rings to be used (111. Parameters for 10 to 100 TeV collider have been discussed Estimated sizes of the 

accelerator complexes required for 0.1- TeV, 0.5-TeV and 4-TcV muon colliders are compared with the sizes of other 
possible future colliders, and with the FNAL and BNL sites in Fig. |l|. Note that muon colliders with CoM energies 
up to « 4 TeV would fit on these existing laboratory sites. The cost of building a muon collider is not yet known. 
However, since muon colliders are relatively small, they may be significantly less expensive than alternative machines. 

Since muons decay quickly, large numbers of them must be produced to operate a muon collider at high luminosity. 
Collection of muons from the decay of pions produced in proton-nucleus interactions results in a large initial phase 
volume for the muons, which must be reduced (cooled) by a factor of 10 6 for a practical collider. This may be 
compared with the antiproton stochastic cooling achieved in the Tevatron. In this case the 6-dimensional (6-D) phase 
space is reduced by approximately a factor of 10 6 , while with stacking the phase space density ||,|7j is increased by a 
factor of 10 10 . The technique of ionization cooling is proposed for the fx + pT collider p|-|ll"|]. This technique is uniquely 
applicable to muons because of their minimal interaction with matter. 

Muon colliders also offer some significant physics advantages. The small radiative losses permit very small beam- 
energy spreads to be achieved. For example, momentum spreads as low as AP/P = 0.003% are believed to be 
possible for a low-energy collider. By measuring the time-dependent decay asymmetry resulting from the naturally 
polarized muons, it has been shown [lj] that the beam energy could be determined with a precision of AE/E = 10 -6 . 
The small beam-energy spread, together with the precise energy determination, would facilitate measurements of 
the masses and widths of any new resonant states scanned by the collider. In addition, since the cross-section for 
producing a Higgs-like scalar particle in the s-channel (direct lepton-antilepton annihilation) is proportional to m , 
this extremely important process could be studied only at a muon collider and not at an e + e~ collider Finally, 
the decaying muons will produce copious quantities of neutrinos. Even short straight sections in a muon-collidcr 
ring will result in neutrino beams several orders of magnitude higher in intensity than presently available, permitting 
greatly extended studies of neutrino oscillations, nucleon structure functions, the CKM matrix, and precise indirect 
measurements of the IV-boson mass |l4| (see section II. I). 

The concept of muon colliders was introduced by G. I. Budker [§,9, and developed further by A. N. Skrinsky et 
al. |t5 -22 and D. Neuffcr jl3| , p3| B5fl . They pointed out the significant challenges in designing an accelerator complex 



that can make, accelerate, and collide /i + and pT bunches all within the muon lifetime of 2.2 /is (ct = 659 m) 
A concerted study of a muon collider design has been underway in the U.S. since 1992 [p6|-p2[. By the Sausalito 
workshop |3(if| in 1995 it was realized that with new ideas and modern technology, it may be feasible to make muon 
bunches containing a few times 10 12 muons, compress their phase space and accelerate them up to the multi-TcV 
energy scale before more than about 3/4 of them have decayed. With careful design of the collider ring and shielding 
it appears possible to reduce to acceptable levels the backgrounds within the detector that arise from the very large 
flux of electrons produced in muon decays. These realizations led to an intense activity, which resulted in the muon- 
collider feasibility study report |l3],[l4j prepared for the 1996 DPF/DPB Summer Study on High-Energy Physics (the 
Snowmass'96 workshop). Since then, the physics prospects at a muon collider have been studied extensively [|45| [47[ , 
and the potential physics program at a muon collider facility has been explored in workshops [|3S|] and conferences 

Encouraged by further progress in developing the muon-collider concept, together with the growing interest and 
involvement of the high-encrgy-physics community, the Muon Collider Collaboration became a formal entity in May 
of 1997. The collaboration is led by an executive board with members from Brookhaven National Laboratory (BNL), 
Fermi National Accelerator Laboratory (FNAL), Lawrence Berkeley National Laboratory (LBNL), Budker Institute 
for Nuclear Physics (BINP), University of California at Los Angeles (UCLA), University of Mississippi and Princeton 
University. The goal of the collaboration is to complete within a few years the R&D needed to determine whether a 
Muon Collider is technically feasible, and if it is, to design the First Muon Collider. 

Table | gives the parameters of the muon colliders under study [[l8-55|, which have CoM energies of 0.1 TeV, 



0.4 TeV and 3 TeV and Figs. || and || show possible outlines of the 0.1 TeV and 3 TeV machines. In the former 
case, parameters are given in the table for operation with three different beam-energy spreads: Ap/p = 0.12, 0.01, 
and 0.003%. In all cases, proton bunches containing 2.5-5 x 10 13 particles are accelerated to energies of 16 GeV. The 
protons interact in a target to produce O(10 13 ) charged pions of each sign. A large fraction of these pions can be 
captured in a high-field solenoid. Muons are produced by allowing the pions to decay into a lower-field solenoidal 
channel. To collect as many particles as possible within a useful energy interval, rf cavities are used to accelerate the 
lower-energy particles and decelerate the higher-energy particles (so-called phase rotation). With two proton bunches 
every accelerator cycle, the first used to make and collect positive muons and the second to make and collect negative 
muons, there are about 10 13 muons of each charge available at the end of the decay channel per accelerator cycle. 
If the proton accelerator is cycling at 15 Hz, then in an operational year (10 7 s), about 10 21 positive and negative 
muons would be produced and collected. 

As stated before, the muons exiting the decay channel populate a very diffuse phase space. The next step in the 
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TABLE I. Baseline parameters for high 
a Higgs width F = 2.7 MeV; 1 year = 10 7 


- and low-ener; 
s. 


jy muon colliders. 


Higgs/year assumes a cross section a 


= 5 x 10 4 fb; 


CoM energy 


TeV 


3 


0.4 




0.1 




p energy 


GeV 


16 


16 




16 




p's/bunch 




Z.O X 1U 


Z.O X 1U 




5 x 10 13 




Bunches/fill 




i 

4 


i 

4 




2 




Rep. rate 


Hz 


15 


15 




15 




p power 


MW 


4 


4 




4 




/i/bunch 




z X 1U 


Z X 1U 




4 x 10 12 




(i power 


MW 


28 


4 




1 




Wall power 


MW 


204 


120 




81 




Collider circum. 


m 


6000 


1000 




350 




Ave bending field 


T 


5.2 


4.7 




3 




Rms Ap/p 


% 


0.16 


0.14 


0.12 


0.01 


0.003 


6-D £6,iV 


(vrm) 3 


1. ; X 1U 


1. ( X 1U 


1.7 x 10" 10 


1.7 x 10~ 10 


1. / X 1U 


Rms t n 


7r mm-mrad 


50 


50 


85 


195 


290 


P 




0.3 


2.6 


4.1 


9.4 


14.1 




cm 


U.o 


IX) 


4.1 


9.4 


14.1 


ay spot 


/im 


3.2 


26 


86 


196 


294 


o e IP 


nirad 


1.1 


1.0 


2.1 


2.1 


2.1 


Tune shift 








0.051 


0.022 


U.U10 


"-turns (effective) 




785 


700 


450 


450 


450 


Luminosity 


cm~ 2 s _1 


7 x 10 34 


10 33 


1.2 x 10 32 


2.2 x 10 31 


10 31 


Higgs/year 








1.9 x 10 3 


4 x 10 3 


3.9 x 10 3 




muon-collider complex is to cool the muon bunch, i.e., to turn the diffuse muon cloud into a very bright bunch with 
small longitudinal and transverse dimensions, suitable for accelerating and injecting into a collider. The cooling must 
be done within a time that is short compared to the muon lifetime. Conventional cooling techniques (stochastic 
cooling [5f|_and electron cooling take too long. The technique proposed for cooling muons is called ionization 
cooling iE0,O], and will be discussed in detail in sec. V. Briefly, the muons traverse some material in which they 
lose both longitudinal and transverse momentum by ionization losses (dE/dx). The longitudinal momentum is then 
replaced using an rf accelerating cavity, and the process is repeated many times until there is a large reduction in the 
transverse phase space occupied by the muons. The energy spread within the muon beam can also be reduced by 
using a wedge-shaped absorber in a region of dispersion (where the transverse position is momentum dependent). The 
wedge is arranged so that the higher-energy particles pass through more material than lower-energy particles. Initial 
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calculations suggest that the 6-D phase space occupied by the initial muon bunches can be reduced by a factor of 
10 5 -10 6 before multiple Coulomb scattering and energy straggling limit further reduction. We reiterate that ionization 
cooling is uniquely suited to muons because of the absence of strong nuclear interactions and electromagnetic shower 
production for these particles at energies around 200 McV/c. 

Rapid acceleration to the collider beam energy is needed to avoid excessive particle loss from decay. It can be 
achieved, initially in a linear accelerator, and later in recirculating linear accelerators, rapid-cycling synchrotron, or 
fixed-field-alternating-gradient (FFAG) accelerators. Positive and negative muon bunches are then injected in opposite 
directions into a collider storage ring and brought into collision at the interaction point. The bunches circulate and 
collide for many revolutions before decay has depleted the beam intensities to an uninteresting level. Useful luminosity 
can be delivered for about 800 revolutions for the high-energy collider and 450 revolutions for the low-energy one. 

There are many interesting and challenging problems that need to be resolved before the feasibility of building 
a muon collider can be demonstrated. For example, (i) heating from the very intense proton bunches may require 
the use of of a liquid-jet target, and (ii) attaining the desired cooling factor in the ionization-cooling channel may 
require the development of rf cavities with thin beryllium windows operating at liquid-nitrogen temperatures in high 
solenoidal fields. In addition, the development of long liquid-lithium lenses may be desirable to provide stronger radial 
focusing for the final cooling stages. 

This article describes the status of our muon-collider feasibility studies, and is organized as follows. Section II gives 
a brief summary of the physics potential of muon colliders, including physics at the accelerator complex required for 
a muon collider. Section III describes the proton-driver specifications for a muon collider, and two site-dependent 
examples that have been studied in some detail. Section IV presents pion production, capture, and the pion-decay 
channel, and section V discusses the design of the ionization-cooling channel needed to produce an intense muon beam 
suitable for acceleration and injection into the final collider. Sections VI and VII describe the acceleration scenario 
and collider ring, respectively. Section VIII discusses backgrounds at the collider interaction point and section IX 
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deals with possible detector scenarios. A summary of the conclusions is given in section X. 

II. THE PHYSICS POTENTIAL OF MUON COLLIDERS 
A. Brief overview 

The physics agenda at a muon collider falls into three categories: First Muon Collider (FMC) physics at a machine 
with center-of-mass energies of 100 to 500 GcV; Next Muon Collider (NMC) physics at 3-4 TeV center-of-mass 
energies; and front-end physics with a high-intensity muon source. 

The FMC will be a unique facility for neutral Higgs boson (or techni-resonance) studies through s-channel resonance 
production. Measurements can also be made of the threshold cross sections for production of W + W~ , tt, Zh, and 
pairs of supersymmetry particles — XiXi > X2X1; and vv — that will determine the corresponding masses to high 
precision. A — > Z a factory, utilizing the partial polarization of the muons, could allow significant improvements 
in sin 2 6 W precision and in J3-mixing and CP-violating studies. In Fig. ^, we show the cross sections for SM processes 
versus the CoM energy at the FMC. For the unique s-channel Higgs boson production, where \/*uu = niH, results 
for three different beam energy resolutions are presented. 
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FIG. 4. Cross sections for SM processes versus the CoM energy at the FMC. a pt = a(fi + fi — > 7* — > e + e ). For the 
s-channel Higgs boson production, three different beam energy resolutions of 0.003%, 0.01% and 0.1% are presented. 

The NMC will be particularly valuable for reconstructing supcrsymmctric particles of high mass from their complex 
cascade decay chains. Also, any Z' resonances within the kinematic reach of the machine would give enormous event 
rates. The effects of virtual Z' states would be detectable to high mass. If no Higgs bosons exist below ~1 TeV, then 
the NMC would be the ideal machine for the study of strong WW scattering at TeV energies. 

At the front end, a high-intensity muon source will permit searches for rare muon processes sensitive to branching 
ratios that are orders of magnitude below present upper limits. Also, a high-energy muon-proton collider can be 
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constructed to probe high— Q 2 phenomena beyond the reach of the HERA ep collider. In addition, the decaying 
muons will provide high-intensity neutrino beams for precision neutrino cross-section measurements and for long- 
baseline experiments |57] Plus, there are numerous other new physics possibilities for muon facilities Q,|39| that 
we will not discuss in detail in this document. 



B. Higgs boson physics 

The expectation that there will be a light (mass below 2Mw) SM-likc Higgs boson provides a major motivation for 
the FMC, since such a Higgs boson can be produced with a very high rate directly in the s-channel. Theoretically, the 
lightest Higgs boson h° of the most general supcrsymmetric model is predicted to have mass below 150 GeV and to 
be very SM-like in the usual decoupling limit. Indeed, in the minimal supcrsymmetric model, which contains the five 
Higgs bosons h°, H°, A Q 7 H ± , one finds m h o < 130 GeV and the h° is SM-like if m A o > 130 GeV. The light SUSY 
h is regarded as the jewel in the SUSY crown. Experimentally, global analyses of precision electroweak data now 
indicate a strong preference for a light SM-like Higgs boson; this could well be the smoking gun for the SUSY Higgs 
boson. The goals of the FMC for studying the SUSY Higgs sector via s-channel resonance production are: to measure 
the light Higgs mass, width, and branching fractions with high precision, in particular sufficient to differentiate the 
MSSM h from the SM hsM', and, to find and study the heavier neutral Higgs bosons H° and A . 

The production of Higgs bosons in the s-channel with interesting rates is a unique feature of a muon collider |||7). 
The resonance cross section is 



4irT(h -> pLp) T(h -> X) 



(1) 



Gaussian beams with root-mean-square (rms) energy resolution down to R = 0.003% are realizable. The corresponding 
rms spread cr^ in CoM energy is 



100 GeV 



(2) 



The effective s-channel Higgs cross section convolved with a Gaussian spread, 



°h{Vs) = 



2ir a r- 



exp 



dV§, 



(3) 



is illustrated in Fig. | for m h = 110 GeV, T h = 2.5 MeV, and resolutions R = 0.01%, 0.06% and 0.1%. A resolution 
T/,, is needed to be sensitive to the Higgs width. The light Higgs width is predicted to be 
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(4) 



for 80 GeV < 



m h 



< 



120 GeV, where the smaller values apply in the decoupling limit of large m^o. We note that, 



in the MSSM, m^o is required to be in the decoupling regime in the context of mSUGRA boundary conditions in 
order that correct electroweak symmetry breaking arises after evolution of parameters from the unification scale. In 
particular, decoupling applies in mSUGRA at tan/3 ~ 1.8, corresponding to the infrared fixed point of the top quark 
Yukawa coupling. 

At yfs — mh-, the effective s-channel Higgs cross section is 



07, 



47T 

2 

77) — 



BF(/i -> fip,)BF(h -> X) 



-,1/2 



(5) 



BF denotes the branching fraction for h decay; also, note that <Jh ot 1/ct^- for cr^ > T 
the bb rates are 
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tot- 
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Effective Cross Sections: rn h =110 GeV 
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FIG. 5. Effective s-channel Higgs cross section obtained by convoluting the Breit-Wigner resonance formula with a 
Gaussian distribution for resolution R. From Ref. kql. 



signal ?s 10 events x L(fb ) , 
background ~ 10 events x L(fb _1 ) , 



(6) 
(7) 



assuming a 6-tagging efficiency e ~ 0.5. The effective on-resonance cross sections for other nih values and other 
channels (ZZ*, WW*) are shown in Fig. |6| for the SM Higgs. The rates for the MSSM Higgs are nearly the same as 
the SM rates in the decoupling regime of large m^o . 

The important factors that make s-channel Higgs physics studies possible at a muon collider are energy resolutions 
a r of order a few MeV, little bremsstrahlung and no beamstrahlung smearing, and precise tuning of the beam energy 

to an accuracy AE ~ 10~ 6 E through continuous spin-rotation measurements JT3] . As a case study, we consider a 
SM-likc Higgs boson with mt ~ 110 GeV. Prior Higgs discovery is assumed at the Tevatron (in Wh,tih production 
with h — > bb decay) or at the LHC (in gg — > h production with h — > 77, A£ decays with a mass measurement of 
Am?,, ~ 100 MeV for an integrated luminosity of L = 300 fb _1 ) or possibly at a NLC (in Z* — > Zh, h — > bb giving 
Am/, ~ 50 MeV for L = 200 fb _1 ). A muon collider ring design would be optimized to run at energy *Js = nih- For 
an initial Higgs-mass uncertainty of Am^ ~ 100 MeV, the maximum number of scan points required to locate the 
s-channel resonance peak at the muon collider is 



2Am h 



100 



(8) 



for a R = 0.003% resolution of <t^ w 2 MeV. The necessary luminosity per scan point (i s .p.) to observe or eliminate 

the /i-resonance at a significance level of S/y/B = 3 is L s p . ~ 1.5 x 10 _3 fb _1 . (The scan luminosity requirements 
increase for mh closer to Mz\ at nih ~ Mz the L s p . needed is a factor of 50 higher.) The total luminosity then 
needed to tunc to a Higgs boson with = 110 GeV is L tot = 0.15 fb _1 . If the machine delivers 1.5 x 10 31 cm" 2 s _1 
(0.15 fb _1 /year), then one year of running would suffice to complete the scan and measure the Higgs mass to an 
accuracy Arrih ~ 1 MeV. Figure illustrates a simulation of such a scan. 

Once the /i-mass is determined to ~ 1 MeV, a 3-point fine scan 0] can be made across the peak with higher 
luminosity, distributed with L\ at the observed peak position in ^fs and 2.5L\ at the wings (s/s = peak ± 2ct^). 

Then, with L tot = 0.4 fb _1 the following accuracies would be achievable: 16% for T% ot , 1% for crBF(bb) and 5% 
for crBF(WW*). The ratio r = BF(WW*)/BF(bb) is sensitive to m A o for m A a values below 500 GeV. For example, 
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FIG. 6. The SM Higgs cross sections and backgrounds in bb, WW* and ZZ*. Also shown is the luminosity needed for a 
5 standard deviation detection in bb. From Ref. p5[ . 

''mssmAsm = 0.3, 0.5, 0.8 for m^o = 200, 250, 400 GeV [H. Thus, using s-channel measurements of the h, it may be 
possible not only to distinguish the h° from the SM hsM but also to infer m^o. 

The study of the other neutral MSSM Higgs bosons at a muon collider via the s-channel is also of major interest. 
Finding the H° and A may not be easy at other colliders. At the LHC the region m^o > 200 GeV is deemed to 
be inaccessible for 3 < tan/3 < 5-10 p§. At an NLC the e+e~ -> H°A° production process may be kinematically 
inaccessible if H° and A are heavy (mass > 230 GeV for s/s = 500 GeV). At a 77 collider, very high luminosity 
(~200 fb _1 ) would be needed for 77 — > H°, A studies. 

At a muon collider the resolution requirements for s-channel H° and A studies are not as demanding as for the 
h° because the H°,A° widths are broader; typically V ~ 30 MeV for m A o < 2m t and r ~ 3 GeV for m^o > 2m*. 
Consequently R ~ 0.1% (<7^ ~ 70 MeV) is adequate for a scan. This is important, since higher instantaneous 

luminosities (corresponding to L ~ 2 — 10 fb _1 /yr) are possible for R ~ 0.1% (as contrasted with the L ~ 0.15 fb _ /yr 
for the much smaller R ~ 0.003% preferred for studies of the h°). A luminosity per scan point L s p . ~ 0.1 fb _1 probes 
the parameter space with tan/3 > 2. The y^-range over which the scan should be made depends on other information 
available to indicate the ^4° and H mass range of interest. A wide scan would not be necessary if r is measured with 
the above-described precision to obtain an approximate value of m^o . 

In the MSSM, m^a m H o sa m H ± at large m^o (as expected for mSUGRA boundary conditions), with a very close 
degeneracy in these masses for large tan (3. In such a circumstance, only an s-channel scan with the good resolution 
possible at a muon collider may allow separation of the A and H° states; see Fig. ||. 

C. Light particles in technicolor models 

In most technicolor models, there will be light neutral and colorless technipion resonances, 7r^ and ttj!, with masses 
below 500 GeV. Sample models include the recent top-assisted technicolor models , in which the technipion masses 
are typically above 100 GeV, and models J7(| in which the masses of the neutral colorless resonances come primarily 
from the one-loop effective potential and the lightest state typically has mass as low as 10 to 100 GeV. The widths 
of these light neutral and colorless states in the top-assisted models will be of order 0.1 to 50 GeV |71j| . In the one- 
loop models, the width of the lightest technipion is typically in the range from 3 to 50 MeV. Neutral technirho and 
tcchniomega resonances are also a typical feature of technicolor models. In all models, these resonances are predicted 
to have substantial Yukawa-like couplings to muons and would be produced in the s-channcl at a muon collider, 

AiV- -> 4, p° T , co°, (9) 
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FIG. 7. Number of events and statistical errors in the 6b final states as a function of ^/s in the vicinity of m/, SM = 110 GeV, 
assuming R = 0.003%. From Ref. H. 
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with high event rates. The peak cross sections for these processes are estimated to be « 10 4 -10 7 fb 71 1. The dominant 
decay modes depend on eigenstate composition and other details but typically are fnl 



n T ~ * 99i bb, tt, cc, tt , 



n T ~ * 99 1 bb, cc, tt, t + t~ 
Pt — > tt t itt, Witt, WW , 

U1 T — > CC, bb, TT, tt, 77Ty, Z~k\ . 



(10) 

(11) 
(12) 
(13) 



Such resonances would be easy to find and study at a muon collider. 
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D. Exotic narrow resonance possibilities 



There are important types of exotic physics that would be best probed in s-channel production of a narrow resonance 
at a muon collider. Many extended Higgs sector models contain a doubly-charged Higgs boson A~~ (and its A ++ 
partner) that couples to (jT(i~ via a Majorana coupling. The s-channel process — * A~~ has been shown |72] ] 
to probe extremely small values of this Majorana coupling, in particular values naturally expected in models where 
such couplings are responsible for neutrino mass generation. In supcrsymmetry, it is possible that there is R-parity 
violation. If R-parity violation is of the purely leptonic type, the coupling \^ T ^ for — > v T is very possibly the 
largest such coupling and could be related to neutrino mass generation. This coupling can be probed down to quite 
small values via s-channel v T production at the muon collider [[73). 



E. ^-factory 

A muon collider operating at the Z-boson resonance energy is an interesting option for measurement of polarization 
asymmetries, B®-B® mixing, and of CP violation in the -B-meson system j74|. The muon collider advantages are 
the partial muon beam polarization, and the long B-decay length for i?-mcsons produced at this y/s. The left- 
right asymmetry A^r is the most accurate measure of sin 2 w , since the uncertainty is statistics dominated. The 
present LEP and SLD polarization measurements show standard deviations of 2.4 in A n LR , 1.9 in A°p B and 1.7 in 
A°p T B |fr5|] . The CP angle (3 could be measured from B° — > K s J/i(j decays. To achieve significant improvements over 
existing measurements and those at future B- facilities, a data sample of 10 s Z-boson events/year would be needed. 
This corresponds to a luminosity > 0.15 fb _1 /year, which is well within the domain of muon collider expectations; 
R ~ 0.1% would be more than adequate, given the substantial ~ 2.4 GeV width of the Z. 



F. Threshold measurements at a muon collider 



With 10 fb 1 integrated luminosity devoted to a measurement of a threshold cross-section, the following precisions 
on particle masses may be achievable f7(|: 



W+W- AM W = 20 MeV , 
tt Am* = 0.2 GeV , 

Zh Am h = 140 MeV 



(14) 



(if m/j = 100 GeV) . Precision M\y and m* measurements allow important tests of electroweak radiative corrections 
through the relation 



AL 
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(15) 



where Sr represents loop corrections. In the SM, Sr depends on mf and logm^. The optimal precision for tests of this 
relation is AMw ~ jJo^ m t, so the uncertainty on M\y is the most critical. With AAIw = 20 MeV the SM Higgs 
mass could be inferred to an accuracy 



Am/,. 



30 GeV 



rrih 



100 GeV 



(16) 



Alternatively, once nih is known from direct measurements, SUSY loop contributions can be tested. 

In top-quark production at a muon collider above the threshold region, modest muon polarization would allow 
sensitive tests of anomalous top quark couplings |77|. 

One of the important physics opportunities for the First Muon Collider is the production of the lighter chargino, 
7§| ]. Fine-tuning arguments in mSUGRA suggest that it should be lighter than 200 GeV. A search at the 
upgraded Tevatron for the process qq — > xtx*2 with Xi ~* Xi^ +l/ an d X2 ~~ * Xi^ + ^~ decays can potentially reach 
masses m~+ ~ m^o ~ 170 GeV with 2 fb -1 luminosity and ~ 230 GeV with 10 fb -1 [[79|. The mass difference 

M(x°) — -^(Xi) can be determined from the t + t~ mass distribution. 

The two contributing diagrams in the chargino pair production process arc shown in Fig. ^; the two amplitudes 
interfere destructively. The Xi an d Vy. masses can be inferred from the shape of the cross section in the threshold 
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FIG. 9. Diagrams for production of the lighter chargino. 



region |8C|| . The chargino decay is Xi ~ * ff'Xi- Selective cuts suppress the background from W + W~ production and 
leave ~ 5% signal efficiency for 4jets + ^Jt events. Measurements at two energies in the threshold region with total 
luminosity L = 50 fb _1 and resolution R = 0.1% can give the accuracies listed in table |j] on the chargino mass for 
the specified values of m-+ and ma . 



TABLE II. Achievable uncertainties with 50 fb 1 luminosity on the mass of the lighter chargino for representative m-,+ and 
masses. From Ref. J8c| ]. 



Am-+ (MeV) 
x i 


m + (GeV) 


m c „ (GeV) 


35 


100 


500 


45 


100 


300 


150 


200 


500 


300 


200 


300 



G. Heavy particles of supersymmetry 



The requirements of gauge coupling unification can be used to predict the mean SUSY mass scale, given the value 
of the strong coupling constant at the Z-mass scale. Figure 10 shows the SUSY GUT predictions versus a s (Mz). For 
the value a s (Mz) = 0.1214 ± 0.0031 from a new global fit to precision electroweak data [[75), a mean SUSY mass of 
order 1 TeV is expected. Thus, it is likely that some SUSY particles will have masses at the TeV scale. Large masses 
for the squarks of the first family are perhaps the most likely in that this would provide a simple cure for possible 
flavor changing neutral current difficulties. 



O 



go 




0.08 0.1 0.12 0.14 

oc s (M z ) in MS scheme with DR 



0.16 



FIG. 10. a s prediction in supersymmetric GUT with minimal particle content in the Dimensional Regularization scheme. 

At the LHC, mainly squarks and gluinos will be produced; these decay to lighter SUSY particles. The LHC will 
be a great SUSY machine, but some sparticle measurements will be very difficult or impossible there | pl| , p2| , namely: 
(i) the determination of the LSP mass (LHC measurements give SUSY mass differences); (ii) study of sleptons of 
mass > 200 GeV because Drell-Yan production becomes too small at these masses; (iii) study of heavy gauginos xf 
and X3 4; which are mainly Higgsino and have small direct production rates and small branching fractions to channels 
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usable for detection; (iv) study of heavy Higgs bosons i? ± , H°, A a when the MSSM tan/3 parameter is not large and 
their masses are larger than 2mt, so that cross sections are small and decays to tt are likely to be dominant (their 
detection is deemed impossible if SUSY decays dominate). 

Detection and study of the many scalar particles predicted in supersymmetric models could be a particularly 
valuable contribution of a high energy lepton collider. However, since pair production of scalar particles at a lepton 
collider is P-wave suppressed, energies well above threshold are needed for sufficient production rates; see Fig. [□]. For 
scalar particle masses of order 1 TeV a collider energy of 3 to 4 TeV is needed to get past the threshold suppression. A 
muon collider operating in this energy range with high luminosity (L ~ 10 2 to 10 3 fb~ 1 /year) would provide sufficient 
event rates to reconstruct heavy sparticles from their complex cascade decay chains [82 84]. 



10° 



io~ 



10 



b 10' 



1 



10 



-1 



10" 



e + e or Collisions 



-''W + W ->h° 



SUGRA: tan)S=2; ,u<:0 



-'/ \Zh m = 2m 1/2 =0.5 TeV 




A 10 
b 

10 



1-1 



e^e or ju + jU Collisions 







SUGRA: tan/3 = 2; n<0 






i + Xi" 










m =2m 1/a = 0.5 TeV 






/ 

i 

1 






- 


Mr, 






_ 


l 4lAl 








1 
1 







2000 4000 6000 8000 10000 2000 4000 6000 8000 10000 

Vs (GeV) ^ (GeV) 

FIG. 11. Cross sections for pair production of Higgs bosons and scalar particles at a high-energy muon collider. From 
Ref. iH. 



In string models, it is very natural to have extra Z bosons in addition to low-energy supersymmetry. The s-channcl 
production of a Z' boson at the resonance energy would give enormous event rates at the NMC. Moreover, the s- 
channcl contributions of Z' bosons with mass far above the kinematic reach of the collider could be revealed as contact 
interactions [B5fl. 



H. Strong scattering of weak bosons 

The scattering of weak bosons can be studied at a high-energy muon collider through the process in Fig. [l^. The 
amplitude for the scattering of longitudinally polarized VF-bosons behaves like 

A{W L W L -► WlWl) ~ m 2 H /v 2 (17) 

if there is a light Higgs boson, and like 

A{W L W L -> W L W L ) ~ s ww /v 2 (18) 

if no light Higgs boson exists; here s ww is the square of the WW CoM energy and v = 246 GeV. In the latter 
scenario, partial-wave unitarity of WlWl — * WlWl requires that the scattering of weak bosons becomes strong at 
energy scales of order 1 to 2 TeV. Thus, subprocess energies ^/s ww > 1.5 TeV are needed to probe strong WW 
scattering effects. 

The nature of the dynamics in the WW sector is unknown. Models for this scattering assume heavy resonant 
particles (isospin scalar and vector) or a non-resonant amplitude based on a unitarized extrapolation of the low- 
energy theorem behavior A ~ s ww /v 2 . In all models, im pre ssive signals of strong WW scattering are obtained at 
the NMC, with cross sections typically of order 50 fb _1 |86[ . Event rates are such that the various weak-isospin 
channels (I = 0, 1, 2) could be studied in detail as a function of sww- After several years of operation, it would even 
be possible to perform such a study after projecting out the different final polarization states (WlWl, WlWt and 
WtWt), thereby enabling one to verify that it is the WlWl channel in which the strong scattering is taking place. 
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FIG. 12. Symbolic diagram for strong WW scattering. 



I. Front end physics 

New physics is likely to have important lepton flavor dependence and may be most apparent for heavier flavors. The 
intense muon source available at the front end of the muon collider will provide many opportunities for uncovering 
such physics. 

1. Rare muon decays 

The planned muon flux of ^10 14 muons/scc for a muon collider dramatically eclipses the flux, ~10 8 muons/sec, of 
present sources. With an intense source, the rare muon processes fi — > ej (for which the current branching fraction 
limit is 0.49 x 10~ 12 ), fiN — > eN conversion, and the muon electric dipole moment can be probed at very interesting 
levels. A generic prediction of supersymmetric grand unified theories is that these lepton flavor violating or CP- 
violating processes should occur via loops at significant rates, e.g. BF(// — ► ej) ~ 10~ 13 . Lepton-flavor violation can 
also occur via Z' bosons, leptoquarks, and heavy neutrinos |8jJ . 

2. Neutrino flux 

The decay of a muon beam leads to neutrino beams of well defined flavors. A muon collider would yield a neutrino 
flux 1000 times that presently available [Q. This would result in ~10 6 vN and vN events per year, which could be 
used to measure charm production (^6% of the total cross section) and measure sin 2 8 W (and infer the ly-mass to an 
accuracy AM W ~ 30-50 MeV in one year) J57|^6f|. 

3. Neutrino oscillations 

A special purpose muon ring has been proposed J58j to store ~T0 21 /i + or fjT per year and obtain ~10 20 neutrinos 
per year from muon decays along ~75-m straight sections of the ring, which would be pointed towards a distant 
neutrino detector. The neutrino fluxes from /i~ — > v^v^er or from /i + — > D^u e e + decays can be calculated with little 
systematic error. Then, for example, from the decays of stored ^~'s, the following neutrino oscillation channels could 
be studied by detection of the charged leptons from the interactions of neutrinos in the detector: 
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The detected or have the "wrong sign" from the leptons produced by the interactions of the D e and flux. 
The known neutrino fluxes from muon decays could be used for long-baseline oscillation experiments at any detector 
on Earth. The probabilities for vacuum oscillations between two neutrino flavors are given by 
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P{v a -> v b ) = sin 2 26 sin 2 (1.27 5m 2 L/E) 



(19) 



with 5m 2 in eV 2 and L/E in km/GeV. In a very long baseline experiment from Fermilab to the Gran Sasso laboratory 
or the Kamioka mine (L = O(10 4 ) km) with i/-energies E v = 20 to 50 GeV (L/E = 500-200 km/GeV), neutrino 
charged-current interaction rates of ^10 3 /year would result. In a long baseline experiment from Fermilab to the 
Soudan mine (L=732 km), the corresponding interaction rate is ~10 4 /year. Such an experiment would have sensitivity 
to oscillations down to 8m 2 ~ 10~ 4 - 10~ 5 cV 2 for sin 2 26 = 1 @. 

4- up collider 

The possibility of colliding 200-GeV muons with 1000-GeV protons at Fermilab is under study. This collider would 
reach a maximum Q 2 ~ 8 x 10 4 GeV 2 , which is ~90 times the reach of the HERA ep collider, and deliver a luminosity 
~10 33 cm -2 s _1 , which is ~300 times the HERA luminosity. The fip collider would produce ^10 6 neutral-current 
deep- inelastic-scattering events per year at Q 2 > 5000 GeV 2 , which is more than a factor of 10 3 higher than at 
HERA. In the new physics realm, lcptoquark couplings and contact interactions, if present, are likely to be larger 
for muons than for electrons. This fip collider would have sufficient sensitivity to probe leptoquarks up to a mass 
Mlq ~ 800 GeV and contact interactions to a scale A ~ 6-9 TcV @] . 

J. Summary of the physics potential 

The First Muon Collider offers unique probes of supersymmetry (particularly s-channel Higgs boson resonances) 
and technicolor models (via s-channel production of techni-resonances), high-precision threshold measurements of 
W, t and SUSY particle masses, tests of SUSY radiative corrections that indirectly probe the existence of high-mass 
squarks, and a possible Z° factory for improved precision in polarization measurements and for B-physics studies of 
CP violation and mixing. 

The Next Muon Collider guarantees access to heavy SUSY scalar particles and Z' states or to strong WW scattering 
if there are no Higgs bosons and no supersymmetry. 

The Front End of a muon collider offers dramatic improvements in sensitivity for flavor- violating transitions (e.g., 
fj, — > e~f), access to high-Q 2 phenomena in deep- inelastic muon- proton and neutrino-proton interactions, and the 
ability to probe very small Sm 2 via neutrino-oscillation studies in long-baseline experiments. 

The muon collider would be crucial to unraveling the flavor dependence of any type of new physics that is found at 
the next generation of colliders. 

Thus, muon colliders are robust options for probing new physics that may not be accessible at other colliders. 

III. PROTON DRIVER 

The overview of the required parameters is followed by a description of designs that have been studied in some 
detail. The section concludes with a discussion of the outstanding open issues. 

A. Specifications 

The proton driver requirements are determined by the design luminosity of the collider, and the efficiencies of muon 
collection, cooling, transport and acceleration. The baseline specification is for a 4-MW, 16-GeV or a 7-MW, 30-GeV 
proton driver, with a repetition rate of 15 Hz and 10 14 protons per cycle in 2 bunches (for the 100-GeV machine) or 
4 bunches (for the higher energies) of 5 x 10 13 or 2.5 x 10 13 protons, respectively. Half the bunches are used to make 
fi~ and the rest for fi + p(| . 

The total beam power is several MW, which is larger than that of existing synchrotrons. However, except for bunch 
length, these parameters are similar to those of Kaon factories |n| and spallation neutron sources p2[ . As in those 
cases, the proton driver must have very low losses to permit inexpensive maintenance of components. 

The rms bunch length for the protons on target has to be about 1 ns to: 1) reduce the initial longitudinal omittance 
of muons entering the cooling system, and 2) optimize the production of polarized muons. Although bunches of up 
to 6 x 10 13 protons per cycle have been accelerated, the required peak current is 2000 A, which is unprecedented. 



19 



Since the collection of highly polarized /i's is inefficient (see section IV. G), the proton driver should eventually 
provide an additional factor of two or more in proton intensity to permit the luminosity to be maintained for polarized 
muon beams. 



B. Possible options 

Accelerator designs are site, and to some extent, time dependent, and there have been three studies at three 
different energies (30 GeV [Q, 16 GcV 94 and 24 GeV |)5|-|)8j; see also In general, if the final energy is higher, 



the required currents are lower, bunch manipulation and apertures are easier, and the final momentum spread and 
space-charge tune shifts are less. Lowering the final energy gives somewhat more 7r's/Watt, a lower rf requirement 
(V r f ~ E 2 ) and perhaps a lower facility cost. 

In the low-energy muon collider, where two bunches of protons of 5 x 10 13 are required on target, two bunches can 
be merged outside the driver. These two bunches would be extracted simultaneously from two different extraction 
ports, and fed by different transmission lines to the same target. By arranging the path lengths of the two lines 
appropriately, the two bunches can be exactly merged. 



1. A generic design 

A 7-MW collider-driver design based on parameters originally proposed in the Snowmass Feasibility study 
consists of a 600-MeV linac, a 3.6-GcV booster and a 30-GeV driver. Both linac and booster are based on the BNL 
Spallation Neutron Source design p^ ], using a lower repetition rate and a lower total number of protons per pulse. 
For the 4-bunch case (2.5 x 10 13 protons per bunch), the (95%) bunch area is assumed to be 2 eV-s at injection and 
< 4.5 eV-s at extraction. The driver lattice is derived from the lattice of the JHF driver using 90° FODO cells with 
missing dipoles in every third FODO cell, allowing a transition energy that is higher than the maximum energy or, 
perhaps, imaginary. 



2. FNAL study 

If a muon collider is built at an existing laboratory, then possibilities abound for symbiotic relationships with the 
other facilities and programs of that laboratory. For example, the proton driver for a muon collider might result from 
an upgrade of existing proton-source capabilities, and such an upgrade could then also enhance other future programs 
that use the proton beams. 

Fermilab has conceived such a proton-source development plan |100[ with three major components: an upgraded 
linac and two rapid-cycling (15 Hz) synchrotrons: a prebooster and a new booster. The two synchrotrons operate in 
series; the four proton bunches for the muon collider are formed in the prebooster and then accelerated sequentially 
in the prebooster and the booster. The plan could be implemented in stages, and other programs would benefit from 
each stage, but all three components are required to meet the luminosity goals of the muon colliders that have been 
considered so far. 

TABLE III. Baseline proton-driver parameters of the FNAL study. 



Linac Booster Driver 

Energy range (GeV) 1 3 16 

Rep. rate (Hz) 15 15 15 

RF voltage per turn (MV) 0.15 1.5 

Circumference (m) 158 474 

Protons per bunch (xlO 13 ) 2.5 2.5 

Beam emittance [95%] (tt mm-mrad ) 200 240 

Bunch area [95%] (eV-s) 1.5 < 2.0 

Incoherent tune shift @ Inj. 0.39 0.39 



Table III presents the major parameters of the two rings. Whenever the needs of the muon collider itself allow 
some flexibility, the parameters have been chosen to optimize the resulting facility as a proton source for the rest of 
the future program at Fermilab. For example, the machine circumferences and rf-harmonic numbers result in bunch 
trains that are compatible with the existing downstream proton machines. 
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A muon collider requires proton bunches that are both very intense and, at the pion-production target, very short. 
Strong transverse and longitudinal space-charge forces might disrupt such bunches in the synchrotrons unless measures 
to alleviate those effects arc incorporated in the design. The Laslctt incoherent-space-charge tune shift quantifies the 
severity of the transverse effects. A useful approximation for the space-charge tune shift Av sc at the center of a round 
Gaussian beam is 

a 3r p Ntot , s 

2e„/?7 2 6 

In this expression r p = 1.535 x 10~ 18 m is the so-called electromagnetic radius of the proton, N to t is the total number 
of protons in the ring, e„ is the 95% normalized transverse emittance, (3 and 7 are the usual Lorentz kinematical 
factors, and b < 1 is the bunching factor, defined as the ratio of the average beam current to the peak current. 

The approximation (|2^) implies that for a given total number of protons, here 10 14 , the factors in the denominator 
are the only ways to reduce the tune shift to a specified maximum tolerable value, taken as 0.4. The bunching factor 
can be raised somewhat by careful tailoring of beam distributions, but here a typical value of 0.25 is conservatively 
assumed. Achieving the desired beam intensity then requires a combination of high injection energy, here taken as 
1 GeV into the first ring, and large transverse normalized omittances, here assumed to be about 2007T mm-mrad. 
The corresponding required aperture is about 13 cm in the first ring and about 10 cm in the second ring. With such 
large apertures in rapid-cycling synchrotrons, careful design of the beam pipes for both rings is required to manage 
eddy-current effects. Two approaches are under consideration. One is a thin Inconel pipe with water cooling and 
eddy-current coil corrections integrated on the pipe, as in the AGS Booster. The other is a ceramic beam pipe with 
a conductor inside to carry beam-image currents, as in ISIS. 

The Fermilab linac presently delivers a 40 0-Me V beam, and is capable with modest modifications of accelerating as 



many as 3 x 10 protons per cycle at 15 Hz [101]. A significant upgrade is required in order to deliver 10 protons at 



1 GeV. The energy can be raised by appending additional side-coupled modules to the downstream end of the linac. 
Increasing the linac beam intensity probably means increasing both the beam current and the duration of the beam 
pulse. Injection into the first ring is by charge stripping of the H~ beam; the incoming beam will be chopped and 
injected into pre-existing buckets to achieve high capture efficiency. 

The circumference of the second ring is set equal to that of the existing Fermilab Booster. This choice provides 
several advantages. First, the new booster could occupy the same tunnel as a relocated Booster; secondly, the beam- 
batch length from a full second ring matches that of the present Booster, which simplifies matching to downstream 
machines for other programs. The output energy of 16 GeV then results from an assumed dipole packing fraction 
of 0.575 and a peak dipole field of 1.3 T, which is the highest dipole field that is consistent with straightforward, 
nonsaturating design of magnets having thin silicon-steel laminations. Driving such magnets into saturation would 
cause significant heating of the magnet yoke as well as potential problems with tracking between the dipoles and 
quadrupoles. 

The prebooster also has 1.3-T dipole fields, and its circumference is one third that of the new booster; it operates 
at an rf harmonic number h = 4. The strategy for achieving the required short bunches at the target while alleviating 
space-charge effects in the rings is to start with four bunches occupying most of the circumference of the first small 
ring in order to keep the bunching factor large, and to do a bunch-shortening rotation in longitudinal phase space 
just before extraction from the second synchrotron. The four bunches are accelerated in the first ring to 3 GeV, then 
transferred bunch-to-bucket into the second ring with its harmonic number h = 12. At that energy, the kinematic 
factor in the tune-shift formula (^) is large enough to compensate for the smaller bunching factor in the second ring. 
The transfer energy of 3 GeV between the two rings roughly equalizes their space-charge tune shifts. 

Both rings employ separated-function lattices with flexible momentum compaction in order to raise their transition 
energies above their respective extraction energies. This not only avoids having to accelerate beam through transition 
but also provides other advantages. Intense beams are not subject to certain instabilities such as the negative- 
mass instability below transition and empirically seem less susceptible to other instabilities such as the microwave 
instability. Also, the negative natural chromaticity is beneficial for stabilizing the beam below transition, thereby 
perhaps obviating the need for sextupole correctors, especially in the first ring. Having transition not too far above 
extraction also provides substantial bucket area in which to accomplish beam-shortening rf manipulations. 



Several potential sources of instabilities in the rings have been examined [ 102 1 , as well as the possibility of com- 
pensating the latter effect by inductive inserts in the rings. Space charge is the main factor affecting the stability of 
the beams; the rings appear to be safe from longitudinal- and transverse-microwave instabilities. Of course, standard 
stabilizing methods such as active dampers are necessary to counteract some of the instabilities. Flexible momentum- 
compaction lattices would be useful not only to raise the transition energy above the extraction energy, but also to 
allow fast changes in the slip factor to facilitate bunch-narrowing manipulations at extraction time. Compensation of 
longitudinal space-charge effects by means of ferrite-loaded inductive inserts would be useful, especially for the first 
ring. 



21 



AGS-RHIC Complex 




10 14 -> 2 x 10 14 ppp + 2.5 Hz AGS 

FIG. 13. The AGS-RHIC accelerator complex and a summary of possible intensity upgrades for the AGS. 



The magnet-power-supply circuit for each ring is a 15-Hz resonant system like that of the existing Booster, with 
dipoles and quadrupoles electrically in series. This implies that the second ring will accelerate only one batch at a 
time from the first ring, which is all that the muon collider needs. Adding about 15% of second harmonic to the 
magnet ramp reduces the required peak accelerating voltage by about 25%, which is probably worth doing, especially 
for the second ring with its large voltage requirement. 

One of the advantages of a two-ring system is that the two rings divide the work of accelerating the beam. The 
rf system of the first ring is relatively modest because of its small circumference and small energy gain; that of the 
second ring is simplified because its high injection energy means a small rf- frequency swing [ 103 1 . 

ESME simulations of longitudinal motion show that the rms bunch length is 2 nsec as desired after the bunch 
rotation that occurs just before extraction from the second synchrotron. The bunch rotation creates momentum spreads 
of about 2% with longitudinal emittances of about 2 eV-s per bunch. Such spreads would contribute a few cm in 
quadrature to the beam size for a short period before extraction. This is thought to be tolerable, given the large 
apertures that are required in any case. High injection energies help to alleviate these longitudinal effects, which 
result from space-charge voltages having the same l/(3j 2 kinematic dependence as the transverse tune shifts. 



3. AGS upgrade 

The third study p5|-p8[ is of an upgrade to the BNL AGS, which should produce bunches larger than those required 
for the muon collider, but at a lower repetition rate. The AGS presently produces 6 x 10 13 protons in eight bunches 
at 25 GeV and 0.6 Hz. A 2.5-GeV accumulator ring in the AGS tunnel and AGS power-supply upgrade to 2.5-Hz 
operation would match the repetition rate to the 10-Hz repetition rate of the booster. This would generate 1 MW 
beam power. With an additional upgrade of the linac energy to 600 MeV, an intensity of 2 x 10 14 protons per pulse 
in four bunches of 5 x 10 13 at 25 GeV and 2.5 Hz could be reached, raising the power to 2 MW. The upgrades to the 
AGS accelerator complex are summarized in Fig.[H| Other options are also under consideration, such as the addition 
of a second booster and 5-Hz operation, that would reach the baseline specification of 4 MW. 

The AGS momentum acceptance of ±3% requires that the longitudinal phase space occupied by one bunch be less 
than 4.5 eV-s. This high bunch density in turn generates stringent demands on the earlier parts of the accelerator 
cycle. In particular, Landau damping from the beam momentum spread may guard against resistive wall instabilities 
during injection and longitudinal microwave instabilities after transition. Beam stability can be restored with a morc- 
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powerful transverse-damping system and possibly a new low-impedance vacuum chamber. The transverse microwave 
instability is predicted to occur after transition crossing unless damped by Landau damping from incoherent tune 
spread or possibly high-frequency quadrupoles. 



C. Progress and open issues 



Conventional rf manipulations appear able to produce 1- to 2-ns proton bunches if enough rf voltage to overcome 
the space-charge forces is used, and the beam energy is far enough from transition so the final bunch rotation is fast. 
Both simulations and experimental work have been directed at demonstrating that a short pulse can be produced 
easily. 

An experiment at the AGS has s hown tha t bunches with a z = 2 ns can be produced near transition from bunches 
with a z ~ 8 ns by bunch rotation |104, 105| . In this experiment, the AGS was flattoped near transition (~ 7 GeV) 
while the 7t-jump system was used to bring the transition energy suddenly to the beam energy letting the bunch- 
energy spread expand and bunch length contract. The experiment also demonstrated that bunches are stable over 
periods of 0.1 s. In addition, the data were used to measure the lowest two orders of the momentum compaction 
factor. 

The AGS bunch area, 1.5 eV-s, was comparable to that expected in the proton driver, but the bunch charge (though 
as large as 3-5 x 10 12 protons) was only about one tenth of that required by the muon target. The proton driver 
would use a flexible momentum compactio n la ttice which would give much better control of the transition energy, 
permitting a very fast final bunch rotation [ 106 1 . In addition, the rf frequency would be higher than that of the AGS 
so the buckets (and bunches) would initially be only half as long. Thus bunch rotation could be expected to be easier 
with the new machine, which should compensate for the larger charge. 

Simulations with the ESME code have also shown that 1-2 ns bunches of 5 x 10 13 can be produced at extraction in 



a 16-GcV ring with the rf and emittance shown in Table III 



The efficiency of capturing and accelerating beam may be increased by compensation of the space-charge forces in 
the proton driver. The use of tunable inductive inserts in the ring vacuum chamber may permit active control and 
compensation of the longitudinal space charge below transition (since the inductive impedance is of the opposite sign 
to the capacitive space charge). Initial experiments at the KEK proton synchrotron and Los Alamos PSR [ 107 1 with 
short fcrritc inserts appear to show a reduction in the synchrotron oscillation frequency shift caused by space charge 
and a decrease in the necessary rf voltage to maintain a given bunch intensity. Further experiments are needed to 
demonstrate this technique fully. 

The high rf voltage and beam power and the relatively small size of the machine require high-gradient, high-power 
rf cavities. Fcrmilab, BNL and KEK are collaborating on research and development of such type of cavities. This 
work includes the study of magnet alloys and hybrid cavities using both ferrite and new magnet alloys, high-power 
amplifiers and beam-loading compensation. 

The employment of barrier-bucket [ 108 1 rf cavities can effectively generate and manipulate a gap in the beam and 
reduce the space-charge effect. A successful test of this scheme has recently been completed ]109| |, and two 40— kV 
barrier cavities have been built by BNL and KEK and are being installed on the AGS. Another high-gradient barrier 
cavity using magnet alloys is under study at Fermilab. 



IV. PION PRODUCTION, CAPTURE AND PHASE ROTATION CHANNEL 

This section first discusses the choice of target technology and optimization of the target geometry, and then 
describes design studies for the pion capture and phase rotation channel. Prospects for polarized muon beams are 
discussed in detail. The section concludes with an outline of an R&D program for target and phase rotation issues. 

Figure |TJ gives an overview of the configuration for production of pions by a proton beam impinging on a long, 
transversely thin target, followed by capture of low-momentum, forward pions in a channel of solenoid magnets with 
rf cavities to compress the bunch energy while letting the bunch length grow. This arrangement performs the desired 
rotation of the beam. 



A. Pion production 

To achieve the luminosities for muon colliders presented in Table |, 2 x 10 12 (or 4 x 10 12 in the 100 GeV CoM case) 
muons of each sign must be delivered to the collider ring in each pulse. We estimate that a muon has a probability 
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FIG. 14. Schematic view of pion production, capture and initial phase rotation. A pulse of 16-30 GeV protons is incident on 
a skewed target inside a high-field solenoid magnet followed by a decay and phase rotation channel. 



of only 1/4 of surviving the processes of cooling and acceleration, due to losses in beam apertures or by decay. Thus, 
0.8 x 10 13 muons (1.6 xlO 13 at 100 GeV) must exit the phase rotation channel each pulse. For pulses of 2.5 x 10 13 
protons (5 x 10 13 for 100 GeV), this requires 0.3 muons per initial proton. Since the efficiency of the phase rotation 
channel is about 1/2, this is equivalent to a capture of about 0.6 pions per proton: a very high efficiency. 

The pions are produced in the interaction of the proton beam with the primary target. Extensive simulations have 
bee n perform ed for pion production from 8-30 Ge V proton beams on different target materials in a high-field solenoid 
,110 113| . Three different Monte Carlo codes [114 117j ] predict similar pion yields despite significant diff eren ces 



in their physics models. Some members of the Collaboration are involved in an AGS experiment BNL E-910 [115] to 
measure the yield of very low momentum pions, which will validate the codes in the critical kinematic region. This 
experiment ran for 14 weeks during the Spring of 1996 and has collected over 20 million events, of which about a 
quarter are minimum bias triggers for inclusive cross section measurements. The targets were varied in material (Be, 
Cu, Au, U) and thickness (2-100% interaction length (A/)) and three different beam momenta were used (6, 12.5, 
18 GeV/c). Presently, the E910 collaboration is doing a careful analysis of the large data sample obtained. Figure |l5| 
shows the dE/dx energy vs. momentum for reconstructed tracks in the TPC; there is clear particle species separation 

The pion yield is greater for relatively high Z materials, and for these, the pion yield is maximal for longitudinal 
momenta of the same order as the average transverse momentum 200 MeV/c). Targets of varying composition 
(6 < Z < 82), ra dii (0 .2-3 cm) and thicknesses (0.5-3 nuclear interaction lengths) have been explored using a Monte 
Carlo simulation [111]. For a fixed number of interaction lengths, the pion yield per proton rises almost linearly with 
proton energy, and hence almost proportional to the energy deposited in the target. The yield is higher for medium- 
and high-Z target materials, with a noticeable gain at Z > 26 for 30 GeV proton beams, but with only a minor effect 



for E < 16 GeV. This is shown in Fig. [16| where results of detailed MARS13(98) [115] simulations are presented. The 
curves show the meson yield (ir + K) from the targets in the momentum interval 0.05 < P < 0.8 GeV/c (labeled Y) 
and the number of mesons that are both captured in the high field solenoid and transported into the decay channel 
(labeled YC). The typical statistical error is a few percent. 



B. Target 

The target should be 2-3 interaction lengths long to maximize pion production. A high-density material is favored 
to minimize the size and cost of the capture solenoid magnet. Target radii larger than about 1 cm lead to lower pion 
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FIG. 15. dE/dx curve in arbitrary units for low momentum tracks; the ionization energy loss is for tracks with 30 or more 
hits in the TPC. The incident beam momentum is 18 GeV/c. From left to right the bands correspond to muons, pions, kaons, 
protons and deuterium, respectively. Note the overlap of the (nearly horizontal) electron band with other species. 



rates due to reabsorption, while smaller diameter targets have less production from secondary interactions. Tilting the 
target b y 100-150 mrad minimizes loss of pions by absorption in the target after one turn on their helical trajectory 
|30, 120 1 . Another advantage of the tilted target geometry is that the high energy and neutral components of the 
shower can be absorbed in a water-cooled beam dump to the side of the focused beam. 

About 30 kJ of energy is deposited in the target by each proton pulse (10% of the beam energy). Hence, the target 
absorbs 400 kW of power at the 15-Hz pulse rate. Cooling of the target via contact with a thermal bath would lead 
to unacceptable absorption of pions, and radiative cooling is inadequate for such high power in a compact target. 
Therefore, the target must move so as to carry the energy deposited by the proton beam to a heat exchanger outside 
the solenoid channel. 

Both moving solid metal and flowing liquid targets have been considered, with the latter as the currently preferred 
solution. A liquid is relatively easy to move, easy to cool, can be readily removed and replaced, and is the preferred 
target material for most spallation neutron sour ces u nder study. A liquid fl owing in a pipe was considered, but 
experience at ISOLDE with short proton pulses [121] as well as simulations [122,123| suggest serious problems in 
shock damage to the pipe. An open liquid jet is th us p roposed. 

A jet of liquid mercury has been demonstrated [122] but not exposed to a beam. For our application, safety and 
other considerations favor the use of a low melting point lead alloy rather than mercury. Gallium alloys, though 
with lower density, are also being considered. Experimental and theoretical studies are underway to determine the 
consequences of beam shock heating of the liquid. It is expected that the jet will disperse after being exposed to the 
beam. The target station must survive damage resulting from the violence in this dispersion. This consideration will 
determine the minimum beam, and thus jet, radius. 
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FIG. 16. Meson yield (tt + K) from different targets tilted by angle a in a solenoidal field B of aperture R a as calculated with 
the MARS13(98) code. The target is aligned along the beam. The curves labeled YC show mesons that are transported into the 
decay channel, (a) Yield from a 1.5 Aj, 1 cm radius target irradiated with 8, 16 and 30 GeV proton beams (a x = a y = 4 mm) 
as a function of target atomic mass (B = 20 T, R a = 7.5 cm, a = 0); (b) Yield from a 3 — Xi, 1 cm radius gallium target 
tilted at a = 150 mrad in a 16 GeV proton beam (a x — a y = 4 mm) vs. solenoid field for a fixed adiabatic invariant BR^; (c) 
Yield as a function of radius of a 3 Xi gallium target in a 16 GeV proton beam (<j x — a y = 4 mm, B = 20 T, R a = 7.5 cm, 
a = 100 mrad); (d) Yield from a 3 A/, 1 cm radius gallium target vs. tilt angle between the axis of the capture solenoid and 
the proton beam for a 16 GeV proton beam (a x — a y — 4 mm, B = 20 T, R a — 7.5 cm). 



For a conducting liquid jet in a strong magnetic field, as proposed, strong eddy currents will be induced in the jet, 
causing reaction forces that may disrupt its flow [124,125]. The forces induced are proportional to the square of the 
jet radius, and set a maximum for this radius of order 5-10 mm. If this maximum is smaller than the minimum radius 
set by shock considerations, then multiple smaller beams and jets could be used; e.g., four jets of 5 mm radius with 
four beams with 2.5 x 10 13 protons per bunch. Other alternatives inclu de targets made from insulating materials such 
as liquid PtC>2 or Re203, slurries (e.g., Pt in water), or powders [126]. 

A moving solid metal target is not the current baseline solution, but is a serious possibility. In this case, the target 
could consist of a long fiat band or hoop of copper- nickel that moves along its length (as in a band saw) 1 127 ] . The 
band would be many meters in length, would be cooled by gas jets away from the target area, and would be supported 
and moved by rollers, as shown in Fig. [l7]. 

The choice and parameters of the target are critical issues that need resolution. These can be resolved by experiments 
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FIG. 17. Alternative concept of a solid metal target in the form of a rotating Cu-Ni band. 



with a strong magnetic field and a beam, as discussed in section IV. H. 



C. Capture 



To capture all pions with transverse momenta px less than their typical value of 200 MeV/c, the product of the 
capture solenoid field B and its radius R a must be greater than 1.33 T m. The use of a high field and small radius 
is preferred to minimize the corresponding transverse emittance, which is proportional to BR 2 : for a fixed transverse 
momentum capture, this emittance is thus proportional to R. A field of 20 T and 7.5 cm radius was chosen on the 
basis of simulations described below. This gives BR = 1.5 T m, BR 2 = 0.1125 T m 2 and a maximum transverse 
momentum capture of p x — 225 McV/c. 

A preliminary design |128| of the capture solenoid has an inner 6 T, 4 MW, water cooled, hollow conductor magnet 
with an inside diameter of 24 cm and an outside diameter of 60 cm. There is space for a 4 cm thick, water cooled, 
heavy metal shield inside the coil. The outer superconducting magnet has three coils, with inside diameters of 60 to 
80 cm. It generates an additional 14 T of field at the target and provides the required tapered fiel d to match into the 
decay channel. Such a hybrid solenoid has parameters compatibl e wi th those of existing magnets 12£j. 

The 20 T capture solenoid is matched via a transfer solenoid [110] into a decay channel consisting of a system of 
superconducting solenoids with the same adiabatic invariant BR 2 oc Rpr- Thus, for a 1.25 T decay channel, B drops 
by a factor of 16 between the target and decay channel, R and pr change by factors of 4 and 1/4, respectively. This 
permits improved acceptance of transverse momentum within the decay channel, at the cost of an increased spread 
in longitudinal momentum. Figure ^(b) shows the meson yield as a function of field in the capture solenoid, with 
the radius of the capture solenoid adjusted to maintain the same BR 2 as in the decay channel. The optimum field is 
20 T in the capture solenoid. 

If the axis of the target is coincident with that of the solenoid field, then there is a relatively high probability that 
pions re-enter the target after one cycle on their helical trajectory and are lost due to nuclear int eractions. When 
the target and proton beam are set at an angle of 100-150 mrad with respect to the field axis [111], the probability 
for such pion interactions at the target is reduced, and the overall production rate is increased by 60%, as shown in 
Fig. ||(d). 

In summary, the simulations indicate that a 20 T solenoid of 16 cm inside diameter surrounding a tilted target will 
capture about half of all produced pions. With target efficiency included, about 0.6 pions per proton will enter the 
pion decay channel [111]. 
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D. Phase rotation linac 



The pions, and the muons into which they decay, have a momentum distribution with an rms spread of approximately 
100% and a peak at about 200 MeV/c. It would be difficult to handle such a wide spread in any subsequent system. 
A linac is thus introduced along the decay channel, with frequencies and phases chosen to decelerate the fast particles 
and accelerate the slow ones; i.e., to phase rotate the muon bunch. Several studies have been made of the design of 
this system, using differing ranges of rf frequency, delivering different final muon momenta, and differing final bunch 
lengths. In all cases, muon capture efficiencies of close to 0.3 muons per proton are obtained. Until the early stages 
of the ionization cooling have been designed, it is not yet possible to choose between them. Independent of the above 
choices is a question of the location of the focusing solenoid coils and rf cavity design, as discussed below in the 
section IV.F. 



TABLE IV. Parameters of the Lower-Energy Phase Rotation Linacs 



Linac 


Length 


Frequency 


Gradient 




(m) 


(MHz) ' 


(MeV/m) 


1 


3 


60 


5 


2 


29 


30 


4 


3 


5 


60 


4 


4 


5 


.37 


4 



1. Lower energy, longer bunch examp 



This example captures muons at a mean kinetic energy of 130 McV. Table ^ gives parameters of the linacs used 



Monte Ca rlo simulations ]52|, with the program MUONMC [13C], were done using pion production calculated by 



ARC [114| for a copper target of 1-cm radius at an angle of 150 mrad. A uniform solenoidal field was assumed in the 
phase rotation, and the rf was approximated by a series of kicks. 

Figure [l^ shows the energy vs. ct at the end of the decay and phase rotation channel. The abscissa ct is a measure 
of bunch length at the end of the channel: the total transit time of each n/fj, is multiplied by the velocity of light and 
the total length of the channel is subtracted. Thus a ficticious reference particle at the center of the incident bunch 
at the target arrives at ct = m. A loose final bunch selection was defined with an energy 130 ± 70 MeV and bunch 
ct from 3 to 11 m. With this selection, the rms energy spread is 16.5%, the rms ct is 1.7 m, and there are 0.39 muons 
per incident proton. A tighter selection with an energy 130 ± 35 MeV and bunch ct from 4 to 10 m gave an rms 
energy spread of 11.7%, rms ct of 1.3 m, and contained 0.31 muons per incident proton. 



2. Higher energy, shorter bunch example 
In this example the captured muons have a mean kinetic e nergy close to 320 McV. It is based on a Monte Carlo 



study which uses the updated MARS pion production model |11£] to generate pions created by 16 GeV protons on 
a 36 cm long, 1 cm radius coaxial gallium target. Figure [l9| shows the longitudinal phase space of the muons at the 
end of an 80 m long, 5 T solenoidal decay channel with cavities of frequency in the 30-90 MHz range and acceleration 
gradients of 4-18 MeV/m. A total of 0.33 muons per proton fall within the indicated cut (6 mx300 MeV). The 
rms bunch length inside the cut is 148 cm and rms energy spread is 62 MeV. The normalized six dimensional (6-D) 
cmittancc is 217 cm 3 and the transverse part is 1.86 cm (the normalized 6-D omittance eg. at is defined in section V). 

A sample simulation with lithium hydride absorbers regularly spaced in the last 60 m of a 120 m decay channel and 
with compensating acceleration captures 0.3 muons with mean kinetic energy of about 380 MeV in a (6 mx300 McV) 
window. The longitudinal phase space is about the same as in the previous example but the transverse part shrinks 
to 0.95 cm due to ionization cooling which reduces the 6-D phase space to 73.5 cm 3 . 
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FIG. 18. Energy vs. ct of /x's at end of the lower-energy phase rotation channel. The symbols +, o and 
polarization P > |, — i < P < | and P < — |, respectively. 



denote muons with 



E. Use of both signs 



Protons on the target produce pions of both signs, and a solenoid will capture both, but the subsequent rf systems 
will have opposite effects on each sign. The proposed baseline approach uses two separate proton bunches to create 
separate positive and negative pion bunches and accepts the loss of half the pions/muons during phase rotation. 

If the pions can be charge separated with limited loss before the phase rotation cavities are reached, then h ighe r 
luminosity may be obtained. The separation of charged pions in a curved solenoid decay line was studied in [ 1 10 1 . 
Because of the resulting dispersion in a bent solenoid, an initial beam of radius R with maximum-to-minimum 
momentum ratio F will require a large beampipe of radius (1 + F)R downstream to accommodate the separated 
beams. A septum can then be used to capture the two beams into separate channels. Typically the reduction in 
yield for a curved solenoid compared to a straight solenoid is about 25% (due to the loss of very low and very high 
momentum pions to the walls or septum) , but this must be weighed against the fact that both charge signs are captured 
for an overall net gain. A disadvantage is that this charge separation takes place over several meters of length during 
which time the beam spreads longitudinally. This makes capture in an rf phase rotation system difficult, although 
a large aperture cavity system could be incorporated in the bent solenoid region to alleviate this. The technique 
deserves further study and may be useful to consider as an intensity upgrade to a muon collection system. 



F. Solenoids and rf 



As noted above, capture using higher frequencies appears to be less efficient, and most studies now use frequencies 
down t o 30 MHz. Such cavities, when conventionally designed, are very large (about 6.6 m diameter). In the Snowmass 
study [131 1 a reentrant design reduced this diameter to 2.52 m, but this is still large, and it was first assumed that 
the 5 T focusing solenoids would, for economic reasons, be placed within the irises of the cavities (see Fig. p0| ). 

A study of transmission down a realistic system of iris located coils revealed betatron resonant excitation from the 
magneti c fie ld periodicities, leading to significant particle loss. This was reduced by the use of more complicated coil 
shapes [ 131 [ , smaller gaps, and shorter cavities, but remained a problem. 

An alternative is to place continuous focusing coils outside the cavities as shown in Fig. In this case, cost will be 
minimized with lower magnetic fields (1.25-2.5 T) and correspondingly larger decay channel radii (21-30 cm). Studies 
are underway to determine the optimal solution. 
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FIG. 19. Longitudinal phase space at the end of decay channel with projections onto time and energy axes. The four dashed 
lines delineate the region deemed acceptable for the cooling channel. 




FIG. 20. Schematic of capture and phase rotation using rf cavities with superconducting solenoids (hatched) inside the irises. 
Three groups of three cavities operating at 90, 50, and 30 MHz are shown from left to right, respectively. 



G. Polarization 



Polarization of the muon beams presents a significant physics advantage over the unpolarizcd case, since signal and 
background of electrowcak processes usually come predominantly from different polarization states. 
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1. Polarized muon production 



In the center of mass of a decaying pion, the outgoing muon is fully polarized (P = —1 for /i + and +1 for fi ). 
In the lab system the polarization depends on the decay angle Od and initial pion energy [ 132 - 134 1 . For pion kinetic 



energy larger than the pion mass, the average polarization is about 20%, and if nothing else is done, the polarization 
of the captured muons after the phase rotation system is approximately this value. 

If higher polarization is required, some selection of muons from forward pion decays (cos^ — ► 1) is required. 
Figure [18], above, showed the polarization of the phase rotated muons. The polarization {P> |, — ^ < P < |, and 
P< — i} is marked by the symbols +, o and — respectively. If a selection is made on the minimum energy of the 
muons. then greater polarization is obtained. The tighter the cut, the higher the polarization, but the less the fraction 
-Fsurv of muons that survive. Figure El] gives the results of a Monte Carlo study. 
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FIG. 21. Polarization vs. fraction F SU rv of /x's accepted. 

If this selection is made on both beams, and if the proton bunch intensity is maintained, then each muon bunch 
is reduced by the factor F surv and the luminosity would fall by F^ UIV . But if, instead, proton bunches are merged 
so as to obtain half as many bunches with twice the intensity, then the muon bunch intensity is maintained and the 
luminosity (and repetition rate) falls only as F surv . 

The luminosity could be maintained at the full unpolarized value if the proton source intensity could be increased. 
Such an increase in proton source intensity in the unpolarized case might be impractical because of the resultant 
excessive high energy muon beam power, but this restriction does not apply if the increase is used to offset losses in 
generating polarization. 

Thus, the goal of high muon beam polarization may shift the parameters of the muon collider towards lower 
repetition rate and higher peak currents in the proton driver. 



2. Polarization preservation 



The preservation of muon polarization has been disqusscd in some detail in [ 135 1 
process the muons lose energy in material and have a spin-flip probability V, 
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where 5E/E is the fractional loss of energy due to ionization. In our case, the integrated energy loss is approximately 
3 GeV and the typical energy is 150 MeV, so the integrated spin-flip probability is close to 10%. The change in 
polarization 8V jV is twice the spin-flip probability, so the reduction in polarization is approximately 20%. This loss 
is included in Fig. Ell 

During circulation in any ring, the muon spin, if initially longitudinal, will precess by j(g — 2)/2 turns per revolution, 
where (g — 2)/2 is 1.166 x 10~ 3 . A given energy spread A7/7 will introduce variations in these precessions and cause 
dilution of the polarization. But if the particles remain in the ring for an exact integer number of synchrotron 
oscillations, then their individual average 7's will be the same and no dilution will occur. 

In the collider, bending can be performed with the spin orientation in the vertical direction, and the spin rotated 
into the longitudinal direction only for the interaction region. The design of such spin rotators appears relatively 
straightforward, but long. This might be a preferred solution at high energies but is not practical in the 100 GeV 
machine. An alternative is to use such a small energy spread, as in the Higgs factory, that although the polarization 
vector precesses, the beam polarization does not become significantly diluted. In addition, calibration of the Higgs 
factory collider energy to 1 part in a million |12| requires the spins to precess continuously from turn to turn. 



H. R&D program 



An R&D program is underway to continue theoretical studies (optimization of pion production and capture) and 
to clarify several critical issues related to targetry and phase rotation [136]. A jet of the room temperature eutectic 
liquid alloy of Ga-Sn will be exposed to nanosecond pulses of 1.5 x 10 1 24 GeV protons at the Brookhaven AGS to 
study the effect of the resulting pressure wave on the liquid. The same jet will also be used in conjunction with a 20 T, 
20 cm bore resistive magnet at the National High Magnetic Field Laboratory (Tallahassee, FL) to study the effect of 
eddy currents on jet propagation. Then, a pulsed, 20 T magnet will be added to the BNL test station to explore the 
full configuration of jet, magnet and pulsed proton beam. Also, a 70 MHz rf cavity will be exposed to the intense flux 
of secondary particles downstream of the target and 20 T magnet to determine viable operating parameters for the 
first phase rotation cavity. The complete configuration of the targetry experiment is sketched in Fig. E3. 
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FIG. 22. Plan view of the full configuration of the targetry experiment. 



The first two studies should be accomplished during 1999, and the third and fourth in the years 2000/01. 



V. IONIZATION COOLING 



A. Introduction 



The design of an efficient and practical cooling system is one of the major challenges for the muon collider project. 
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For a high luminosity collider, the 6-D phase space volume occupied by the muon beam must be reduced by a factor 
of 10 5 - 10 6 . Furthermore, this phase space reduction must be done within a time that is not long compared to the 
muon lifetime (/i lifetime « 2 /is). Cooling by synchrotron rad iatio n, conventional stochastic cooling and co nven tional 



electron cooling are all too slow. Optical stochastic cooling [137], electron cooling in a plasma discharge [ 138 1 , and 



cooling in a crystal lattice 13£,14C] are being studied, but ap p ea r tec hnologically difficult. The new method proposed 
for cooling muons is ionization cooling. This technique fl6| , |l8|j20 ,141| is uniquely applicable to muons because of their 
minimal interaction with matter. It is a method that seems relatively straightforward in principle, but has proven 
quite challenging to implement in practice 

Ionization cooling involves passing the beam through some material in which the muons lose both transverse and 
longitudinal momentum by ionization energy loss, commonly referred to as dE/dx. The longitudinal muon momentum 
is then restored by reacceleration, leaving a net loss of transverse momentum (transverse cooling). The process is 
repeated many times to achieve a large cooling factor. 

The energy spread can be reduced by introducing a transverse variation in the absorber density or thickness (e.g. a 
wedge) at a location where there is dispersion (the transverse position is energy dependent). This method results in a 
corresponding increase of transverse phase space and is thus an exchange of longitudinal and transverse emittances. 
With transverse cooling, this allows cooling in all dimensions. 

We define the root mean square rms normalized omittance as 

e i>N = y/{5^){SpD-{SriSpi) 2 /m^c (22) 

where r,; and pi are the beam canonical conjugate variables with i = 1,2,3 denoting the x, y and z directions, and 
(...) indicates statistical averaging over the particles. The operator S denotes the deviation from the average, so that 
Sri = i~i — (r.i) and likewise for Spi. The appropriate figure of merit for cooling is the final value of the 6-D relativistically 
invariant emittance eg.jv, which is proportional to the area in the 6-D phase space (x,y, z,p x ,p yi p z ) since, to a fairly 
good approximation, it is preserved during acceleration and storage in the collider ring. This quantity is the square 
root of the determinant of a general quadratic moment matrix containing all possible correlations. However, until the 
nature and practical implications of these correlations are understood, it is more conservative to ignore the correlations 
and use the following simplified expression for 6-D normalized emittance, 

£6,N ~ £x,N X £y,N X £z,N (23) 

Theoretical studies have shown that, assuming realistic parameters for the cooling hardware, ionization cooling can 
be expected to reduce the phase space volume occupied by the initial muon beam by a factor of 10 5 - 10 6 . A complete 
cooling channel would consist of 20 - 30 cooling stages, each stage yielding about a factor of two in 6-D phase space 
reduction. 

It is recognized that the feasibility of constructing a muon ionization cooling channel is on the critical path to 
understanding the viability of the whole muon collider concept. The muon cooling channel is the most novel part of 
a muon collider complex. Steady progress has been made both in improving the design of sections of the channel and 
in adding detail to the computer simulations. A vigorous experimental program is needed to verify and benchmark 
the computer simulations. 

The following parts of this section briefly describe the physics underling the process of ionization cooling. We 
will show results of simulations for some chosen examples, and outline a six year R&D program to demonstrate the 
feasibility of using ionization cooling techniques. 



B. Cooling theory 

In ionization cooling, the beam loses both transverse and longitudinal momentum as it passes through a material. At 
the same time its emittance is increased due to stochastic multiple scattering and Landau straggling. The longitudinal 
momentum can be restored by reacceleration, leaving a net loss of transverse momentum. 

The approximate equation for transverse cooling in a step ds along the particle's orbit is ]i1||T8| , pO 24,145,142] 



de N 1 dEp ejy /3 ± (0.014GeV) 2 



ds I3 2 ds E„ 2l3 3 E„m^ L R 



(24) 



where (3 is the normalized velocity, En is the total energy, is the muon mass, e^v is the normalized transverse 
emittance, /3j_ is the betatron function at the absorber, dE^/ds is the energy loss per unit length, and Lr is the 
radiation length of the material. The betatron function is determined by the strengths of the elements in the focusing 
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lattice 144 1 . Together with the beam emittance this function determines the local size and divergence of the beam. 
Note that the energy loss dE^/ds is defined here as a positive quantity, unlike the convention often used in particle 
physics. The first term in this equation is the cooling term, and the second describes the heating due to multiple 
scattering. The heating term is minimized if (3±_ is small (strong-focusing) and Lr is large (a low-Z absorber). 

The minimum, normalized transverse emittance that can be achieved for a given absorber in a given focusing field 
is reached when the cooling rate equals the heating rate in Eq. S3 



(3 ± (UMeV) 2 



(25) 



For a relativistic muon in liquid hydrogen with a betatron focusing value of 8 cm, which corresponds roughly to 
confinement in a 15 T solcnoidal field, the minimum achievable emittance is about 340 mm-mrad. 
The equation for energy spread is p8 24 23] 



d(AE^) 
ds 



dE, 
ds 



dE u 



rf(A^) s 2 trag 
ds 



(26) 



where the first term describes the cooling (or heating) due to energy loss, and the second term describes the heating 
due to straggling. AE^ is the rms spread in the energy of the beam. 

Ionization cooling of muons seems relatively straightforward in theory, but will require extensive simulation studies 
and hardware development for its optimization. There arc practical problems in designing lattices that can transport 
and focus the large emittance beam. There will also be effects from space charge and wake fields. 

We have developed a number of tools for studying the ionization cooling process. First, the basic theory was 
used to identify the most promising beam properties, material type and focusing arrangements for cooling. Given 
the practical limits on magnetic field strengths, this gives an estimate of the minimum achievable emittance for a 
given configuration. Next, the differential equations for cooling and heating described above were incorporated into 
a computer code. Allowance for the shifts in the betatron phase advance due to space charge and aberrations was 
included. This code was used to develop an overall cooling scenario, which broke the cooling system into a number of 
stages, and determined the properties of the beam, radio frequency (rf) cavities, and focusing lattice at each stage. 

Finally, several tr ackin g codes were cith er written or modified to study the cooling process in detail. Two new 
codes (SIMUCOOL [145], and ICOOL [146]) use Monte Carlo techniques to track particles one at a time through the 
cooling system. All the codes attempt to include all relevant physical processes to some degree, (e.g. energy loss, 
straggling, multiple scattering) and use Maxwellian models of the focusing fields. They do not yet take into acco unt 
any space charge or wake field effects. In addition, we have also used a modified version of PARM ELA [147] for 



tracking, which does include s pace charge effects, and a double precision version of GEANT [14£,14S]. 

We have recently developed [ 150 1 a model of beam cooling based on a second order moment expansion. A computer 
code solving the equations for transverse cooling gives results that agree with tracking codes. The code is being 
extended to include energy spread and bends. It is very fast and is appropriate for preliminary design and optimization 
of the cooling channel. All of these codes are actively being updated and optimized for studying the cooling problem. 



C. Cooling system 

The cooling is obtained in a series of cooling stages. Each stage consists of a succession of the following components: 

1. Transverse cooling sections using materials in a strong focusing (low (3±) environment alternated with linear 
accelerators. 

2. Emittance exchange in lattices that generate dispersion, with absorbing wedges to reduce momentum spread. 

3. Matching sections to optimize the transmission and cooling parameters of the following section. 

In the examples that follow it is seen that each such stage lowers the 6-D emittance by a factor of about 2. Since 
the required total 6-D cooling is O(10 6 ), about 20 such stages are required. The total length of the system would 
be of the order of 600 m, and the total acceleration required would be approximately 6 GeV. The fraction of muons 
remaining at the end of the cooling system is estimated to be « 60%. 

The baseline solution for transverse cooling involves the use of liquid hydrogen absorbers in strong solenoid focusing 
fields, interleaved with short linac sections. The solenoidal fields in successive absorbers are reversed to avoid build up 
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of the canonical angular momentum [151]. The focusing magnetic fields are small (~ 1 T) in the early stages where 
the emittances are large, but must increase as the emittance falls. 

Three transverse cooling examples have been designed and simulated. The first uses 1.25 T solenoids to cool 
the very large emittance beam coming from the phase rotation channel. The muon beam at the end of the decay 
channel is very intense, with approximately 7.5 x 10 12 muons/bunch, but with a large normalized transverse emittance 
(ea:,A r (i'nis) w 15 X 10 3 mm-mrad) and a large normalized longitudinal omittance (e Zj jv(rms) f=a 612 mm). The second 
would lie toward the end of a full cooling sequence and uses 15 T solenoids. The third, using 31 T solenoids, meets 
the requirements for the Higgs factory and could be the final cooling stage for this machine. 

The baseline solution for emittance exchange involves the use of bent solenoids to generate dispersion and wedges 
of hydrogen or LiH to reduce the energy spread. A simulated example is given for exchange that would be needed 
after the 15 T transverse cooling case. 

A lithium lens solution may prove more economical for the final stages, and might allow even lower emittances 
to be obtained. In this case, the lithium lens serves simultaneously to maintain the low /3j_, and provide dE/dx for 
cooling. Similar lenses, with surface fields of 10 T, were develop ed at No vosibirsk (BINP) and have been used, at 
low repetition rates, as focusing elements at FNAL and CERN 152 - 155f] . Lenses for the cooling application, which 



would operate at 15 Hz, would need to employ flowing liquid lithium to provide adequate thermal cooling. Higher 
surface fields would also be desirable. 

Studies have simulated cooling in multiple lithium lenses, and have shown cooling through several orders of mag- 
nitude |ll| . But these studies have, so far, used ideal matching and acceleration. Cooling is also being studied in 
beam recirculators, which could lead to reduction of costs of the cooling section [ 156 , 157 1 , but full simulations with 
all higher order effects have not yet been successfully demonstrated. 



D. 15 T solenoid transverse cooling example 

The lattice consists of 11 identical 2 m long cells. In each cell there is a liquid hydrogen absorber (64 cm long, 
10 cm diameter) in the 15 T solenoid focusing magnet (64 cm long, 12 cm diameter). The direction of the fields in 
the magnets alternates from one cell to the next. Between the 15 T solenoids there are magnetic matching sections 
(1.3 m long, 32 cm inside diameter) where the field is lowered and then reversed. Inside the matching sections are 
short, 805 MHz, high gradient (36 MeV/m) linacs. Figure shows the cross section of one cell of such a system, 
together with the betatron function, and the magnetic field along the axis. For convenience in modeling, the section 
shown in Fig. p3|(a) starts and ends symmetrically in the middle of hydrogen absorber regions at the location of the 
peaks in the axial magnetic field. In practice each cell would start at the beginning of the hydrogen region and extend 
to the end of the rf module. 

A GEANT simulation of muons travers ing a sec tion of the cooling channel is shown in Fig. 



Additional simulations were performed [ 151 , 158 1 using the program ICOOL. The only likely significant effects which 
are not yet included arc space charge and wakefields. Analytic calculations for particle bunches in free space indicate 
that these effects should, for the later stages, be significant but not overwhelming. A full simulation must be done 
before we are assured that no problems exist. Particles are introduced with transverse and longitudinal omittance 
(186 MeV/c, 1400 tt mm-mrad transverse, and 1100 tt mm longitudinal), together with a number of naturally occurring 
correlations. Firstly, the particles are given the angular momentum appropriate for the starting axial magnetic field. 
Secondly particles with large initial radius r and/or divergence o have longer pathlengths in a solenoidal field and 
tend to spread out with time. This can be parameterized by defining an initial transverse amplitude 
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A' = (27) 
P± 

The temporal spreading can be minimized by introducing an initial correlation between p z and A 2 that equalizes the 
forward velocity of the initial particles. This correlation causes the average momentum of the beam to grow from 
the reference value of 186 MeV/c to i=a 195 MeV/c. Lastly, a distortion of the longitudinal bunch distribution can be 
introduced to reflect the asymmetric nature of the "alpha" -shaped rf bucket. 

Figure p5|(a) shows the average momentum of the beam as a function of distance along the channel. The momentum 
drops as the beam crosses the liquid hydrogen absorbers. The gradient and phase of the rf cavities have been adjusted 
so that the reaccelcration given to the reference particle equals the mean energy loss. This causes the average 
momentum of the beam to remain in a narrow band around 195 MeV/c. Figure p5|(b) shows the mechanical and 
canonical angular momenta as a function of distance along the channel. The mechanical angular momentum shows 
the rotational motion of the beam around the axial solenoidal field. It periodically reverses sign when the solenoids 



alternate direction. The canonical angular momentum is defined such that it removes the axial field dependence [158 
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FIG. 23. (a) Cross section of one half period of an alternating solenoid cooling lattice; (b) axial magnetic field vs. z; (c) 
(3± function vs. z. 
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FIG. 24. GEANT simulation of muons traversing a section of the alternating solenoid cooling channel. The variation of the 
magnetic field B z is shown for l| cells of the figure. 
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Without the absorbers, the beam would have a constant (0) value for the canonical angular momentum. However, 
the presence of absorbers causes the canonical angular momentum to grow and would lead to severe emittance growth 
by the end of a long channel. This growth is stopped by alternating the direction of the solenoid field, as shown 
in Fig. [1(b) 

Simulations have shown that 2 m is a reasonable (half) period for the field, since the net growth in 
canonical angular momentum is small. In addition synchrobetatron resonances are avoided since the periodicity of 
the field forces the average betatron wavelength to be 2 m, whereas the synchrotron oscillation wavelength seen in 
the simulations for this arrangement is ~ 14 m. 
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FIG. 25. a) Average momentum vs. z; b) Average angular momentum: mechanical (solid curve) and canonical (dashed 
curve), vs. z. 

Figure |2^(a) shows the rms and maximum radius of any particle in the beam distribution as a function of distance 
along the channel. The rms radius shows that most of the beam is confined to within 2 cm of the axis. The peak 
rms radius decreases towards the end of the channel as a result of the cooling. The maximum particle radius is about 
8 cm, which determines the radius of the windows required in the rf cavities. Figure |26|(b) shows the rms momentum 
spread corrected for the correlation between p z and transverse amplitude imposed on the initial particle distribution. 
The momentum spread grows as a function of distance since the alternating solenoid system only cools the transverse 
emittance. Figure ^6|(c) shows the rms bunch length s a function of distance along the channel. Again this grows 
with distance since this channel does not cool longitudinally. 

Figure |2?](a) shows the decrease in transverse normalized emittance as a function of distance along the channel. 
The system provides cooling by a factor of w 2 in both the x and y transverse phase spaces. From the changing slope 
of the curve we note that the rate of cooling is dropping. This sets w 22 m as the maximum useful length for this 
type of system. It must be followed by a longitudinal emittance exchange region to reduce the momentum spread and 
bunch length approximately back to their starting values. Figure ^(b) shows the increase in longitudinal normalized 
emittance in the channel due to the increase in momentum spread and bunch length. Finally, Fig. |27](c) shows the 
decrease in the 6-D normalized emittance as a function of distance along the channel. There is a net decrease in 6-D 
emittance by a factor of f=a 2 in the channel. Table ^ gives the initial and final beam parameters. 



Th is simulation has been confirmed, with minor differences, by double precision GEANT [ 149 1 and PARMELA 
[1471 codes. 
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FIG. 27. Emittance vs. z: a) transverse emittance; b) longitudinal emittance; and c) 6-D emittance. 



E. 31 T solenoid transverse cooling example 



As in the preceding example, the lattice consists of 11 identical 2 m long cells with the direction of the fields in the 
solenoids alternating from one cell to the next. The maximum solcnoidal field is higher (31 T) than in the previous 
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TABLE V. Initial and final beam parameters in a 15 T transverse cooling section. 







initial 


final 


final/initial 


Particles tracked 




1000 


980 


0.98 


Reference momentum 


MeV/c 


186 


186 


1.0 


Transverse Emittance 


7r mm-mrad 


1400 


600 


0.43 


Longitudinal Emittance 


7r mm-mrad 


1100 


2300 


2.09 


6-D emittance xlCT 12 


(-7T m-rad) 


2000 


800 


0.40 


rms beam radius in hydrogen 


cm 


U.o 


U.00 


n fin 


rms beam radius in linac 


cm 


2.0 


1.4 


0.70 


max beam radius in linac 


cm 


7.0 


7.0 


1.0 


rms bunch length 


cm 


1.5 


2.2 


1.5 


max bunch full width 


cm 


13 


19 


1.5 


rms Ap/p 


% 


3.8 


5.6 


1.5 



example, but the bore is smaller (8 cm), and the liquid hydrogen absorber has a smaller diameter (6 cm). Between the 
31 T solenoids there are 1.3 m long matching sections with an inside diameter of 32 cm, superimposed on a 36 MeV/m 
reacceleration linac operating at 805 MHz. 

Table VI gives the initial and final parameters for the 31 T example, together with the required emittances for a 
Higgs factory. In setting these requirements a dilution of 20% during acceleration is assumed in each of the three 
emittances. 



TABLE VI. 


Initial and final beam parameters in 


a 31 T transverse coolin 


g section. 








initial 


final 


final/initial 


required 


Particles tracked 




4000 


3984 


0.99 




Reference momentum 


MeV/c 


186 


186 






Transverse Emittance 


7r mm-mrad 


460 


240 


0.52 


240 


Longitudinal Emittance 


7r mm-mrad 


850 


1600 


1.9 




6-D emittance xl0~ 12 


(tt m-rad) 3 


150 


95 


0.63 


98 


rms beam radius in hydrogen 


cm 


0.44 


0.33 


0.75 




rms beam radius in linac 


cm 


.4 


1.1 


0.80 




max beam radius in linac 


cm 


6.0 


6.0 


1.0 




rms bunch length 


cm 


1.5 


1.8 


1.2 




max bunch full width 


cm 


11 


19 


1.7 




rms Ap/p 


% 


3.5 


5.0 


1.4 





F. Bent solenoid emittance exchange example 



We have been considering using a system that exchanges longitudinal and transverse emittance by exploiting 
dispersion in a large acceptance channel, with a low-Z wedge absorber in the region of dispersion. 

In a bent solenoid, in the absence of any dipolc field, there is a drift perpen dicu lar to the bend plane of the center of 
the Larmor circular orbit, which is proportional to the particle's momentum [159]. In our example we have introduced 
a uniform dipolc field over the bend to cancel this drift exactly for particles with the reference momentum. Particles 
with momenta differing from the reference momentum then spread out spatially, giving the required dispersio n (0. 4 
m). The dispersion is removed, and the momentum spread reduced, by introducing liquid hydrogen wedges [ 160 1 . 
The hydrogen wedges would be contained by thin beryllium or aluminum foils, but these were not included in this 
simulation. 

After one bend and one set of wedges, the beam is asymmetric in cross section. Symmetry is restored by a following 
bend and wedge system rotated by 90 degrees with respect to the first. Figure ^8] shows a representation of the two 
bends and wedges. The total solenoid length was 8.5 m. The beam tube outside diameter is 20 cm, and the minimum 
bend radii is 34 cm. 

Figure ^9|(a) shows the magnetic fields (B z , B y , and B x ) as a function of the position along the cell. The solenoid 
bend curvature is exactly that given by the trajectory of a reference particle (equal in momentum to the average 
momenta given in Fig. ^(b)) in the given transverse fields. The actual shape of the bend turns out to be very 
important. Discontinuities in the bend radius can excite perturbations which increase the transverse emittance. We 
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FIG. 28. Representation of a bent solenoid longitudinal emittance exchange section 



have shown, for example, that the transverse emittance growth in a bent solenoid depends on discontinuities of the 
bend radius as a function of distance, and its first and second derivatives, the size and tilt of the solenoidal coils, 
auxiliary fields and the 6-D phase space of the beam. Thus optimization is not straightforward. One solution to this 
problem is to have long, adiabatic bends. However, this adds undesirable length to the emittance exchange section. 
We are studying options with coupling sections to tight bends roughly half a Larmor length long, which seems to 
minimize transverse emittance growth while also minimizing the length of the section. Due to similar problems, 
the length and longitudinal distribution of the wedge material has also been found to affect emittance growth. For 
example, the growth can be minimized when the vector sum of the Larmor phases at the absorber elements is small 
or zero. 

The simulations were performed using the program ICOOL. The maximum beam radius is 10 cm. Transmission was 
100%. Figure |3^(a) shows the rms longitudinal momentum spread relative to the reference momentum as a function 
of the position along the cell. The fractional spread decreases from an initial value of approximately 5%, to a final 
value of approximately 2.2%. At the same time, since this is an emittance exchange, the transverse beam area grows, 
as shown in Fig. ^0|(b). One notes that the area increases not only in the regions of bends (region 1 and 8), but also 
in the regions of wedges (2-6 and 9-11). This is probably due to failures in matching that have yet to be understood. 

Figure ^l] shows scatter plots of the transverse particle positions against their momenta. The dispersion is clearly 
observed in Fig. pl|(b) (after the first bend) and in Fig. 0(e) (after the second). It is seen to be removed, with a 



corresponding decrease in momentum spread, in Fig. |^(cj (after the first set of wedges) and Fig. [5l](f) (after the 
second set of wedges). 

Figure [32] shows a scattcrplot of the square of the particle radii vs. their longitudinal momenta, (a) at the start, 
and (b) at the end of the omittance exchange section. The decrease in momentum spread and rise in beam area are 
clearly evident. 



The initial and final beam parameters are given in table VII. Although this example demonstrate a factor of 



reduction in the longitudinal momentum spread, there is a 37% increase in the 5-D phase space. The simulations 
must be extended to include rf so that the 6-D emittance can be studied and the emittance exchange section can be 
optimized. 

Emittance exchange in solid LiH wedges, with ideal dispersion and matching, has also been successfully simulated 
using SIMUCOOL [161]. Dispersion generation by weak focusing spectrometers [156 and dipolcs with solenoids [162] 
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FIG. 29. a) Axial B z and dipole B y , B x , magnetic fields; b) average momentum; both as a function of the position along 
the cell. 
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have also been studied. 
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a) region 1 
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b) region 2 
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c) region 7 
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d) region 7 
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e) region 9 
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f) region 12 




120 140 160 180 200 220 

Pz (MeV/c) 



FIG. 31. y vs. p z plots: a) at the start, b) after the first bend, c) after the first set of wedges, x vs. p z plots: d) after the 
first wedges, e) after the second bend, and f) at the end of the emittance exchange section, following the second set of wedges. 
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FIG. 32. Scatterplot of squared radii vs. longitudinal momentum: a) at the start, and b) at the end of the emittance exchange 
section. 
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TABLE VII. Initial and final beam parameters in a longitudinal omittance exchange section. 



initial final final/initial 

Longitudinal Momentum spread MeV/c 9.26 3.35 0.36 

Ave. Momentum MeV/c 180 150 0.83 

Transverse size cm 1.33 2.26 1.70 

Transverse Momentum spread MeV/c 6.84 7.84 1.15 

Transverse Emittance tt mm-mrad 870 1694 1.95 

Emit^ rans x Apj (tt m-mrad) 2 MeV/c 7.0 9.6 1.37 



G. RF for the cooling systems 



The losses in the longitudinal momentum of the muon beam from the cooling media have to be restored using rf 
acceleration sections. These rf structures are embedded in solenoidal fields that reverse direction within each section. 
In the two transverse cooling examples above, the rf frequency is 805 MHz and the peak gradient is 36 MeV/m. The 
magnetic fields that extend over the cavities vary from to 10 T, reversing in the center. It should be pointed out 
that in the earlier stages, the bunches are longer, and lower frequencies will be required. 




FIG. 33. Two full cell sections plus two half cell sections of the interleaved tt/2 mode accelerating cavities. The volumes 
labeled C are powered separately from the volumes labeled D. 

In order to realize maximum accelerating gradients within the acceleration cavities, we take advantage of the 
penetrating properties of a muon beam by placing thin windows between each rf cell, thereby creating an accelerating 
structure closely approximating the classic pill-box cavity. This permits operating conditions in which the axial 
accelerating field is equal to the maximum wall field and gives a high shunt impedance. 

The windows in the 15 T example are 16 cm diameter, 125 /mi thick Be foils. In the 31 T case, they are 10 cm 
diameter and 50 /jm thick. 

For these rf structures, we will use an interleaved cavity design in which two parts are independently powered 
(Fig. [33]). The mode of the system will be referred to as ir/2 interleaved. Each section supports a standing wave tt 
mode, with each acceleration cell tt/2 long, giving a good transit time factor. To reduce the peak rf power requirements 
(by a factor of 2), we are considering operating the cells at liquid nitrogen temperatures. 



The characteristics of the rf systems currently being studied are summarized in Table VIII . Figure t54j shows a full 



1.3 m section with interleaved cavities. Each cell is 8.1 cm in length and the 1.3 m section consists of 16 cells. 
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FIG. 34. A 1.3 m acceleration section with quarter section cut away for viewing the inter-cavity windows. 



TABLE VIII. Characteristics of the rf system. 



RF frequency [MHz] 
Cavity Length [cm] 
Cavity Inner Radius [cm] 
Cavity Outer Radius [cm] 
Q/1000 

Peak Axial Gradient [MV/m] 

Shunt Impedance [Mfi/m] 

Zt 2 [Mfi/m] 

Fill Time (3 r) \ps] 

RF Peak Power [MW/m] 

Ave. Power (15Hz) [KW/m] 

Be window aperture [cm] 

Be window thickness [/im] 



805 

8.1 

14.6 

21 
2 x 20 

36 
2 x 44 
2 x 36 
2 x 12 
i x29 

5.3 

16 (10 for 31 T case) 
127 (50 for 31 T case) 



H. The liquid lithium lens 



The final cooling clement ultimately determines the luminosity of the collider. In order to obtain smaller transverse 
emittance as the muon beam travels down the cooling channel, the focusing strength must increase, i.e. the /3j_'s 
must decrease. A current within a conductor produces an active lens absorber, which can maintain the beam at small 
/3j_ throughout an extended absorber length, while simultaneously attenuating the beam momentum. An active lens 
absorber, such as a lithium lens, may prove to be the most efficient cooling clement for the final stages. 

The cooling power of a Li lens is illustrated in Figure [35[ where the x vs. p x phase space distributions at the 
beginning and at the end of the absorber are shown. This example corresponds to a 1 m long lens, with 1 cm 
radius, and a surface field of 10 T. The beam momentum entering the lens was 267 MeV/c, with Gaussian transverse 
spatial and momentum distributions: a x = a y = 2.89 mm, a Vi = a Py = 26.7 McV/c, and a normalized emittance 
of e x ,N = 710 mm-mrad. The normalized emittance at the end of the absorber was e x ^ = 450 mm-mrad (cooling 
factor ~ 1-57), and the final beam momentum was 159 MeV/c. The results were obtained using a detailed GEANT 
simulation of a single stage. 

An alternative cooling scheme under study uses a series of Li lenses. The lens parameters would have to vary to 
match the changing beam emittance along the section and in addition, acceleration of the beam between the lenses 
has to be included. 

Lithium lenses have been used with high reliability as focusing elements at FNAL and CERN p, 153| ,154|. Although 
these lenses have many similar properties to those required for ionization cooling, there are some very crucial differences 
which will require significant advances in lens technology: ionization cooling requires longer lenses (~ 1 meter), higher 
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FIG. 35. x-p x phase space distribution at the beginning and at the end of the absorber described in the text. 

fields (~ 10 T), and higher operation rates (15 Hz). The last requirement calls for operating the lenses with lithium 
in the liquid phase. 

A liquid Li lens consists of a small diameter rod-like chamber filled with liquid Li through which a large current is 
drawn. 

The azimuthal magnetic field focuses the beam to give the minimum achievable emittance e x ,N ~ Cf3± where the 
constant C depends on the properties of the material, for example, Cl% — 79 mm-mrad/cm. The focusing term can 
be written as ~ 0.08[cm] y/p/ J withp is the muon momentum in MeV/c, and J is the current density in MA/cm 2 . 
Increasing J is obviously desirable. Decreasing p can also be useful. However, below about 250 MeV/c the slope 
of ^L(E) tends to increase the longitudinal emittance. The requirement for the highest current density causes large 
ohmic power deposition. The current density will be limited by the maximum tolerable deposited energy, which will 
produce instantaneous heating, expansion, and pressure effects. Understanding these effects is part of the ongoing 
liquid Li lens R&D. 

The structural design of the lithium lens is determined by how the pressure pulse and heat deposition are handled. 
We assume that the Li will be flowing rapidly under high pressure, confined by electrical insulators radially and by 
fairly thick Be windows longitudinally. Operation at 15 Hz for long periods poses severe challenges. Shock, fatigue and 
other failure modes are being evaluated, in addition to studies of material compatibilities, corrosion and degradation 
to insure safe operation over long periods. It seems that the minimum required radius of the lens may be the most 
important parameter to determine, since mechanical problems rise while losses decrease as a function of radius. 

Transferring the beam from one lens to another, with linacs to reaccelerate and provide longitudinal focusing, is 
also a challenging problem, because of the multiple scattering introduced in the windows, straggling and the large 
divergence of the beams. We are in the process of evaluating a number of designs for this transfer channel, using 
detailed tracking simulations that include solenoids, quadrupolcs and other focusing elements together with Li lenses. 

A group from BINP has designed, and is constructing, a 15 cm long liquid lithium lens prototype that will eventually 
be tested at FNAL. It is planned to extend this R&D program to design, construct, and test longer lenses. The design 
of two lenses, whose behavior will be tested at first on a bench and then with muon beams at the Ionization Cooling 



Demonstration Facility, will then follow [ 163 1 
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I. Ionization cooling experimental R&D 



An R&D program has been proposed to design and prototype the critical sections of a muon ionization cooling 
channel. The goal of this experimental R&D program is to develop the muon ionization cooling hardware to the point 
where a complete ionization cooli ng ch annel can be confidently designed for the First Muon C ollid er. Details can be 



found in Fermilab proposal P904 [162]. A summary of the R&D program can be found in ref. [164 



The proposed R&D program consists of: 

• Developing an appropriate rf re-acceleration structure. It is proposed to construct a 3-cell prototype rf cavity 
with thin beryllium windows, which will be tested at high power and within a high-field solenoid. 

• Prototyping initially a 2 m section, and eventually a 10 m section, of an alternating solenoid transverse cooling 
stage. It is proposed to test the performance of these sections in a muon beam of the appropriate momentum. 

• Prototyping an emittance exchange (wedge) section and measuring its performance in a muon beam of the 
appropriate momentum. 

• Prototyping and bench testing ~lm long liquid lithium lenses, and developing lenses with the highest achievable 
surface fields, and hence the maximum radial focusing. 

• Prototyping a lithium lens-rf-lens system and measuring its performance in a muon beam of the appropriate 
momentum. 

• Developing, prototyping, and testing a hybrid lithium lens/wedge cooling system. 

The measurements that are needed to demonstrate the cooling capability and optimize the design of the alternating 
solenoid, wedge, and lithium lens cooling stages will require the construction and operation of an ionization cooling 
test facility. This facility will need 

1. a muon beam with a central momentum that can be chosen in the range 100-300 MeV/c, 

2. an experimental area that can accommodate a cooling and instrumentation setup of initially ~ 30 m in length, 
and eventually up to ~ 50 m in length, and 

3. instrumentation to precisely measure the positions of the incoming and outgoing particles in 6-D phase space 
and confirm that they are muons. 

In the initial desig n sho wn in Fig. [}6|, the instrumentation consists of identical measuring systems before and after the 
cooling apparatus [165]. Each measuring system consists of (a) an upstream time measuring device to determine the 
arrival time of the particles to one quarter of an rf cycle (~ ±300 ps), (b) an upstream momentum spectrometer in 
which the track trajectories arc measured by low pressure TPC's on either side of a bent solenoid, (c) an accelerating 
rf cavity to change the particles momentum by an amount that depends on its arrival time, (d) a downstream 
momentum spectrometer, which is identical to the upstream spectrometer, and together with the rf cavity and the 
upstream spectrometer forms a precise time measurement system with a precision of a few ps. The measuring systems 
are 8 m long, and are contained within a high-field solcnoidal channel to keep the beam particles within the acceptance 
of the cooling apparatus. 

It is proposed to accomplish this ionization cooling R&D program in a period of about 6 years. At the end of this 
period we believe that it will be possible to assess the feasibility and cost of constructing an ionization cooling channel 
for the First Muon Collider, and if it proves feasible, begin a detailed design of the complete cooling channel. 



VI. ACCELERATION 



A. Introduction 



Following cooling and initial bunch compression, the beams must be rapidly accelerated. In this section some of 
the options in accelerator design will be described and examples of acceleration scenarios presented. 

Separate acceleration scenarios are given here for a low momentum spread 100 GcV First Muon Collider (Higgs 
factory), and for a high luminosity 3 TeV collider. Ideally, though more difficult, this accelerator designed for the low 
energy machine should be extendable to the ~ 250 GeV beam energy and from there to the ~ 2 TeV beam energy 
needed for a very high-energy collider. 
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FIG. 36. Schematic of the cooling test apparatus arrangement. 



While acceleration of muons to high energy is clearly possible, an optimal and cost-effective acceleration complex 
is needed. In the scenarios described below, a low-frequency linac would take the beam from the end of cooling to 
an energy of wl GeV followed by recirculating- linac systems to take the beam to 50-70 GeV. The multi-TeV energy 
regime can be reached through a series of very rapid cycling synchrotrons. Variations on the acceleration model and 
potential difficulties are discussed, including the use of Fixed Field Alternating Gradient (FFAG) accelerators in place 
of, or together with, the recirculating-linacs. Finally, topics for further study and research are described. 



B. Accelerator options 



The acceleration time is limited by muon decay (r M = 2.2/^ s at rest) and requires that: 

eV' rf » 0.16MeV/m 

where eV^ is the acceleration rate. An acceleration rate whose value of 0.16 MeV/m is low for a linac, but very high 
for a conventional synchrotron. 

At the lowest energies (< 700 MeV), the momentum spread and beam sizes are so large that only a linac is feasible, 
and acceleration to full energy in a single-pass linac would be good, but it would be very expensive. 

Thus, following the initial linac, some form of recirculating acceleration is preferred. A synchrotron would be 
possible, in principle, but the acceleration must occur so rapidly that conventional magnet ramping is unlikely to be 
practical. Two alternative multi pass methods are being considered: recirculating linac accelerators similar to those 
used at TJNAF and Fixed Field Alternating Gradient (FFAG) accelerators. 

In a recirculating linac accelerator, the beam is circulated through the same linac for several passes, with separate, 
energy matched, fixed-field return paths for each pass. Each return path is optically independent and can be separately 
designed. In the initial lattice design for the muon recirculating linac accelerator, the return arcs are similar alternating 
gradient (AG) systems with the same dipole layout, but with differing quadrupole strengths to allow separate tuning 
and chronicity matching in each arc. Multiple aperture superconducting magnets have also been designed which would 
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reduce the- diameter of the recirculating linac, lowering muon loss from decay and possibly being more 1 economical (see 
Fig 7.12 44 1 ) . In cither case, both the linac and return transports must accommodate large transverse emittances 
(rms) of w 300 it mm-mrad. Strong focusing is required not only to keep apertures down, but also to minimize orbit 
deviations due to the large momentum spreads, which in the initial stages of acceleration, can be as large as 10% rms. 

More recently, an adap tation of the FFAG (Fixed-Field Alternating Gradient) accelerator concept has been proposed 
for (i+ fx— acceleration [ 166 , 167 1 . In this variation, return transports are designed with a very large (factor of 5-10) 
energy acceptance, so that separate energy turns can pass through the same fixed-field elements. More acceleration 
turns are possible than with a recirculating linac accelerator which reduces the rf requirements; but the orbits and 
focusing properties are now energy-dependent. Such FFAG configurations require strong superconducting magnets 
with large apertures to accommodate the energy-dependent spread in closed orbits. The extra cost associated with 
increased magnet apertures must be evaluated against potential savings in the number of magnets and reduced rf/turn 
requirements. 

For the higher energy stages, the muon life time is greater and the needed rate of energy increase is less. Thus, 
above a few hundred GeV, rapid accelerating synchrotrons become possible. In a rapid accelerating synchrotron, the 
beam is also multi-pass accelerated through an rf system, but the beam returns in a single arc, as the magnetic field is 
ramped to match the increase in beam energy. As above, more acceleration turns are possible than with a recirculating 
accelerator, but we now have a single moderate aperture return transport. But rapid cycling synchrotrons have higher 
power costs and the technical challenges associated with the rapid acceleration time needed to conserve muons. 

Thus a complete system would likely include an initial linac followed by a sequence of recirculating linac acceler- 
ators and/or Fixed Field Alternating Gradient machines. Depending on the final collider energy, one or more rapid 
acceleration synchrotrons would follow. Each system increases the beam energy by a factor of 5-10. 



C. Scenario examples 



Several scenarios have been discussed earlier [ 168 , 169 1 , see for instance the parameters (table XI ) used in a simulation 
of longitudinal motion discussed below. The ones given here are more recent, and more detailed, but they should 
not be taken to be definitive. They are examples that were derived to probe the design problems and to show that 
solutions should be possible. 



1. Acceleration for Higgs collider 



Table IX 



gives an example of a sequence of accelerators for a 100 GeV Higgs Factory, i.e. a machine with very low 
momentum spread (0.003%, see table [j]) and relatively large rms transverse emittance 300 tt mm-mrad). 

Following initial linacs, recirculating accelerators are used. The number of arcs in each recirculating accelerator 
is about 10. In this example, conventional fixed field 2 T magnets are used, but the effective ramp frequencies that 
would be needed if pulsed magnets were used are given for reference. 

In this example, all the accelerating cavities are room temperature copper structures, and the accelerating gradients 
are modest (< 10 MeV/m). Nevertheless, the acceleration is rapid enough that the total losses from decay are only 
30%. The heating from these decays is also modest (w 10 W/m) because of the small number of turns and relatively 
low energy. Since no superconducting magnets or rf are used in this example, this heating should cause no problem. 

In this machine, the transverse emittances are large and strong focusing is thus required, but the maximum mo- 
mentum spread is moderate (up to 1.37 rms in the first recirculator) and is thus not likely to be a problem. 

If the same machine is to also run at a high luminosity, with larger momentum spread, then although the six 
dimensional omittance is the same as in the Higgs collider discussed above, the transverse emittance is smaller 
(~ 907T mm-mrad rms ) , and the longitudinal omittance larger (by about a factor of 4) and the momentum spread in 
the first recirculating accelerator would be nearly 6% rms; or about 50% full width. This is a very large momentum 
spread that could only be accepted in FFAG like lattices as discussed below. 

Similarly, if the same acceleration is to be usable as the front end for a 250 + 250 GeV or higher energy machine 
then the transverse emittance will again be less, the longitudinal emittance even larger, and the problem of very large 
momentum spread will be worse. Clearly, although not absolutely needed for a Higgs Factory, it is desirable to solve 
this problem even in that First Collider. 

A separate parameter set for the high luminosity 50 + 50 GeV collider and a 250 + 250 GeV collider could have 
been presented, but their parameters arc very similar to those of the front end of the 3 TeV machine given below, 
and arc thus omitted here. 
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TABLE IX. Accelerator parameters for a Higgs Factory (100 GeV) 



Acc. type 




linac 


linac 


recirc 


recirc 


recirc 


sums 


Magnet type 








warm 


warm 


warm 




rf type 






Ou 




LiU 


Ou 






(GeV) 


0.10 


0.20 


0.70 


2 


7 




^final 


(GeV) 


0.20 


0.70 


2 


7 


50 




Circ. 


(km) 


0.04 


0.07 


0.07 


0.19 


1.74 


2.11 


Turns 




1 


1 


8 


10 


11 




Loss 


(%) 


2.31 


3.98 


7.27 


7.91 


13.94 


31.06 


Decay heat 


(W/m) 


0.89 


1.98 


11.04 


12.99 


12.44 




^fixed 
Ramp freq. 


(T) 

/ITT \ 

(kHz) 






2 

281 


2 

79.83 


2 

8.00 




Disp 


(m) 






1 


1.50 


3 




Anax 


(m) 


0.83 


1.42 


3.00 


5.31 


21.08 




oi nlt 


(cm) 


2.71 


2.22 


1.42 


1.64 


0.90 




Ap/p init 


(%) 


3.58 


2.80 


1.64 


0.56 


0.32 




"is 


(cm) 


1.09 


1.14 


1.01 


0.85 


0.94 




0~x 


(cm) 






1.93 


1.19 


1.34 




Pipe full height 


(cm) 


10.92 


11.42 


10.14 


8.49 


9.39 




Pipe full width 


(cm) 


10.92 


11.42 


19.28 


11.94 


13.35 




rf freq. 


(MHz) 


200 


200 


200 


200 


400 




Acc/turn 


(GeV) 


0.20 


0.40 


0.17 


0.50 


4 




Acc time 


(fis) 






1 


6 


62 




V 


(%) 


5.15 


5.36 


6.36 


2.84 


6.92 




Acc. Grad. 


(MV/m) 


8 


8 


8 


10 


10 




Synch, rot's 




0.62 


0.63 


0.62 


3.92 


23.16 




Cavity rad. 


(cm) 


54.37 


54.88 


54.88 


60.47 


38.26 




Beam time 


(ms) 


0.00 


0.00 


0.00 


0.01 


0.06 




rf time 


H 


0.17 


0.17 


0.17 


0.18 


0.13 




Tot. peak rf 


(GW) 


0.05 


0.10 


0.05 


0.26 


4.71 


5.17 


Ave. rf power 


(MW) 


0.14 


0.24 


0.13 


0.68 


9.54 


10.73 


rf wall 


(MW) 


0.64 


1.16 


0.46 


2.42 


28.06 


32.75 



2. Acceleration for 3 TeV collider 

For a high energy machine, the muon accelerators are physically the largest component and are also probably 
the most expensive. More work is nccclcd on its design. Table [x| gives an early example of a possible sequence of 
accelerators for a 3 TcV collider. 

Linacs are used up to 700 MeV, followed by recirculating linac accelerators. In the first of these, because of the 
very large longitudinal cmittance, the momentum spread as the beam enters the first recirculating linac accelerator 
is 8.5% rms, which is very large. The lattice must have very strong focusing, small dispersion and large aperture. If 
this is not possible, higher energy linacs or lower frequenc y rf could relieve the requirement. 



For the final three stages, pulsed magnet synchrotrons [17C ] are used. In the 200 GeV ring, all the magnets in the 



ring arc pulsed, but in the last two rings a superconducting-pulsed hybrid solution is used. In these cases, if only 
pulsed magnets were used, then the power consumed would be too high, and because only low pulsed fields could 
be used, the circumferences would also be ve ry la rge. It is thus proposed to use rings with alternating warm pulsed 



magnets and superconducting fixed magnets [171] (sec figure |37|). The fixed magnets are superconducting at 8 T; the 
pulsed magnets are warm with fields that swing from - 2 T to + 2 T. The effective ramp frequency is given in the 
table. Both of these rings are in the same tunnel, with the fraction of magnet length pulsed (vs. fixed) being different 
(73% and 43%). 

In all the final three rings superconducting rf is employed to minimize the peak power requirements and to obtain 
high wall to beam efficiency and thus keep the wall power consumptio n rea sonable. In the final two rings the frequency 



and cavity designs are chosen to be the same as that in the TESLA |172| proposal 
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FIG. 37. Schematic of hybrid superconducting-pulse magnet accelerator ring 



D. Design issues 



1. Recirculating linac accelerator lattice issues 



Beam transport R&D for recirculating linac accelerators follows the model of the Thomas Jefferson National Accel- 
erator Facility (TJNAF). The layout is a racetrack with linacs in the straight sections and multi-pass return arcs. At 
the ends of the linacs the multi-pass beam lines are recombined. A pulse magnet at each separation/recombination 
point is used to guide the beam into the energy matched return arc. Some initial lattice design concepts for recir- 
culating linac accelerators are being developed. The basic return arc unit would be a FODO lattice, but with the 
quadrupole strengths varied in order to perturb the arc dispersion function and obtain nearly isochronous motion 
around the arcs. The arcs are dispersion matched by setting the arc phase advance to a multiple of 2n. Arc designs 
based upon the flexible momentum compaction (FMC) module can also be used. 

In the special case of the very low momentum spread Higgs factory, the transverse emittances are very large 
(~ 300 7r mm-mrad rms), and will require strong focusing in the lattices. Momentum acceptance in the rings is, in 
this case, not a problem. But the longitudinal phase space of the muons in the other machines is much larger and 
requires, at low energies, cither long bunches, or large momentum spreads. The requirement of high accelerating 
gradients argues for high frequencies, and thus short bunches. One therefore needs accelerators with large momentum 
acceptances. 

In the 3 TeV example above, the acceptance at injection into the first recirculating accelerator is 8 .5% rms. This is 
very large by conventional standards, but far less than that in the FFAG lattices being studied [ 166 1 . Thus the early 
return arcs of such a recirculating linac accelerator would have to have very strong focusing, and be FFAG-like. Of 
course, if a true FFAG accelerator were to be used for its avoidance of the switchyards and multi- aperture magnets, 
then the specified momentum spread would certainly not be a problem. 

Permanent, ferric or superferric («2 T), or high field superconducting magnets could be used for recirculating linac 
accelerators. The lower field magnets may be economic for initial turns, while high field magnets minimize particle 
travel times, and therefore decay losses. Designs for multi-aperture superconducting magnets suitable for recirculating 
linac accelerators have been developed |44j] , and superconducting magnets with as many as 18 apertures with 0.7-7 T 
fields have been designed. A variety of magnet configurations can be developed; cost/performance optimization will 
be needed in developing a final choice. 
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TABLE X. Parameters of Acceleration for 3 TeV Collider 
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2. RF peak power requirements 

Because of the need for rapid acceleration, the peak rf powers are high, and the resulting numbers of power 
sources large. For the linacs and early recirculating accelerators, the powers are high because of the high gradients 
and low frequencies needed to accelerate the long bunches. At these frequencies («200 MHz), currently available 
sources (triodes and tetrodes) have relatively low maximum output p ower s and are expensive. Low temperature 
operation of the cavities, and superconducting or conventional SLED |l73fl systems, which would reduce the peak 
power requirements, are being considered. 

Study of the example suggests that in the first two recirculators (up to 7 GeV) there is no hope for the rf to keep 
up with the beam loading. The cavity can only be filled in a suitable filling time (twice the time constant in these 
examples), and the rf voltage allowed to sag as the beam makes its multiple passes. If excessive sensitivity to beam 
current is to be avoided, then the stored energy must be large compared to that used, which is somewhat inefficient. 

In the final recirculating accelerator, continuous filling (cw) is just possible, but requires yet higher peak power 
(~ 5 GW total at 400 MHz) because of the high acceleration rate. The use of superconducting cavities can reduce 
losses, and thus reduce this peak rf power source requirement, and was included in the above 3 TeV example. At this 
frequency (400 MHz), klystrons are available with greater power (« 20 MW) than that of the sources at the lower 
frequencies, but a yet higher power klystron (50-100 MW) could probably be developed and would be desirable. 
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3. Pulsed magnet systems 



A pulsed current 4 T magnet has been designed for acceleration to 250 GeV in 360 /is [174 



but efficiency favors 

use of ferric materials in rapid acceleration magnets, although this would limit peak magnetic fields to «2 T. The 
average field can be increased by interleaving magnets swinging from -2 T to 2 T with fixed field 8 T superconducting 
magnets. 

Faster pulsing magnets would require special materials to minimize energy losses from eddy currents. Options 
include silicon steel, metglass laminations, or finemct laminated tape or powdered solid. A 30um mctglass lamination 
suitable for several kHz cycling has been developed. A design of suitable pulsed magnets [ 171 1 (see Fig. [38]) has been 
shown to have sufficiently low losses for this application. The magnets employ cables made of many fine i nsulated 
strands (litz cable) and the yokes are made of very thin (0.28 mm), 3% Si-Fc laminations, possibly of metglas [175,176| 
for the higher rate cases. Detailed designs must be developed and prototypes constructed and the practical limits of 
recycling scenarios should be determined. 
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FIG. 38. A 2-D picture of an H frame magnet lamination with grain oriented 3% Si-Fe steel.The arrows show both the 
magnetic field and the grain direction. 



4- Superconducting linacs 

While the gradients needed in the acceleration systems are not excessive, they are larger than previous experience 
at the lower frequencies. The high peak power pulsed operation poses power handling difficulties at lower energies and 
high peak current presents collective effect (wakcficld) difficulties at higher energies. Higher gradients and efficiencies 
in all sections would improve performance. 

The superconducting rf would operate in pulsed mode, matched to the acceleration time of up to a few ms. This 
pulse structure is similar to the multibunch acceleration mode planned for TESLA (25 MV/m at 1300 MHz designs), 
and studies indicate that this design could be adapted to fi + — ji~ acceleration. At lower frequencies, structures such 
as the CERN 350 MHz superconducting rf cavities could be used. These cavities ha ve b een tested in pulsed mode 



operation, and tests indicate that pulsed acceleration fields > 10 MV/m are possible [177]. 

The high single bunch intensities required for high intensities imply large higher order mode losses and large wakefield 
effects from the short, high intensity bunches. Higher order mode (HOM) load designs adapted from superconducting 
rf experience could be used. HOM loads and wakefields are expected to vary as a~ 2 and A~ 2 and <7 -1 / 2 , where a 



is the cavity aperture, A is the acceleration wavelength and a is the bunch length [17S ,17£ ] . Calculations indicate 
that the wakefields would limit bunch intensities to ss 2 x 10 12 with 1300 MHz superconducting rf in a recirculating 
linac accelerator scenario. The longitudinal dynamics is microtron-like or synchrotron-like and off-c rest acceleration 



enables compensation of the linear part of the wakefields, with synchrotron- like phase stability [168] 
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E. Simulations 



TABLE XI. Parameters of Acceleration for 4 TeV Collider 



Linac 



RLA1 



RLA2 



RCS1 



RCS2 



E (GeV) 
f r/ (MHz) 

N turns 

V r/ (GV/turn) 
C tu ™(km) 
Beam time (ms) 

0"z,6eam (cm) 
OE, beam (GeV) 

Loss (%) 



0.1-> 1.5 
30 -» 100 
1 

1.5 
0.3 
0.0013 
50 ^ 8 
0.005 -> 0.033 
5 



1.5 -> 10 
200 
9 

1.0 
0.16 
0.005 
4 -> 1.7 
0.067 -» 0.16 
7 



10 -» 70 
400 
11 
6 

1.1 

0.04 
1.7 -> 0.5 
0.16 -> 0.58 
6 



70 -> 250 
800 
33 
6.5 
2.0 
0.22 
0.5 -> 0.25 
0.58 -> 1.14 
7 



250 -> 2000 
1300 
45 
42 
11.5 
1.73 
0.25 -> 0.12 
1.14 -> 2.3 
10 



A study J 169] followed the longitudinal motion of particles thro ugh a similar sequence of recirculating accelerators 
(see table XI). Cavities similar to those proposed for TESLA [172| were assumed. Figure |3^ shows, after optimization 
of parameters, the final longitudinal phase space distributions corresponding to wakefields estimated for four different 
bunch charges: a) very small, b) 0.83 x 10 12 muons, c) 2.08 x 10 12 muons, and d) 4.17 x 10 12 muons. 

For the design beam charge of 2 x 10 12 muons (approximately as for Fig. |3^(c)) the wakcficld amplitude was 
estimated to be 2.5 MV/m, the accelerating phase was 35°, and rf voltage depression 26%. The simulation used an 
initial longitudinal phase space of 20 eV-s. It gave negligible particle loss, a final longitudinal phase space of 21.6 eV-s, 
resulting in an increase of longitudinal omittance of only 8%. 



F. Acceleration research needed 



As discussed above, possible acceleration configurations have been developed, and critical longitudinal motion 
simulations have been performed. These calculations support the general feasibility of acceleration of muons from 
cooling to collider energies. However the designs of acceleration systems have not been fully detailed and much work 
would be needed to obtain a buildable design. Complete transport lattices for linacs and return arcs have not yet been 
derived, and 6-D phase space tracking of beams through the accelerators has not been attempted. Also the geometry 
of combining and separating multi pass beams has not been worked out and optimized. 

The rf requirements and systems have been specified at only the rudimentary requirements level, and have not 
been developed to a constructible level. Optimal configurations and choices of normal or superconducting rf must be 
developed, as well as more optimal choices in acceleration frequencies. The simple wakcficld models used in the initial 
simulations should be expanded to obtain more realistic systems, and more precise calculations of wakcficld effects 
must be developed. 

Rapid accelerating systems have only been outlined at the simplest conceptual level. Prototype magnet design and 
testing are needed to test the limits of cycling rate and field strengths. Successful magnet concepts must then be 
specified in terms of stable beam transport configurations, including focusing and transport matching. While beam 
is stored for only a few turns, the individual bunch intensities are large enough that the possibility of single bunch 
instabilities must be considered and calculated. The larger number of passes in a recirculating linac accelerator places 
greater demands on the rf systems and higher order mode (HOM) loads, particularly for superconducting systems. 



VII. COLLIDER STORAGE RING 



A. Introduction 



After one fi + bunch and one fi~ bunch have been accelerated to collision energy, the two bunches are injected into 
the collider ring, which is a fixed field storage ring. Parameters for several possible collider storage rings are given in 
table H Collider ring lattices have been developed for two of the collision energies in this table: 100 GeV and 3 TeV 
in the center of mass. 

Three operational modes are proposed in the above table for the 100 GeV collider, each requiring different machine 
optics. The following sections discuss a 100 GeV collider lattice for two of the modes, the broad momentum spread 
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FIG. 39. Recirculating linac accelerator simulation results with wakefields, with beam accelerated from 200 to 2000 GeV in 
a 10-turn recirculating linac accelerator. Longitudinal phase space plots for different bunch charges: A) very small number; B) 
0.83 x 10 12 ; C) 2.08 x 10 12 and D) 4.17 x 10 12 muons in a bunch. 



case (Ap/p of 0.12%, rms) and the narrow momentum spread case {Ap/p of 0.003%), as well as a 3 TeV collider 
lattice. 



B. Collider Lattices 



1. Design criteria 



Stringent criteria have been imposed on the collider lattice designs in order to attain the specified luminosities. The 
first and most difficult criterion to satisfy is provision of an Interaction Region (IR) with extremely low (3* values at 
the collision point consistent with acceptable dynamic aperture. The required (3* values for the 100 GeV collider are 
4 cm for the broad momentum spread case and 14 cm for the narrow momentum spread case. For the 3 TeV machine, 
0* is only 3 mm. These j3* values were tailored to match the longitudinal bunch lengths in order to avoid luminosity 
dilution from the hour-glass effect. Achieving this requirement in the 3 TeV lattice is complicated by the high peak 
beta function values in the final focus quadrupoles requiring 8-10 cm radial apertures. The correspondingly weakened 
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gradients combined with the ultra-high energy make for a long final focus structure. (In contrast, the lower energy 
and larger (3* values in the 100 GeV collider lead to an efficient, compact final focus telescope.) Compounding the 
problem, particularly for the 3 TeV design, is the need to protect the superconducting coils from the decay products 
of the muons. Placing a tungsten shield between the vacuum chamber and the coils can increase the radial aperture 
in the 3 TeV quadrupoles by as much as 6 cm, lowering available gradients still further. Final focus designs must 
also include collimators and background sweep dipoles, and other provisions for protecting the magnets and detectors 
from muon decay electrons. Effective schemes have been incorporated into the current lattices. 

Another difficult constraint imposed on the lattice is that of isochronicity. A high degree of isochronicity is required 
in order to maintain the short bunch structure without excessive rf voltage. In the lattices presented here, control 
over the momentum compaction is achieved through appropriate design of the arcs. 

A final criterion especially important in the lower energy colliders is that the ring circumference be as small as 
feasible in order to minimize luminosity degradation through decay of the muons. Achieving small circumference 
requires high fields in the bending magnets as well as a compact, high dipole packing fraction design. To meet 
the small circumference demand, 8 T pole tip fields have been assumed for all superconducting magnets, with the 
exception of the 3 TeV final focus quadrupoles, whose pole tips are assumed to be as high as 12 T. In addition, design 
studies for still higher field dipoles are in progress. 



2. rf system 



The rf requirements depend on the momentum compaction of the lattice and on the parameters of the muon 
bunch. For the case of very low momentum spread, synchrotron motion is negligible and the rf system is used solely 
to correct an energy spread generated through the impedance of the machine. For the cases of higher momentum 
spreads, there are two approaches. One is to make the momentum compaction zero to high order through lattice 
design. Then the synchrotron motion can be eliminated, and the rf is again only needed to compensate the induced 
energy spread correction. Alternatively, if some momentum compaction is retained, then a more powerful rf system 
i s needed to maintain the specified short bunches. In either case, rf quadrupoles will be required to generate BNS 
[180, 181 1 damping of the transverse head-tail instability. 



3. 3 TeV CoM lattice 



The 3 TeV ring has a roughly racetrack design with two circular arcs separated by an experimental insertion on 
one side, and a utility insertion for injection, extraction, and beam scraping on the other. The experimental insertion 
includes the interaction region (IR) followed by a local chromatic correction section (CCS) and a matching section. 
The chromatic correction section is optimized to correct the ring's linear chromaticity, which is mostly generated by 
the low beta quadrupoles in the IR . In designs of e + e~ colliders, it has been found that local chromatic correction of 
the final focus is essential [ 182| - 185] , as was found to be the case here. The 3 TeV IR and CCS are displayed in Fig. |4f]. 
The accompanying 3 TeV arc module in Fig |4l] is an example of a module which controls momentum compaction (i.e. 
isochronicity) of the entire ring. 



4. 100 GeV CoM lattices 



For the 100 GeV CoM collider 186], two operating modes are contemplated: a high luminosity case with broad 
momentum acceptance to accommodate a beam with a Ap/p of ±0.12% (rms), and one with a much narrower 
momentum acceptance and lower luminosity for a beam with Ap/p of ±0.003% (rms). For the broad momentum 
acceptance case, j3* must be 4 cm and for the narrow momentum acceptance case, 14 cm. In either case, the 
bunch length must be held comparable to the value of (3* . The 100 GeV ring geometry is highly compact and more 
complicated than a racetrack, but the lattice has regions with the same functions as those of the 3 TeV ring. 

Two independent 100 GeV lattice designs have evolved; these are described below in separate sections and denoted 
Example (a) and Example (b), respectively. The first design described is a lattice which has two optics modes. In 
the high luminosity mode, the (3* value is 4 cm with a transverse and momentum aperture sufficient to accept a 
normalized beam emittance of 907r (rms) and a Ap/p of ±0.12% (rms). The second, lower luminosity mode has a 
(3* value of 14 cm with a very large transverse acceptance, but small, approximately monochromatic, momentum 
acceptance. 
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FIG. 40. Example (a): 3 TeV IR and chromatic correction. 




The second 100 GeV lattice described is another collider design with a 4 cm /3* optics mode. Although the number 
of magnets differ between the two lattices, the most important optics difference between the two is in the modules 



■ r j(i 



used in the arcs. 



5. 100 GeV CoM. Example (a) 



The need for different collision modes in the f 00 GeV machine led to an Interaction Region design with two optics 
modes: one with broad momentum acceptance (Ap/p of 0.12%, rms) and a collision (3* of 4 cm, and the other 
basically monochromatic (Ap/p of 0.003%, rms) and a larger collision j3* of 14 cm. The first lattice design, denoted 
Example (a), shown in Fig. ^ and Fig. [l3| has a total circumference of about 350 m with arc modules accounting for 
only about a quarter of the ring circumference. 
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FIG. 42. Example (a): 4 cm j3* Mode showing half of the IR, local chromatic correction, and one of three arc modules. 



The low beta function values at the IP are mainly produced by three strong superconducting quadrupolcs in the 
Final Focus Telescope (FFT) with pole tip fields of 8 T. The full interaction region is symmetric under reflection 
about the interaction point (IP). Because of significant, large angle backgrounds from muon decay, a background 
sweep dipole is in cluded in the final focus telescope and placed near the IP to protect the detector and the low (3 
quadrupolcs [187]. It was found that this sweep dipole, 2.5 m long with an 8 T field, provides sufficient background 
suppression. The first quadrupole is located 5 m away from the interaction point, and the beta functions reach a 
maximum value of 1.5 km in the final focus telescope, when the maxima of the beta functions in both planes are 
equalized. For this maximum beta value, the quadrupole apertures must be at least 11 cm in radius to accommodate 
5 <t of a 90 7r mm-mrad, 50 GeV muon beam (normalized rms emittance) plus a 2 to 3 cm thick tungsten liner [ 188 1 . 
The natural chromaticity of this interaction region is about —60. 

Local chromatic correction of the muon collide r int eraction regi on i s required to achieve broad momentum accep- 
tance. The basic approach developed by Brown [ 183 1 and others [ 189 1 , is implemented in the Chromatic Correction 



Region (CC). The CC contains two pairs of sextupoles, one pair for each transverse plane, all located at locations 
with high dispersion. The sextupoles of each pair are located at positions of equal, high beta value in the plane 
(horizontal or vertical) whose chromaticity is to be corrected, and very low beta waist in the other plane. Moreover, 
the two sextupoles of each pair are separated by a betatron phase advance of near n, and each sextupole has a phase 
separation of (2n + l)f from the IP, where n is an integer. The result of this arrangement is that the geometric 
aberrations of each sextupole is canceled by its companion while the chromaticity corrections add. 



57 




FIG. 43. Example (a): 14 cm f3* Mode showing half of the IR, local chromatic correction, and one of three arc modules. 



The sextupoles of each pair are centered about a minimum in the opposite plane ((3 m in < Ito), which provides 
chromatic correction with minimal cross correlation between the planes. A further advantage to locating the opposite 
plane's minimum at the center of the sextupole, is that this point is \ away from, or "out of phase" with, the source 
of chromatic effects in the final focus quadrupolcs; i.e the plane not being chromatically corrected is treated like the 
IP in terms of phase to eliminate a second order chromatic aberration generated by an "opposite-plane" sextupole. 

In this lattice example, the CC (Fig. |4|) was optimized to be as short as possible. The /?max is only 100 m and 
the /3 m j n = 0.7 m, giving a /? ra tio between planes of about 150, so the dynamic aperture is not compromised by a 
large amplitude dependent tuneshift. 

This large beta ratio, combined with the opposite plane phasing, allows the sextupoles for the opposite planes to 
be interleaved, without significantly increasing the nonlinearity of the lattice. In fact, interleaving improved lattice 
performance compared to that of a non-interleaved co rrect ion scheme, due to a shortening of the chromatic correction 
section, which lowers its chromaticity contribution |190fl . The use of somewhat shallower beta minima with less 
variation in beta through the sextupoles was made soften the chromatic aberrations, although this caused a slight 
violation of the exact 7r phase advance separation between sextupole partners. The retention of an exact 7r phase 
advance difference between sextupoles was found to be less important to the dynamic aperture than elimination of 
minima with /3 mm < 0.5 m. 

The total momentum compaction contributions of the IR, CC, and matching sections is about 0.04. The total length 
of these parts is 173 m, while that of the the momentum compaction correcting arc is 93 m. From these numbers, 
it follows that this arc must have a negative momentum compaction of about —0.09 in order to offset the positive 
contributions from the rest of the ring. 

The arc module is shown in Fig. ^5|. It has the small beta functions characteristic of FODO cells, yet a large, 
almost separate, variability in the mome ntum co mpaction of the module which is a characteristic associated with the 
flexible momentum compaction module [191,192]. The small beta functions are achieved through the use of a doublet 
focusing structure which produces a low beta simultaneously in both planes. At the dual minima, a strong focusing 
quadrupolc is placed to control the derivative of dispersion with little impact on the beta functions. Negative values 
of momentum compaction as low as a — —0.13 have been achieved, and jt = 2 i, has been achieved with modest 
values of the beta function. 

This arc module was able to generate the needed negative momentum compaction with beta functions of 40 m or 
less. 



58 



nfflrTTUTTti 




Dx (i 



FIG. 44. Example (a): The chromatic correction module. 
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FIG. 45. Example (a): A flexible momentum compaction arc module. 



A very preliminary calculation of the dynamic aperture [ 190 1 without optimization of the lattice or inclusion of 
errors and end effects is given in Fig. ^6|. One would expect that simply turning off the chromatic correction sextupoles 
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in the 4 cm (3* mode would result in a linear lattice with a large transverse aperture. With only linear elements, the 
4 cm (3* optics was found to be strongly nonlinear with limited on-momcntum dynamic acceptance. 

A normal form analysis using COSY INFINITY [ 193 1 showed that the variation of tune shift with amplitude was 
large, which was the source of the strong nonlinearity in the seemingly linear lattice. To locate the source of this 
nonlinearity, a lattice consisting of the original IR and arcs only (no CC), was studied. Numerical studies confirmed 
similar dynamic aperture and variation of tune shift with amplitude. This ruled out the possibility that the dynamic 
aperture was limited by the low beta points in the local chrom atic correction section and points to the IR as the 
source of the nonlinearity. These findings were also verified |194| using a Runge-Kutta integrator to track through the 
IR and a linear matrix for the rest of the lattice. Further analytical study using perturbation theory showed that the 
first order contribution to the tune shift with amplitude is proportional to 7^ and 7x7y, which are large in this IR. 
These terms come from the nonlinear terms of p x /po and Py/po, which, to first order, equal the angular divergence of 
a particle. As a demonstration, a comparison to the LHC low beta IR was done. Taking into account only the drift 
from the IP to the first quadrupole, the horizontal detuning at 10cx of the present IR (/?* = 4 cm) is 0.01, whereas 
the detuning of the entire LHC lattice is below 1E-4. This also explains the fact that the on-momentum aperture of 
the wide momentum spread mode remains roughly constant despite various versions and correction attempts. 

It was therefore concluded and later shown that the dynamic aperture of the more relaxed f3* of 14 cm would 
not have the same strong nonlincarities due to the reduced angular terms. In fact, the variation of tunc shift with 
amplitude was less by an order of magnitude; hence the large transverse acceptance shown in Fig. Bq (dashed line). 



6. 100 GeV CoM. Example (b) 

The second lattice design, Example (b), is shown in Fig. [I?] starting from the IP. The 1.5 m background clearing 
dipole is 2.5 m away from the IP and is followed by the triplet quadrupoles with the focusing quadrupole in the 
center. The interaction region (IR) stops at about 24 m from the IP. Because of the small low betatron functions 
in both transverse planes, the betatron functions at the final focusing triplets increase to ~ 1550 m. The natural 
chromaticities, of order ~ —40, are high, requiring local correction. Due to the size limitation of the collider ring, it 
appears that we have room for only two pairs of interleaved sextupoles on each side of the IP, each pair correcting 
chromaticity in one transverse plane. The correction section on each side of the IP spans a distance of roughly 61.3 m. 

The SXl's are the two horizontal correction sextupoles. They should be placed at positions with the same betatron 
functions and dispersion function, and separated horizontally and vertically by phase advances Aip x an d = ir so 
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FIG. 47. Example (b): Lattice structure of the IR including local chromaticity corrections. The maximum and minimum j3 x 
are 1571.74 and 0.040 m, the maximum and minimum /3 y are 1550.94 and 0.040 m, while the maximum and minimum dispersions 
are 4.31 and —3.50 m. The natural horizontal and vertical chromaticities are —41.46 and —39.90, giving a transition gamma 
of 7t = 5.52. The total module length is 85.32 m with a total bend angle of 1.307 rad. 



that their nonlinear effect will be confined in the region between the two sextupoles. Their horizontal phase advances 
should also be integral numbers of it from the triplet focusing F-quadrupole so that the chromaticity compensation 
for that quadrupole will be most efficient [189]. The SX2's are the two vertical correction sextupoles which should 
be placed similarly at designated locations. In general, it will be difficult to satisfy all the requirements mentioned; 
especially in this situation, luminosity arguments limit the l atti ce size. For this lattice, the Twiss properties at the 
centers of the four correction sextupoles are listed in table XI], where all the figures given by the lattice code are 
displayed. An attempt was made to satisfy all the requirements at the expense of having Aip y /(2ir) = 0.60 instead of 
0.50 for the SXl's. This trade-off is explained below. 

TABLE XII. Twiss properties of the IR correction sextupoles. 





Distance 


Phase Advances 


Betatron Functions (m) 


Dispersion 




(m) 


ipx/(2n) 


^v/(2jt) 


A* 


Py 


(m) 


SX2 


33.5061 


0.48826 


0.74953 


1.00000 


100.00012 


2.37647 


SX2 


62.3942 


0.98707 


1.24953 


1.00000 


100.00009 


2.37651 


SX1 


49.3327 


0.74892 


0.87703 


100.00023 


1.00000 


2.66039 


SX1 


74.6074 


1.24892 


1.47987 


99.99967 


0.99992 


2.65817 



The second order effects of the sextupoles contribute to the amplitude dependent tune spreads, which, if too large, 
can encompass resonances leading to dynamical aperture limitation. For example, in this lattice, 

1^ = 8.126337- 100 e x - 4140 e y , 

u y = 6.239988- 4140 e x - 54.6 e y , (28) 

where e x and e y are the horizontal and vertical unnormalized omitta nces in 7rm. In order to eliminate these tune 
spreads due to the sextupole nonlinearity, the sufficient conditions are |l95| : 
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where for the fcth thin normal sextupole with strength SNk = [B"£/(Bp)]k, 
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V+fc = (2^1^)^, and = (2i/ y ±i/ x )fc. The 5 requirements come about because there are 5 first order resonances 
driven by the sextupoles when the residual tunes of the ring satisfy [iv x ] = 0, [v±] = and two [v x ] = 0. The nominal 
tunes shown in Eq. ( p8|) are far from these resonances. Therefore, the sines in the denominators of Eq. (^9|) can be 
omitted in this discussion. Since the strengths of SX1 and SX2 are similar, we have Ssx2 <C Ssxi *C Ssx2 <C S'sxi- 
In fact, they are roughly in the ratios of 1 : (/3 max //3min) 1/2 : /3ma X //3min : (/3max//3 m in) 3/2 , which amount roughly to 
1:10:100:1000 in this lattice. In above, /3 max represents either B x at the SXl's or By at the SX2's, and B m - m represents 
either (3 y at the SXl's or /3 X at the SX2's. Thus, the first two restrictions in Eq. ( f29|) are the most important, implying 
that all /?max and /3 m ; n for each pair of SXl's must be made equal and Aip x = tt between them must be strictly 
obeyed. The third restriction is the next important one, for which S$x2 must be made equal for each pair of SX2's 
and their horizontal phase difference must equal tt. The only two parameters left are Ai/j y between a pair of SXl's and 
Aijjy between a pair of SX2's. They affect the restrictions for the v± resonances only, where the effective sextupole 
strengths S'sxi and Ssx2 are involved. Thus if we allow one restriction to be relaxed, the relaxation of Aip y = n for 
the SXl's will be least harmful. 



Flexible momentum compaction (FMC) modules 192 are used in the arc. The momentum compaction of the arc 
has to be made negative in order to cancel the positive momentum compaction of the IR, so that the whole ring 
becomes quasi-isochronous. This is accomplished in three ways: 1) removing the central dipole of the usual FMC 
module; 2) increasing the length of the first and last dipoles, and 3) increasing the negative dispersion at the entrance. 
Two such modules will be required for half of the collider ring, one of which is shown in Fig. To close the ring 
geometrically, there will be a ~ 72.0 m straight section between the two sets of FMC modules. The total length of the 
collider ring is now only C = 354.3 m. This is a nice feature, since a small ring allows a larger number of collisions 
before the muons decay appreciably. Note that the IR and local correction sections take up 48.2% of the whole ring. 
The momentum compaction factor of this ring is now ao = —2.77 x 10~ 4 . The rf voltage required to maintain a bunch 
with rms length and rms momentum spread a g is V r { = \rj\EC 2 <Tg / (2irhaf), where r\ is the slippage factor and E 
the muon energy. On the other hand, if the bucket height is taken as k times the rms momentum spread of the bunch, 
the rf harmonic is given by h = C/(kTra g ). Thus, for a ( = 4 cm and a s = 0.0012, this lattice requires an rf voltage of 
Vtf « 88fc kV. Since ao is negative already, its absolute value can be further lowered easily if needed. However, we 
must make sure that the contributions from the higher order momentum compaction are small in addition. 



The dynamical aperture of the lattice is computed by tracking particles with the code COSY INFINITY [ 193 1 . 
Initially 16 particles with the same momentum offset and having vanishing x' and y' are placed uniformly on a circle 
in the x-y plane. The largest radius that provides survival of the 16 particles in 1000 turns is defined here as the 
dynamical aperture at this momentum offset and is plotted in solid in Fig. ^ in units of the rms radius of the beam. 
(At the 4 cm low beta IP, the beam has an rms radius of 82 /xm.) As a reference, the 7 a aperture spanning ±6 sigmas 
of momentum offset is also displayed as a semi-ellipse in dashdot. To maximize the aperture, first, the tunes must 
be chosen to avoid parametric resonances. The on-momentum amplitude dependent horizontal and vertical tunes are 
given in Eq. (p8|). With the designed rms e x = e y = 0.169 x 10~ 6 7rm, the on-momentum tune variations are at most 
0.0007. Second, the chromaticity variations with momentum must be as small as possible. This is shown in Fig. ^ 
(right hand side plot). Note that there arc no families of sextupoles to correct for the higher order chromaticities 
in this small ring with only four FMC modules. As the momentum spread varies from —1 to 0.9%, v x varies from 
8.16698 to 8.07459, and v y from 6.28305 to 6.22369 for the center of the beam. 

During aperture tracking we notice that particle loss occurs mostly in the horizontal direction. We are convinced 
that the small momentum aperture is a result of the large dispersion swing in the lattice from +4.5 to —3.5 m. For 
example, 4.5 m dispersion and 0.6% momentum offset translates into a 2.7 cm off-axis motion. The nonlinearity of the 
lattice will therefore diminish the dynamical ape rtur e. A resonant strength study using, for example, swamp plots and 



normalized-resonance-basis-coefficient analysis [196] actually reveals that this lattice and some of its variations are 
unusually nonlinear. Recently, we make a modification of the FMC arc modules which have a smaller dispersion swing 
from —2.6 to +2.0 m only. The IR has not been changed except for the matching to the arc modules. The aperture 
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FIG. 48. Example (b): Lattice structure of the flexible momentum compaction module. The maximum and minimum (3 X are 
19.57 and 0.29 m, the maximum and minimum /3 y are 23.63 and 7.80 m, the maximum and minimum dispersions are 1.35 and 
—3.50 m. The natural horizontal and vertical chromaticities are —1.77 and —0.92, giving a transition gamma of jt = i4.43. 
The total module length is 27.91 m with a total bend angle of 0.917 rad. 



has been tracked with TEAPOT [ 197 1 in the same way as COSY and is plotted as dashes in Fig. ^ (left hand side 
plot). We see that the momentum aperture has widened appreciably. The dynamical aperture near on-momentum, 
however, is one sigma less than the lattice presented here. Nevertheless, it is not clear that this decrease is significant 
because all tracking has been performed in steps of one sigma only. However this type of aperture is still far from 
satisfactory, because so far we have been studying a bare lattice. The aperture will be reduced when fringe fields, 
field errors, and misalignment errors are included. 

We suspe ct th at the aperture for small momentum spread is limited by the dramatic changes in betatron functions 
near the IP 194 1. These changes are so large that Hill's equation would no longer be adequate and the exact equation 
for beam transport must be used. This equation brings in nonlinearity and limits the aperture, which can easily be 
demonstrated by turning off all the scxtupoles. In other words, although the momentum aperture can be widened by 
suitable deployment of sextupoles, the on-momentum dynamical aperture is determined by the triplet quadrupoles 
and cannot be increased significantly by the sextupoles. Some drastic changes in the low beta design may be necessary. 



C. Scraping 



It has been shown [ 198 1 that detector backgrounds originating from beam halo can exceed those from decays in 
the vicinity of the interacti on p oint (IP). Only with a dedicated beam cleaning system far enough from the IP can 
one mitigate this problem [188]. Muons injected with large momentum errors or betatron oscillations will be lost 
within the first few turns. After that, with active scraping, the beam halo generated through beam-gas scattering, 
resonances and beam-beam interactions at the IP reaches equilibrium and beam losses remain constant throughout 
the rest of the cycle. 



Two beam cleaning schemes have been designed [188], one for muon colliders at high energies, and one for those at 
low energies. 



The studies [ 188 1 showed that no absorber, ordinary or magnetized, will suffice for beam cleaning at 2 TeV; in fact, 
the disturbed muons are often lost in the IR, but a simple metal collimator was found to be satisfactory at 100 GcV. 
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FIG. 49. Example (b): Left hand side plot is dynamical aperture of the lattice vs. momentum offset. COSY calculation in 
solid, 7 a in dot-dashes, and TEAPOT calculation with modified FMC modules in dashes. Right hand side plot is chromaticities 
vs. momentum offset. 



1. Scraping for high energy collider 

At high energies, a 3 m long electrostatic deflector (Fig. [5(]) separates muons with amplitudes larger than 3 a and 
deflects them into a 3 m long Lambertson magnet, which extracts these downwards through a deflection of 17 mrad. 
A vertical septum magnet is used in the vertical scraping section instead of the Lambertson to keep the direction of 
extracted beam down. The shaving process lasts for the first few turns. To achieve practical distances and design 
apertures for the separator/Lambertson combinations, (3 functions must reach a kilometer in the 2 TeV case, but only 
100 m at 50 GeV. The complete system consists of a vertical scraping section and two horizontal ones for positive 
and negative momentum scraping (the design is symmetric about the center, so scraping is identical for both /x + and 
fx~). The system provides the scraping power of a factor of 1000; that is, for every 1000 halo muons, one remains. 



2. Scraping for low energy collider 



At 50 GeV, collimating muon halos with a 5 m long steel absorber (Fig. ^l]) in a simple compact utility section 
does an excellent job. Muons lose a significant fraction of their energy in such an absorber (8% on average) and have 
broad angular and spatial distributions. Almost all of these muons are then lost in the first 50-100 m downstream of 
the absorber with only 0.07% of the scraped muons reaching the low (3 quadrupoles in the IR, i.e. a scrap ing power 
is 1500 in this case, which is significantly better than with an earlier septum scraping system design [188] similar to 
that developed for the high energy collider. 
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FIG. 50. Schematic view of a collider beam halo extraction. 
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FIG. 51. Scraping muon beam halo with a 5 m steel absorber. 



D. Beam-beam tune shift 



Several studies have considered beam emittance grow th due to the beam-beam tune shift and none have observed 



significant luminosity loss. For instance, a study [199], using the high energy collider parameters (see table in 
which particles were tracked assuming Gaussian beam field distributions, and no muon decay, showed a luminosity 
l oss of only 4%. With muon decay included, the loss contribution from beam-beam effects is even less. Another study 
J200fl using a particle in cell approach with no assumptions about field symmetry obtained a similar result. Collisions 
between beams displaced by 10% of their radius also gave little loss. But all these studies assumed an ideal lattice, 
and none considered whether small losses due to nonlinearities give rise to an unacceptable background. 

E. Impedance/ wakefleld considerations 



A study [201 1 has examined the resistive wall impedance longitudinal instabilities in rings at several energies. At 
the higher energies and larger momentum spreads, solutions were found with small but finite momentum compaction 
and moderate rf voltages. For the special case of the Higgs Factory, with its very low momentum spread, a solution 
was found with no synchrotron motion, but rf was provided to correct the first order impedance generated momentum 
spread. The remaining off-momentum tails which might generate background could be removed by a higher harmonic 
rf correction without affecting luminosity. Solutions to the higher energy and larger momentum spread cases without 
synchrotron motion are also being considered. 

Given the very slow or nonexistent synchrotron oscillations, the transverse beam breakup instability is significant. 



This instability can be stabilized using rf quadrupolc 1 18 1 ] induced BNS damping. For in stan ce, the required tune 
shift with position in the bunch, calculated using the two particle model approximation 202 1, is only 1.58 x 10~ 4 
for the 3 TeV case using a 1 cm radius aluminum pipe. This stabilizes the resistive wall instability. However, this 
application of BNS damping to a quasi-isochronous ring, and other head-tail instabilities due to the chromaticities £ 
and rji, needs more study. 

F. Bending magnet design 

The dipole field assumed in the 100 GcV collider lattices described above was 8 T. This field can be obtained using 
1.8° niobium titanium (NbTi) cos theta superconducting magnets similar to those developed for the LHC. The only 
complication is the n eed for a tungsten shield between the beam and coils to shield the latter from beam decay 
heating. 

The /x's decay within the rings — > e~T7^v^) 1 producing electrons whose mean energy is approximately 0.35 that 
of the muons. With no shielding, the average power deposited per unit length would be about 2 kW/m in the 4 TeV 
machine, and 300 W/m in the 10 GeV Higgs factory. Figure ^2] shows the power penetrating tungsten shields of 



different thickness [44,187,188 203 1 . One sees that 3 cm in the low energy case, or 6 cm at high energy would reduce 
the power to below 10 W/m, which can reasonably be taken by superconducting magnets. 
Figure |5^ shows the cross section of a baseline magnet suitable for the 100 GcV collider. 

The quadrupoles could use warm iron poles placed as close to the beam as practical. The coils could then be either 
superconducting or warm, placed at a greater distance from the beam and shielded from it by the poles. 
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FIG. 52. Power penetrating tungsten shields vs. their thickness for a) 4 TeV, and b) 0.1 TeV, colliders. 




FIG. 53. Cross Section of a baseline dipole magnet suitable for the 100 GeV collider. 

The collider ring could be made smaller, and the luminosity increased, if higher field dipoles were used. In the 
low energy case, the gain would not be great since less than half the circumference is devoted to the arcs. For this 
reason, and to avoid yet another technical challenge, higher field magnets are not part of the baseline design of a 100 
GeV collider. But they would give a significant luminosity improvement for the higher energy colliders, and would be 
desirable there. There have been several studies of possible designs, three of which (two that are promising and one 
that appears not to work) are included below. 
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1. Alternative racetrack NbzSn dipole 



A higher field magnet based on Nl^Sn conductor and racetrack coils is presently being designed. The NbaSn 
conductor allows higher fields and provides a large temperature margin over the operating temperature, but being 
brittle and sensitive to bending or other stress, presents a number of engineering challenges. 

In this design, the stress levels in the conductor are reduced by the use of a rectangular coil block geometry and 
end support problems are reduced by keeping the coils flat. In the more conventional cos theta designs, the conductor 
is distributed around a cylinder and the forces add up towards the midplanc; in addition, the ends, as they arc over 
the cylinder, are relatively hard to support. 

The geometry of the cross section is shown in Fig. [54|. It uses all 2-D flat racetrack coils. Each quadrant of the 
magnet aperture has two blocks of conductors. The block at the pole in the first quadrant has a return block in the 
second quadrant, similar to that in a conventional design. The height of this block is such that it completely clears 
the bore. In a conventional design, the second block, the midplane block, would also have a return block in the second 
quadrant. That would, however, require the conductor block to be lifted up in the ends to clear the bore and thus 
would lose the simple 2-D geometry. In the proposed design, the return block retains the 2-D coil geometry, as it is 
returned on the same side (see Fig. |54|) and naturally clears the bore. Since the return block does not contribute to 
the field, this design uses 50% more conductor. This, however, is a small penalty to pay for a few magnets where the 
performance and not the cost is a major issue. The field lines ar e also shown in Fig. ^4]. 

Preliminary design parameters for two cases are given in table Kill . The first case is one where the performance of 
the cable used is the same that is in the LBL D20 magnet, which created a central field of 13.5 T. The second case is 
the one where the cable is graded and two types of cable are used, and it is assumed that a reported improvement in 
cable performance is realized. It is expected to produce a central field of 14.7 T when operated at 4.2 °K. 



TABLE XIII. Preliminary design parameters for a racetrack NbsSn dipole with two different types of cable. 
Case 1 : Same conductor as in LBL 13.5 T D20 magnet without grading 



Central field at quench 


13 T at 4.2 °K 


Coil dimensions 


25 mm x 70 mm 


Total number of racetrack coils in whole magnet 


6 


Total number of blocks per quadrant in aperture 


2 (+1 outside the aperture) 


Yoke outer radius 


500 mm (same as in D20) 


Field harmonics 


a few parts in 10 at 10 mm 


Midplane gap (midplane to coil) 


5 mm (coil to coil 10 mm) 


Minimum coil height in the end 


45 mm (Note: coils are not lifted up.) 


Case 2: Newer conductor and graded 


Central field at quench 


14.7 T at 4.2 °K 


Grading 


70 mm divided in two 35 mm layers 


Overall current densities 


370 A/mm 2 and 600 A/mm 2 


Peak fields 


16 T and 12.5 T 


Copper current density 


1500 A/mm 2 


Other features are the same as in Case 1 





2. Alternative Cos Theta Nb^Sn dipole 



In this case the problem with the brittle and sensitive conductor is solved by winding the coil inside many separate 
slots cut in metal support cylinders. There is no build up of forces on the coil at the mid-plane. The slots continue 
around the ends, and thus solve the support problem there too. 

Figure pq sho ws th is alternative NbaSn dipole cos theta design. It is an extension of the concept used to build 



helical magnets [204] for the polarized proton program at RHIC [ 205 1 . The magnet is wound with pre- reacted 



kapton-insulatcd, B-stage impregnated, low current cable. The build up of forces is controlled by laying the cables 
in machined slots in a metal support cylinder. After winding, the openings of the slots are bridged by metal spacers 
and the coils pre-comprcsscd inward by winding B-stagc impregnated high tensile thread around the spacers. After 
curing, the outside of each coil assembly is machined prior to its insertion into an outer coil, or into the yoke. There 
are 3 layers. The inner bore is 55 mm radius, the outer coil radius approximately 118 mm, and the yoke inside radius 
is 127 mm. The maximum copper current density is 1300 A/mm 2 . 
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FIG. 54. Cross section of alternative high field (sa 15 T) race track coil dipole magnet with Nt>3Sn conductor. 




FIG. 55. Cross section of alternative high field slot dipole made with NbaSn conductor. 



Using the same material specifications as used in the above high field option, a central short sample field of 13.2 
T was calculated. This is somewhat less than the block design discussed above, but could be improved by increasing 
the cable diameters to improve the currently rather poor (64%) cable to cable-plus-insulator ratio. 
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3. Study of C-magnet dipole 
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FIG. 56. Cross section of an unsuccessful alternative high field C magnet with open mid-plane. 



Figure [56] shows the cross section of a high field dipole magnet in which it was hoped to bring the coils closer to 
the beam pipe without suffering excessive heating from beam decay. The coil design [206] appeared reasonable, but 
the required avoidance of coil heating was not achieved. 

Decay electrons are generated at very small angles (« I/7) to the beam, and with an average energy about 1/3 
of the beam. Such electrons initially spiral inward (to the right in Fig. |5^) bent by the high dipole field. In the high 
energy case, these electrons also radiate a significant fraction of their energy as (as 1 GcV) synchrotron gamma rays, 
some of which end up on the outside (to the left in Fig. ^). The concept was to use a very wide beam pipe, allow 
the electrons to exit between the coils, and be absorbed in an external cooled dump. Unfortunately a preliminary 
study found that a substantial fraction of the electrons did not reach the dump. They were bent back outward before 
reaching it by the return field of the magnet coils and the nature of the curved ring geometry. Such electrons were 
then trapped about the null in the vertical field and eventually hit the upper or lower face of the unshielded vacuum 
pipe. They showered, and deposited unacceptable levels of heat in the coils. 

Another idea called for collimators between each bending magnet that would catch such trapped electrons. This 
option has not been studied in detail, but the impedance consequences of such periodic collimators are expected to 
be unacceptable. 

Further study of such options might find a solution, but the use of a thick cylindrical heavy metal shield appears 
practical, adequate, and is thus the current baseline choice. 



G. Energy scale calibration 



In order to scan the width of a Higgs boson of mass around 100 GeV, one needs to measure the energy of the 
individual muon stores to an accuracy of a few parts per million, since the width of a Higgs boson of that mass is 
expected to be a few MeV. Assuming that muon bunches can be produced with modest polarizations of ks 0.25, and 
that the polarization can be maintained from turn to turn in the collider, it is possible to use the precession of the 
polarization in the ring to measure accurately the average energy of the muons fl2| | . The total energy of electrons 
produced by muon decay observed in the calorimeter placed in the ring varies from turn to turn due to the g — 2 
precession of the muon spin, which is proportional to the Lorentz factor 7 of the muon beam. Figure [5^ shows the 
result of a fit of the total electron energy observed in a calorimeter to a functional form that includes muon decay 
and spin precession. Figure ^8| shows the fractional error £7/7 obtained from a series of such fits plotted against the 
fractional error of measurement in the total electron energy that depends on the electron statistics. It has been shown 
that precisions of a few parts per million in 7 are possible with modest electron statistics of sw 100, 000 detected. It 
should be noted that there are 3.2 x 10 6 decays per meter for a muon intensity of 10 12 muons. 

Our current plans to measure the energy due to decay electrons entail an electromagnetic calorimeter that is 
segmented both longitudinally and transversely and placed inside an enlarged beam pipe in one of the straight 
sections in the collider ring. The length of the straight section upstream of the calorimeter can be chosen to control 
the total number of decays and hence the rate of energy deposition. The sensitive material can be gaseous, since the 
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FIG. 57. a) Energy detected in the calorimeter during the first 50 turns in a 50 GeV muon storage ring (points). An average 
polarization value of P = —0.26 is assumed and a fractional fluctuation of 5 x 10~ 3 per point. The curve is the result of a 
MINUIT fit to the expected functional form, b) The same fit, with the function being plotted only at integer turn values. A 
beat is evident, c) Pulls as a function of turn number, d) Histogram of pulls. A pull is defined by (measured value-fitted 
value) /(error in measured-fitted). 



energy resolution is controlled by decay fluctuations rather than sampling error. In order to measure the total number 
of electrons entering the calorimeter, we plan to include a calorimeter layer with little absorber upstream of it as the 
first layer. 

This scheme will enable us to calibrate and correct the energy of individual bunches of muons and permit us to 
measure the width of a low mass Higgs boson. 



VIII. RADIATION AND BACKGROUNDS 



A. Conventional radiation 

The proton source generates a 4 MW proton beam, which is comparable to the proposed spallation source [P2[ . 
This is a very high power and will, as in the spallation source, require great care in reducing unwanted particle losses, 
as well as careful machine shielding, and target and beam dump design. Initial studies of the target and capture 
solenoid region have been performed with the MARS code, and preliminary specifications for shielding determined, 
but more work is needed. 

The cooling and accelerator chain is rather clean, since a relatively small fraction of the muons decay, and their 
energies are low. Power deposited in the accelerators is typically 10-30 W/m (see table [X and table |x|). 



If no muons are lost, then the only sources of radiation are the muon decays yielding electrons and neutrinos. The 
neutrino radiation we discuss below. The electrons shower in the collider beam pipe shields, depositing most of their 
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FIG. 58. a) Fractional error in ^7/7 obtained from the oscillations as a function of polarization P and the fractional error 
in the measurements PERR. b) Fractional error in £7/7 obtained from the decay term as a function of polarization P and the 
fractional error in the measurements PERR. c) The total \ 2 of the fits for 1000 degrees of freedom. PERR is the percentage 
measurement error on the total electron energy in the calorimeter measuring the decay electrons. 



energy there and a relatively small amount in the magnet coils and yoke. Rad ioactivation levels, as calculated by 



MARS [207 1, after five years of 4 TeV collider operation are given in table XIV for the cases immediately after turn 
off and 1 day after turn off. It is seen that the areas in the tunnel that are outside the magnets are relatively free 
of radioactivation. Special procedures will be needed when the shield pipe has to be opened, as for instance when a 
magnet is changed. For the lower energy colliders, the radioactivation levels are proportionally less. 

TABLE XIV. 4 TeV (CoM) collider ring radioactivation levels (mrem/hour) after turn off, for parameters in table | 





immediate 


after 1 day 


Inside face of shield 


9000 


4000 


Outside face of shield 


200 


170 


Outside of coils 


30 


14 


Outside of yoke 


3 


1.4 



If muons are lost either accidentally, by scraping, or deliberately after some number of turns, then the muons 
penetrate to considerable distances in the soil/rock (3.5 km at 2 TeV, 800 m at 250 GeV) and deposit their energy 
directly or through their interaction products. Figure^ and Fig. |6^ show the distribution of radiation levels, assuming 
25% of all muons (4 bunches of 2 x 10 12 at 15 Hz) are dumped into soil/rock with density 2.24 g/cm 3 . The outer 
contours correspond to the federal limits, reaching at maxima of 18 m (2 TeV) and 14.5 m (250 GeV). To confine this 
radiation beneath the ground one can deflect the extracted beams down by 4.5 mrad at 2 TeV and about 10 mrad at 
<250 GeV. If any water were present in the soil/rock, then the first two meters around the tunnel and around the 
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aborted beam axis would require insulation or drainage up to a distance of 2.5 km at 2 TeV or 550 m at 250 GeV. 

-l 




FIG. 59. Isodose contours in the soil/rock (p=2.24g/cm 3 ) for 2 TeV muons extracted at 3 x I0 13 per second. Right scale is 
dose rate in rem/s. 




Depth (cm) 

FIG. 60. Isodose contours in the soil/rock (p=2.24g/cm 3 ) for 250 GeV muons extracted at 3 x 10 13 per second. Right scale 
is dose rate in rem/s. 
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B. Neutrino induced radiation 



It has been shown [63 , 201 211] that the neutrinos created in muon beam decays can generate excessive secondary 



radiation a t la rge distances from a muon collider (see Fig. [n]). The surface radiation dose Db{Sv) in units of 



equivalent [ 2 1 2 1 doses (Sv) over a time t(s), in the plane of a bending magnet of field B(T), in a circular collider with 
beam energy E(TeV), average bending field < B(T) >, at a depth d(m) (assuming a spherical earth), with muon 
current (of each sign) of /^(muons/s/sign) is given by: 

Db „ 4.4 x 10- 24 <*> (31) 

a B 

and the dose D$ at a location on the surface, in line with a high beta straight section of length £(m), is: 

D s » 6.7 x lO" 24 I < B >. (32) 

The equation for D s assumes that the average divergence angles satisfy the condition: ag « -. This condition is 
not satisfied in the straight sections approaching the IP, and these regions, despite their length, do not contribute a 
significant dose. 

diak width = 4 m at 30 km 
from the collider 




FIG. 61. Neutrino radiation disk. For a 3 TeV CoM collider the neutrino radiation width is « 4 m at a distance of 30 km. A 
hot spot produced by 0.1 m straight section in the ring contains roughly twice the number of neutrinos on the disk on average, 
depending on the details of the collider lattice. 

For the 3 TeV parameters given in table | and muon currents 1^ = 6 x 10 20 fi~ /yr, < B >= 6 T, B = 10 T and 
depth = 500 m, and taking the Federal limit on off-site radiation dose/year, -Dp e( j, to be 1 mSv/year (100 mrem/year), 
the annual dose Db (1 year is defined as 10 7 s), in the plane of a bending dipole is, 



D B = 1.07 x 10~ 5 Sv w 1% D 



and for a straight section of length 0.6 m is: 

D s = 9.7 x 10~ 5 Sv 



Fed' 



10% D Fcd , 



(33) 
(34) 



which may be take n to be within a reasonable limit. The general trend of these expressions has been verified by Monte 
Carlo simulations [ 213 1 using MARS. In particular, for the 3 TeV case the needed depth to stay within 1% £>p ec j, is 
300 m instead of 500 m. 

Special care will be required in the lattice design to assure that no field- free region longer than 0.6 m is present. 
This may sound difficult, but it may be noted that the presence of a field of even 1 T, is enough to reduce the d ose to 
a level below the Federal limit. The application of such a field over all rf and other components seems possible [ 2 13 1 . 

For lower energy machines, the requirements rapidly get easier: a 0.5 TeV machine at 100 m depth could have 25 m 
long sections, for the same surface dose. For a 100 GeV machine the doses are negligible. 

For machines above 3 TeV, various strategies can be employed: 
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• The machines could be built at greater depths (mines many km deep are common). 

• The vertical beam orbits in the machine could be varied so as to spread the plane of radiation and thus reduce 
the peak doses. 

• The specific locations in line with straight sections could be purchased and restricted. 

• Straight sections could be shortened further by using continuous combined function magnets. 

• The machines could be built on an island, but this could have difficulties associated with access to power and 
other utilities. 

But for any large increase in energy, to 10 TeV for instance, some reduction in muon beam flux probably will be 
required. The resultant loss of luminosity might be made up in a number of ways H : 



The beam-beam tune shi ft co nstraint could be avoided by introducing a conducting medium (e.g. liquid lithium) 



at the interaction point [214] 



• The focusing strength could be increased by the use of plasma or other exotic focusing method. 



Better cooling could be developed. Optical stochastic cooling [137|, for instance, might reduce the emittances 
by many orders of magnitude, thus greatly reducing the required beam currents. Indeed, such cooling would 
require lower currents to function appropriately. 

Such options will need future study. 

C. Muon decay background 

With 4 x 10 12 muons per bunch in a 2 + 2 TeV collider ring there are approximately 4 x 10 5 muon decays per meter 
giving rise to high energy electrons. These off-energy, off-axis electrons undergo brcmsstrahlung when they traverse 
magnetic fields. When they exit the beam pipe they interact and produce electromagnetic showers and, to a lesser 
extent, hadrons and muons. Much of this debris can be locally shielded, so the primary concern is muon decays near 
the interaction point p4j . This is the background we discuss in some detail below. 

Detailed Monte Carlo simulations of electromagnetic, hadronic and muon components of the background 
20^j213|j215| , plo| have been performed using the MARS (Tl| and GEANT |L4§ codes. The most re- 



[44,187 198 



cent study 21!:] has been done with GEANT. Figure p2| shows the final 130 meters of the 2 + 2 TeV detector region 
in this study. It includes the final four quadrupoles, dipoles and a solenoidal field surrounding the detector. This 
study: 

• followed shower neutrons and photons down to 40 keV and electrons to 25 keV. 

• used a tungsten shield over the beam, extending outward to an angle of 20 degrees from the axis. 

• Inside this shield, the clear radius has a minimum, in the high energy cases, at a distance from the IP of 1.1 m 
(80 cm for 50+50 GcV). At this point, and in an expanding cone beyond it, the clear radius is maintained at 
approximately 4 sigma of the beam size. 

• Between this minimum aperture point and the IP, the clear radius follows an inverse cone, increasing as it 
approaches the IP, with an angle a little greater than the 4 sigma of the beam divergence. These cones are 
designed so that the detector could not 'see' any surface directly illuminated by the initial decay electrons, 
whether in the forward or backward (albedo) direction (see figure |63| ). 

• The resulting open space between the IP and the tip of the cone is approximately 3 cm in the 4 TeV and 
500 GeV CoM cases, and approximately 6 cm in the 100 GeV CoM case. 

• The inner surface of each shield is shaped into a series of collimating steps and slopes to maximize the absorption 
of electron showers from electrons at very small angles to the cone surface, thus reducing the funnelling of low 
energy electrons down the pipes. 

• Further upstream, prior to the first quadrupole (from 2.5 to 4 m in the Higgs case), an 8 T dipole, with 
collimators inside, is used to sweep decay electrons before the final collimation. 
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FIG. 62. Region up to 130 m from the IP, of the 2 + 2 TeV interaction region modeled in GEANT. The triangular blue 
regions represent tungsten shielding. On the right hand of the figure, the red areas represent quadrupoles in the beam line. The 
areas around the IP represent the various detector volumes used in the calculations of particle fiuences. The detector (white 
and green areas) is 10 m in diameter and 20 m long. 



Note that there is currently an inconsistency, in the very low Ap/p Higgs Factory Case, between the short open 
space between shields (+/- 6 cm) and the rms source length (^source = VVT^CTz) of 10 cm. Some modifications to 
the parameters and shielding design will be required for this case. 

Every modern detector will have to be able to identify and reconstruct secondary vertices such as those associated 
with b-quark decays. In order to estimate the viability of a vertex detector we have to show that the occupancy of 
its elements is not higher than about 1%. Figure [34] shows the occupancy as a function of radial distance from the 
interaction point for the three CoM energies studied : 0.1, 0.5 and 4 TeV. The occupancy was calculated for silicon 
pads of 300 fj.ni x 300 fim, and assuming interaction probabilities of 0.003 and 0.0003 for low energy photons and 
neutrons respectively. One can observe that the total occupancy (left figure) is above one percent for small radii. Most 
of the hits is due to conversions of photons. The occupancy due to hits resulting from charged particles is below 1% 
(right hand figure). One can lower the occupancy at small radii by using smaller pixel sizes, as indicated in table XV 
below, a s we ll as by using innovative detector ideas as described in the next section. 

Table XV gives the hit density for the Higgs factory from the various sources and the occupancy of pixels of the 
given sizes; in each case the number is given per bunch crossing. The hit density for the higher energy machines is 
found to be somewhat lower due to the smaller decay angles of the electrons. 



TABLE XV. Detector backgrounds from /x decay 



Radius 


cm 


5 


10 


20 


100 


Photons hits 


cm" 2 


26 


6.6 


1.6 


0.06 


Neutrons hits 


cm" 2 


0.06 


0.08 


0.2 


0.04 


Charged hits 


cm" 2 


8 


1.2 


0.2 


0.01 


Total hits 


cm" 2 


34 


8 


2 


0.12 


Pixel size 


[im x [im 


60 x 150 


60 x 150 


300 x 300 


300 x 300 


Total occupancy 


% 


0.6 


0.14 


0.4 


0.02 


Occupancy charged 


% 


0.14 


0.02 


0.04 


0.002 
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FIG. 63. Detail of the tungsten shielding designed for the 50 + 50 GeV case. It is designed so that the detector is not 
connected by a straight line with any surface hit by decay electrons in forward or backward directions. The picture extends 
out to a radius of 6 cm and, on the right, to a distance 4 m from the IP. The dipole from 2.5-4.0 m is not shown. 



The radiation damage by the neutrons on a silicon detector has also been estimated. In the Higgs case, at 5 cm from 
the vertex, the number of hits from neutrons above 100 keV is found to be 1.8 x 10 13 per year. This is significantly less 
than that expected at the LHC which is now ordering silicon detectors claimed to survive 5 x 10 14 hits, approximately 
three times that assumed here. The damage for silicon detectors in the higher energy machines is of the same order 
(see table [XV]). 



TABLE XVI. Radiation damage by neutrons on silicon detectors. The working assumptions are: 1000 turns, 15 Hz and 
1 year=10 7 s. An acceptable number of hits per year is 1.5 x 10 14 . 



CoM 
(TeV) 
4 

0.5 
0.1 



/i's/bunch 
(10 12 ) 
2 
4 
1 



neutrons / cm 2 / crossing 
(above 100 KeV) 
100 
50 
30 



Hits/year 
(10 13 ) 
3 
3 

1.8 



Lifetime 
(years) 
5 
5 
8 



D. Halo background 

Muon halo refers to those muons which are lost from the beam bunch as it circulates around the collider ring. In 
conventional electron or proton accelerators, beam particles which are lost away from the IP are of little concern as 
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FIG. 64. Occupancy for 300 /imx 300 ^im silicon pads, as a function of the radius for the three energies studied. Left figure 
shows the total occupancy and the right figure shows the occupancy from hits resulting from charged particles. 



they can be locally shielded. However, muons can traverse long distances and therefore have the potential to generate 
background in a detector. The magnitude of this background depends on a detailed knowledge of the injected beam 
profile and a credible model for beam halo and beam losses. More work is needed before these are well enough 
understood. Nevertheless, it is clear that the beam will need careful preparation before injection into the collider, 
and the injection system will have to be precise and free of ripple. 

The collimation system described in the previous subsection was designed to scrape the beam both initially and 
during the 1000 turns, to assure that all loss occurs at t he s craper and not near the IP. That study indicated 
suppressions better than 10 3 of background in the detector [188], 



Beam loss must be limited as far as possible. Gas scattering has been studied |217| and shown to give a negligible 
contribution. The effects of beam-beam scattering are under study and need further work. Momentum spread tails 
from uncorrected wakcficld effects must be controlled. Assuming that the total loss from all causes, after injection 
and the first few turns is less than 10~ 4 in 1000 turns, (i.e. 10 -7 per turn), then the number of background muons 
passing through the detector should be less than 800 (2 x 4 x 10 12 x 10 -7 x 10~ 3 ) per turn. This is a low density of 
tracks per cm 2 and should be acceptable, but lower losses or better scraping would be desirable. 



E. Pair production 



Coherent beam-beam electron pair production (beamstrahlung) has been shown [ 200, 218| to be negligible, but the 
incoherent pair production (i.e. e + e ~) in the 4 TeV collider case is significant. 

The cross section is estimated to be 10 mb [ 2 IS | , which would give rise to a background of w 3 x 10 4 electron pairs 



per bunch crossing. The electrons at production do not have significant transverse momentum but the fields of the 
on-coming 3^m bunch can deflect them towards the detector. A simple program was written to track electrons from 
close to the axis (the worst case) as they are deflected away from the bunch center. Once clear of the opposing bunch 
the tracks spiral under the influence of the experimental solenoid field. Figures |6q shows the radii vs. length of these 
electron tracks for initial momenta from 3.8 to 3000 MeV in geometric steps of v2. Fig. |65|(a) is for a solenoidal field 
of 2 T and Fig. |6|(b) for 4 T. In the 2 T case tracks with initial energy below 30 MeV do not make it out to a detector 
at 10 cm, while those above 100 MeV have too small an initial angle and remain within the shield. Approximately 
10% (3000 tracks) of these are in this energy range and pass through a detector at 10 cm. The track fluences at the 
ends of the detector are less than 10 tracks per cm 2 which should not present a serious problem. At 5 cm, there are 
4500 tracks giving a fluence of 30 per cm 2 , which is also probably acceptable. If the detector solenoid field is raised 
to 4 T, then no electrons reach 10 cm and the flux at 5 cm is reduced by a factor of 2. 
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F. Bethe-Heitler muons 



TABLE XVII. Bethe-Heitler Muons 



CoM Collider Energy (TeV) 


4 


0.5 


0.1 


Assumed source length (m) 


130 


33 


20 


fi [pmuon > 1 GeV/c) per electron 


5.4 x 1(T 4 


8.3 x 10~ 5 


9.6 x 10 -6 


Beam ^i's per bunch 


2 x 10 12 


2 x 10 12 


4 x 10 12 


Bethe-Heitler fi'a per bunch crossing (xlO 3 ) 


28 


17.5 


6.1 


< Pmuon > initial (GeV) 


22 


9.5 


4.4 


/i's entering calorimeter 


220 


160 


25 


< Pmuon > (GeV) 


15.4 


6.3 


1.8 


< E dEp > (GeV) 


2.9 


1.3 


0.4 


Total E dep (GeV) 


640 


210 


10 


Edep pedestal subtracted (GeV) 


50 


2o 


1 


Fluctuation in Edep (GeV) 


55 


15 


1 


Etrans pedestal subtracted (GeV) 


15 


15 


.5 


Fluctuation in Etrans (GeV) 


40 


8 


0.5 



The GEANT/MARS studies JM, 203, 207 also found a significant flux of muons with quite high energies, from fi 
pair production in electromagnetic showers (Bethe-Heitler). Figures ^36| and ^ show the trajectories of typical muons 
from their sources in the shielding around the beam pipe to the detector. Figure |6^ is for a 4 TeV CoM collider, 
where the muons have high energy and long path lengths. A relatively long (130 m) section of beam pipe prior to 
the detector is shown. Figure |(57] is for the 100 GeV CoM collider for which, since the muons have rather short path 
lengths, only a limited length of beam pipe is shown. Note that the scales are extremely distorted: the 20° shielding 
cones on the right hand of the figures appear at steeper angles. 

The most serious effect appears to arise when these muons make deeply inelastic interactions and deposit spikes of 
energy in the electromagnetic and hadronic calorimeters. This is not serious in the Higgs case, for which the fluxes and 
cross sections are low, but at the higher collider energies they generate significant fluctuations in global parameters, 
such as transv erse energy and missing transverse energy. 

Table XVII gives some parameters of the muons for three different machine energies. In the 4 TeV and 500 GeV 
CoM cases, massive lead shielding outside the focus quadrupoles has been included. 

Figures ^ and ^ show energy deposition from Bethe-Heitler muons in a typical bunch crossing. These depositions 
are plotted against the cosine of the polar angle and azimuthal angle in the calorimeter for 4 TeV and for 500 GeV 
CoM, respectively. The massive lead shielding referred to above was not included in this study. Right hand plots in 
Figs, and p9 show the same distributions with a 1 ns timing cut. It is seen that the timing cut, if it is possible, 
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FIG. 66. Trajectories of typical Bethe-Heitler muons from their source in the shielding around the beam pipe to the detector 
for a 4 TeV CoM collider. As indicated in the text the scales are extremely distorted, the total horizontal lengt h is w 130 m 
and the outer edge of the calorimeter is « 4 m. Notice that < 1% of the tracks end in the calorimeter (see table KVII ). 



is effective in removing energy spikes at small rapidity, but has little effect in the forward and backward directions. 
The overall reduction in energy deposition is about a factor of two. 

The energy spikes can cause at least three problems: 1) they affect the triggers and event selections based on overall 
or transverse energy balance; 2) they can generate false jets and 3) they can give errors in the energies of real jets. 
After a pedestal subtraction, the effects on energy balances do not seem serious. The generation of false jets can be 
eliminated by a longitudinal energy distribution cut without introducing significant inefficiency. Energy errors in real 
jets appear to be the most serious problem. They can be reduced by the application of radial energy distribution 
cuts, but such cut s int roduce significant inefficiencies for lower energy jets. More study is needed. 

Earlier studies |216| with MARS, using less sophisticated shielding, gave results qualitatively in agreement with 
those from GEANT. 



IX. DETECTOR SCENARIOS 



The background consists of neutral and charged particles. For neutrons, the longitudinal and radial flucnecs were 
found to be comparable. The photons (average energy about 1 MeV) show a clear radial source. The charged particles 
and the photons do not all point back to the interaction point, but to the general vicinity of the IP, namely to the 
region where the 20 degree tungsten shield becomes thinner. The flux of secondary muons (Bcthe-Hcitlcr pairs) is 
mainly longitudinal. 

We would expect this background to pepper the tracking volume with random hits and produce significant energy 
pedestals in the calorimeter cells. These effects are considered in more detail in the following sections. In general, 
in designing a strawman detector that must operate in a large background flux we will want to employ as many 
detector channels as is practical. A strawman muon collider detector design with a few times 10 6 non-pixel channels 



would seem reasonable [219]. Over the last few years, development of pixel detec tors has resulted in a quantum jump 
in the number of electronic channels. For example, the SLD vertex detector p2Q| contains 300 x 10 6 pixels, and 



similar numbers of pixels are planned for the LHC vertex detectors. Hence, a strawman muon collider vertex detector 
employing 10 8 — 10 9 pixels would seem reasonable. 
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FIG. 67. Trajectories of typical Bethe-Heitler muons from their source in the shielding around the beam pipe to the detector 
for a 100 GeV CoM collider. As indicated in the text the scale is extremely distorted, the total horizontal lengt h is w 20 m 
and the outer edge of the calorimeter is « 4 m. Notice that < 0.5% of the tracks end in the calorimeter (see table XVII). 



A. Silicon vertex detector schemes 



From tabic XVI, it can be seen that the radiation damage to silicon de tect ors is acceptable in terms of the number 
of hits per year and the resultant lifetime of the detector. This prompts ]219(] us to consider the following options for 
silicon vertex detector design for the muon collider: 

• Silicon drift detector. The idea, which is described in the muon collider feasibility study G3IJ33]], is to exploit 



HI 

(see 



the time gap between bunch crossings by using the silicon drift detector technology [221] (see Fig. |7C| ). Using 
50 x 300 [ivc? detectors it should be possible to obtain a resolution of a few microns in the drift direction. This 
would facilitate a very precise vertex detector, although questions of radiation hardness remain to be resolved 
for this option. 



• Columnar pixels |222[| . The idea is to exploit the very well localized primary vertex position by using long thin 
tracking pixels that point at the IP and therefore record large ionization signals only for tracks coming from the 
IP (Fig. |7l|). For example, one can construct 50 x 50 /im 2 pixels that are 300 /zm deep. The pixels are produced 
using controlled feed-through-drilling technology to create a lattice of anodes and cathodes that extend through 
the 300 /im thick wafer. 

• Pixel micro-telescopes [ [223[ . The idea is to replace a single pixel layer with two layers separated by a small 
distance, and read them out by taking the AND between appropriate pairs. The distance between the layers 
is optimized so that soft MeV tracks (which are associated with almost 80% of the predicted background hits) 
produced in one layer curl up in the magnetic field before reaching the second layer. Thus, the pixel micro- 
telescope is blind to the soft background hits and also blind to tracks that do not come from the IP. In the 
example shown in Fig. |7^ the top measurement layer has a finer granularity than the bottom confirmation layer. 
The corresponding rows in the two pixel layers can be read out with different clock speeds to maintain the 
correct correspondence at the input into the AND gate that registers valid hits in the telescope. If the readout 
rows are the ones parallel to the beam direction, then variable clock speeds can be used to maintain the correct 
accepted direction with respect to the IP. 
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The challenge of a high background environment is clearly fruitful ground for new ideas. The above considerations 
suggest that, provided silicon detectors can be used in the inner tracking volume, it should be possible to construct a 
vertex detector able to tag secondary vertices from short lived particles at a muon collider. Detailed simulations are 
currently underway to establish this more concretely. 



B. Outer tracking schemes 

The predicted background fluxes for a Higgs factory detector at a radius of 50 cm are 200 photons/cm 2 , 
350 neutrons/cm 2 , and 0.08 charged tracks per cm 2 . The neutron flux is therefore about the same as the flux in 
the inner tracking volume, whereas the photon and charged particle fluxes are significantly less than those predicted 
at smaller radii. There are two alternative tracking strategies under consideration: 



Low field, large tracking volume drift chamber option. This option, which is described in the muon collider book 



[224 1, uses a TPC to exploit the 20 fjs time between bunch crossings. This option is viable for the very high 
energy muon collider (1.5 x 1.5 TcV). The large neutron flux necessitates choosing a gas that does not contain 
hydrogen. A mixture of 90% neon plus 10% CF4 gives a drift velocity of 9.4 cm//zs, which is close to that required 
to match the bunch crossing time. High-py tracks from the IP embedded in the predicted background flux have 
been simulated for the TPC shown in Fig. [73|. The simulation includes ionization, drift and diffusion of the 
electrons in the gas, multiplication, and other details of the detection process. The majority of the background 
hits arises from low energy Compton recoils yielding very low energy electrons that have a radius of curvature 
of less than 1 mm in the 2 T field. Their projection on the readout plane covers not more than one readout 
pitch (0.3 x 0.4 cm 2 ). These background electrons, together with the nuclear recoils from neutron scatters, yield 
large pulses that can be removed by cutting on the maximum acceptable pulse height. The simulation predicts 
that with an average background flux of 100 photons/cm 2 , reasonable pulse height cuts remove only 1% of the 
effective TPC volume, and yield tracks of high quality. However, it was realized that positive ion build-up may 
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FIG. 69. Left hand plot shows the energy deposition from Bethe-Heitler muons vs. the cosine of the polar angle and azimuthal 
angle in the calorimeter for a 0.5 TeV CoM collider. Right hand plot shows the same distributions with a 1 ns timing cut. 
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FIG. 70. Silicon drift vertex detector. 

be a problem with the design shown in Fig. |73|. If this problem can be overcome, the design shown in the figure 
yields a simulated momentum resolution of about 1.2% for tracks with p^ = 50 GeV/c. 

• High field, compact silicon tracker option. An alternative strategy is to make a compact tracker by using 
silicon in a high field (for example, 4 T). As an example, consider the geometry shown in Fig. [74| in which a 
4-layer pixel vertex detector is embedded in a 4-layer small angle stereo cylindrical silicon microstrip detector 
with a 50 x 300 /jm 2 resolution. We take the inner layer of the vertex detector to consist of a cylinder of 
50 x 300 /zm 2 pixels, and the outer 3 vertex layers to consist of spherical shells of 50 x 50 fim 2 columnar 
pixels or pixel micro-telescopes. The system is assumed to correspond to 15% of a radiation length at 90°. 
Using a parametric calculation of the momentum resolution, including multiple scattering, we obtain a p /p 2 = 
10~ 4 (10~ 2 ) (GeV/c)" 1 for p = 100 GeV/c (1 GeV/c). 
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FIG. 71. Columnar pixel geometry. Courtesy of A. Sill. 
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FIG. 72. Pixel micro-telescope geometry [223], showing trajectories of 0.2 GeV/c, 0.5 GeV/c, and 1 GeV/c tracks coming 
from the IP and bending in a 4 T field. 



Both the low field and high field tracking solutions look interesting and should be pursued with more complete 
simulations. 



C. Electromagnetic calorimeter schemes 



Background particles entering the electromagnetic calorimeter are expected to give rise to significant energy 
pedestals in the calorimeter cells. Consider a 4 m long calorimeter that is 25 radiation lengths deep, has an in- 
ner radius of 120 cm, and is constructed from 2x2 cm 2 cells. This gives a total of 10 5 electromagnetic calorimeter 
towers. The GEANT background calculation predicts that each cell sees on average n 7 = 4 background photons per 
crossing with a mean energy E 1 = 1 — 2 McV. If an electromagnetic shower occupies 9 cells, then the mean background 
pedestal will be about 70 MeV. This pedestal can be subtracted from the measured energies. The precision of the 
resulting electron and photon energy measurements will depend on the fluctuations in the mean background energy 
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FIG. 74. Compact tracker geometry in a 4 T field. 

per cell. For an electromagnetic shower occupying 9 cells, the fluctuations in the pedestals are predicted to be about 
10 MeV. This takes into account the fluctuations in the number and the energies of the background photons. 

D. Hadronic calorimeter schemes 

None of the energy generated by background photons in the electromagnetic calorimeter is expected to penetrate 
into the hadronic calorimeter. GEANT calculations show that the total kinetic energy deposited by neutrons in the 
calorimeter is of the order of 140 TeV with an average energy of 30 MeV per neutron. In order to estimate what 
fraction of the kinetic energy of the neutrons will be visible, we should consider the materials involved. For this 
simulation we have presumed an equal mix by volume of liquid argon (as active medium) and copper (as absorber). 
At 30 MeV we expect only a small fraction of the neutrons to knock off protons and only about 10% of the proton 
ionization to be visible in the liquid. Presuming a hadronic calorimeter with 10 4 towers, with the material composition 
described above, the average energy read in the liquid argon will be of the order of 10 MeV per tower with a fluctuation 
of 5 MeV. In summary, a 50 x 50 GeV collider with 4 x 10 12 muons per bunch, the photons and neutrons are expected 
to generate pedestals of 800 and 100 GeV respectively. The estimates for pedestal fluctuations are at or below the 
level of the expected electronic noise. Therefore we believe that the subtraction of these pedestals would present little 
problem both for the electromagnetic and the hadronic calorimeters. The presence of the high neutron background 
should be taken into account in choosing materials for calorimetry. Liquid argon seems a natural choice for the 
electromagnetic calorimeter. 
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The energy deposited by the Bethc-Hcitler muons in the calorimeter is given in table XVII as a function of the 
center of mass energy of the collider. For low center of mass energies, such as the Higgs factory, the Bethc-Hcitler 
muons are not a problem, since there are fewer of them and they leave less energy by catastrophic brcmsstrahlung in 
the calorimeter. For the higher energy option (4 TeV in the CoM or higher), one should explore ways to correct for the 
energy deposition in the calorimeter, such as pattern recognition of the muon tracks or by using timing information. 



E. Muon detector schemes 



The predicted background flux is expected to be relatively modest beyond a radius of 3 m in the vicinity of the 
muon detector. Several possible technologies for muon detectors at a muon collider were discussed during Snowmass 



• Cathode strip chambers. The idea, which is described in the muon collider book [224], is to use MWPCs with 
segmented cathodes and a short (35 ns) drift time to provide prompt signals for triggering. The precision of the 
coordinate measurements would be expected to be of order 50 ^m x a few mm. 

• Threshold Cherenkov counter. The idea is to use a gas Cherenkov radiator to exploit the directionality of 
Chercnkov radiation in order to select high-p^ muons coming from the IP. The device would also give excellent 
timing resolution (of order 2 ns). 



• Long drift jet chamber with pad readout [ 225 (Fig. [75|). Drift time provides the z coordinate, and pad readout 
provides the r-0 coordinates. Directionality at the trigger level is provided by the pattern of pad hits within a 
limited time window. The drift held is provided by cathode strips on grooved G-10 plates. Using 90% argon 
plus 10% CF4 and a maximum drift distance of 50 cm, the maximum drift time is 5 /xs. 



Vertex Originated Muon 




Background Muon 



PC Pads 



Grooves on G-10 / 
Cathode Plates / 



Field Shaping Hire 
and 

Sense wire Alternate 



FIG. 75. Long drift jet chamber with pad readout for muon detection at a muon collider. 



At high energy in the CoM, the channel — > n + fJ.~ + Higgs boson becomes part icula rly attractive to study 
using the muon collider, if the forward going muons from the interaction can be detected [226]. The method provides 
a capability to search for any missing neutral state such as the Higgs boson via the missing mass technique. We are 
investigating methods to improve our forward muon detection capability. 



X. CONCLUSIONS 



Unlike protons, muons are point-like but, unlike electrons, 
therefore can be accelerated and collided in rings. 



they emit relatively little synchrotron radiation and 
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Another advantage resulting from the low synchrotron radiation is the lack of beamstrahlung and the possibility of 
very small collision energy spreads. A beam energy spread of AE/E of 0.003% is considered feasible for a 100 GeV 
machine. It has been shown that by observing spin precession, the absolute energy could be determined to a small 
fraction of this width. These features become important in conjunction with the large s-channel Higgs production 
(/i + /i~ — ► h, 43000 times larger than for e + e~ — ► h), allowing precision measurements of the Higgs mass, width and 
branching ratios. A higher energy muon collider can also distinguish the nearly degenerate heavy Higgs bosons H° 
and A of the minimal supersymmetric extension of the standard model, since these states can also be produced in the 
s channel. We have also examined the ability of the muon collider to study techni-resonances, do a high luminosity 
study of Z boson physics, scan the W and ti thresholds to make precision mass measurements as well as SUSY and 
strongly interacting W boson physics. The high luminosity proton driver and the cold low energy muons permit the 
study of rare kaon and muon decays. Muon storage rings will permit low-systematics studies of neutrino oscillations 
for a wide range of mixing angle and Sm 2 phase space with hitherto unattainable sensitivity. 

Such machines are clearly desirable. The issues are: 

• whether they can be built and physics done with them 

• what they will cost. 

Much progress has been made in addressing the first question and the answer, so far, appears to be yes. It is too 
early to address the second. 

We have studied machines with CoM energies of 0.1, 0.4 and 3 TeV, defined parameters and simulated many of 
their components. Most recent work has been done on the 0.1 TeV First Muon Collider, the energy taken to be 
representative of the actual mass of a Higgs particle. A summary of progress and challenges follows: 

a. Proton driver The specification of the proton driver for the three machines is assumed the same: 10 14 pro- 
tons/pulse at an energy above 16 GeV and 1-2 ns rms bunch lengths. There have been three studies of how to achieve 
these parameters. The most conservative, at 30 GeV, is a generic design. Upgrades of the FNAL (at 16 GeV) and 
BNL (at 24 GeV) accelerators have also been studied. Despite the very short bunch requirement, each study has 
concluded that the specification is attainable. Experiments arc planned to confirm some aspects of these designs. 

b. Pion production and capture Pion production has been taken from the best models available, but an experiment 
(BNL-E910) that has taken data, and is being analyzed, will refine these models. The assumed 20 T capture solenoid 
will require state-of-the-art technology. Capture, decay and phase rotation have been simulated, and have achieved 
the specified production of 0.3 muons per initial proton. The most serious remaining issues for this part of the machine 
are: 

1. The nature and material of the target: The baseline assumption is that a liquid metal jet will be used, but the 
effects of shock heating by the beam, and of the eddy currents induced in the liquid as it enters the solenoid, 
are not yet fully understood. 

2. The maximum rf field in the phase rotation: For the short pulses used, the current assumptions would be 
reasonably conservative under normal operating conditions, but the effects of the massive radiation from the 
nearby target are not known. 

Both these questions can be answered in a target experiment planned to start within the next two years at the BNL 
AGS. 

Polarization of the muon beams represents a significant physics advantage and is an important feature of a muon 
collider. Polarized muon beams are possible. Muons are produced with 100% polarization in the rest frame of the 
pion, but they travel in all directions. By accepting the forward going muons, it is easy to obtain 25% polarization in 
cither beam easily. The amount of polarization can be increased with an accompanying price in luminosity. 

c. Cooling The required ionization cooling is the most difficult and least understood element in any of the muon 
colliders studied. Ionization cooling is a phenomenon that occurs whenever there is energy loss in a strong focusing 
environment. 

But achieving the nearly 10 6 reduction required is a challenge. Cooling over a wide range has been simulated using 
lithium lenses and ideal (linear matrix) matching and acceleration. Examples of limited sections of solenoid lattices 
with realistic accelerating fields have now been simulated, but the specification and simulation of a complete system 
has not yet been done. Much theoretical work remains: space charge and wakeficlds must be included; lattices at the 
start and end of the cooling sequences must be designed; lattices including liquid lithium lenses must be studied, and 
the sections must be matched together and simulated as a full sequence. The tools for this work are nearly ready, and 
this project should be completed within two years. 

Technically, one of the most challenging aspects of the cooling system appears to be: 
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• High gradient rf (e.g. 36 MV/m at 805 MHz) operating in strong (5-10 T) magnetic field, with beryllium foils 
between the cavities. 

An experiment is planned that will test such a cavity, in the required fields, in about two years time. On an 
approximately six year time scale, a Cooling Test Facility is being proposed that could test ten meter lengths of 
different cooling systems. If they are required, then an urgent need is to develop: 

• Lithium Lenses: (e.g. 2 cm diameter, 70 cm long, liquid lithium lenses with 10 T surface fields and a repetition 
rate of 15 Hz). 

The use of 31 T solenoids could avoid their need, at least in the low energy First Muon Collider, which would ease the 
urgency of this rather long term R&D, but both options would require long-term R&D. Meanwhile a short lithium 
lens is under construction at BINP (Novosibirsk, Russia). 

d. Acceleration The acceleration system is probably the least controversial, although possibly the most expensive, 
part of a muon collider. Preliminary parameters have been specified for acceleration sequences for a 100 GcV and 
a 3 TeV machine, but they need refinement. In the low energy case, a linac is followed by three recirculating or 
FFAG accelerators. In the high energy accelerator, the recirculating or FFAG accelerators are followed by three fast 
ramping synchrotrons employing alternating pulsed and superconducting magnets. The parameters do not appear to 
be extreme, and it does not appear as if serious problems are likely. 

e. Collider The collider lattices are challenging because of the requirement of very low beta functions at the 
interaction point, high single bunch intensities, and short bunch lengths. However, the fact that all muons will decay 
after about 800 turns means that slowly developing instabilities are not a problem. Feasibility lattices have been 
generated for a 4 TeV case, and more detailed designs for 100 GeV machines are been studied. In the latter case, but 
still without errors, 5er acceptances in both transverse and longitudinal phase space have been achieved in tracking 
studies. Beam scraping schemes have been designed for both the low energy (collimators) and high energy (septum 
extractors) cases. 

The short bunch length and longitudinal stability problems are avoided if the rings, as specified, are sufficiently 
isochronous, but some rf is needed to remove the impedance generated momentum spread. Transverse instabilities 
(beam breakup) should be controlled by rf BNS damping. 

The heating of collider ring superconducting magnets by electrons from muon decay can be controlled by thick 
tungsten shields, and this technique also shields the space surrounding the magnets from the induced radioactivity 
on the inside of the shield wall. A conceptual design of magnets for the low energy machine has been defined. 

Although much work is yet to be done (inclusion of errors, higher order correction, magnet design, rf design, etc), 
the collider ring does not appear likely to present a serious problem. 

/. Neutrino radiation and detector background Neutrino radiation, which rises as the cube of the energy, is 
not serious for machines with center of mass energies below about 1.5 TeV. It is thus not significant for the First 
Muon Collider; but above 2 TeV, it sets a constraint on the muon current and makes it harder to achieve desired 
luminosities. However, advances in cooling and correction of tunc shifts may still allow a machine at 10 TeV with 
substantial luminosity (> 10 35 cm _2 s _1 ). 

Background in the detector was at first expected to be a very serious problem, but after much work, shielding 
systems have evolved that limit most charged hadron, electron, gamma and neutron backgrounds to levels that are 
acceptable. Muon background, in the higher energy machines, is a special problem that can cause serious fluctuations 
in calorimeter measurements. It has been shown that fast timing and segmentation can help suppress this background, 
and preliminary studies of its effects on a physics experiment are encouraging. The studies are ongoing. 

g. Detector scenarios We have considered several options for the experimental detector components for various 
CoM energy colliders. Much work needs to be done to optimize the physics reach at each energy by feeding back the 
results of detailed simulations of backgrounds and signal to the detector design. Only then will the feasibility of doing 
physics with a muon collider be fully explored. 

h. Summary Much progress has been made since Snowmass, but much still needs to be done. A time scale of 
two years should allow completion of simulation studies and the experimental testing of crucial technical challenges. 
Prototype construction and testing will take another 4-6 years. The construction of a First Muon Collider by about 
2010 docs seems possible. 
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